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1  Summary 


A  viscous/in  viscid  computational  method  for  the  design  and  analysis  of  quasi-3D  cascades 
has  been  developed.  The  specific  application  targeted  is  supersonic  compressor  blading  with 
splitter  blades,  although  the  method  is  applicable  to  any  type  of  cascade.  The  method  is 
an  extension  of  the  ISES  viscous /inviscid  methodology.  A  streamline-based  inviscid  Eu¬ 
ler  formulation  is  fully  coupled  to  an  integral  boundary  layer  formulation  to  describe  the 
overall  viscous  fiowfield.  Rotation  and  streamtube  contraction  effects  have  been  incorpo¬ 
rated.  The  method  is  accurate  with  regards  to  loading  and  loss,  and  is  substantially  faster 
than  equivalent  Navier-Stokes  solvers.  Inverse  design  and  optimization  capabilities  are  also 
implemented,  giving  an  effective  design/analysis  system. 

Existing  supersonic  splittered  cascades  were  investigated  to  determine  what  character¬ 
istics  are  desirable  or  undesirable  in  this  type  of  blading.  Sensitivity  studies  indicate  that 
tangential  splitter  position  and  splitter  loading  details  have  little  influence  on  performance. 
Substantial  improvements  are  possible  if  the  splitter  is  moved  aft  into  a  tandem-blade 
arrangement  with  the  main  blade.  Comparisons  between  computational  results  and  mea¬ 
sured  data  strongly  suggest  that  traditional  supersonic  cascade  tests  involve  substantial 
three-dimensional  effects  which  are  not  correctable  by  measured  streamtube  contraction. 
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2  Completed  Work 


All  the  goals  of  the  original  proposal  have  been  completed.  The  work  is  documented  in  the 
S.M.  thesis  of  Harold  H.  Youngren,  which  is  included  in  the  Appendix.  The  specific  work 
items  which  were  originally  proposed  are  elaborated  below. 

2.1  Extension  of  method  to  curvilinear  coordinates 

The  existing  ISES  method  has  been  converted  from  a  cartesian  x-y  coordinate  system  to 
an  m'-O  system.  This  is  necessary  for  rotation,  radius  changes,  and  streamtube  thickness 
effects  to  be  introduced.  The  streamtube  thickness  effects  were  extensively  validated  with 
experimental  data.  The  radius  change  and  rotation  effects  were  not  validated  due  to  lack 
of  suitable  data.  However,  checks  were  performed  on  simple  cases  for  which  quantities  such 
as  total  pressure  rise  could  be  analytically  determined. 

2.2  Modification  of  boundary  layer  formulation 

The  existing  integral  method  of  the  ISES  code  has  been  modified  to  account  for  the  rotation, 
radius  change,  and  streamtube  thickness  effects.  In  the  end,  the  modifications  have  turned 
out  to  be  almost  trivial,  with  a  single  term  accounting  for  streamtube  contraction,  and 
several  very  small  terms  accounting  for  rotation  and  radius  changes.  Validation  of  the 
latter  terms  is  difficult  because  of  their  smallness,  and  was  not  performed  due  to  the  lack  of 
suitable  data.  The  more  important  streamtube  contraction  term  was  validated  in  numerous 
comparisons  with  experiments. 

2.3  Investigation  of  effective  cascade  design  procedures 

The  inverse  capabilities  of  the  original  ISES  method  were  adapted  to  turbomachinery  blad¬ 
ing.  The  standard  inverse  problem  where  surface  pressures  are  prescribed  is  ill-posed  if  any 
separation  is  present.  An  alternative  approach  in  which  the  geometry  changes  are  restricted 
to  specified  geometric  modes  was  implemented  as  a  well-posed  alternative.  These  modes 
have  been  generalized  to  permit  the  movement  and/or  rotation  of  the  entire  splitter  blade. 
Also,  optimization  capabilities  readily  fall  out  of  this  formulation,  permitting  loss  mini¬ 
mization,  for  example.  The  overall  method  can  be  viewed  as  a  “toolbox”  which  permits  the 
designer  to  manipulate  the  geometric  and/or  aerodynamic  characteristics  of  the  cascade. 
The  method  is  fast  enough  to  permit  interactive  execution  on  modern  desktop  workstations. 

2.4  Evaluation  of  the  splitter  blade  concept 

The  computational  studies  documented  in  the  Appendix  indicate  that  the  splitter  configu¬ 
ration  initially  targeted  in  the  proposal  has  marginal  benefits  compared  to  a  conventional 
configuration.  The  splitter  increases  the  turning,  but  at  a  significant  loss  penalty.  The 
tangential  position  of  the  splitter  and  the  details  of  the  pressure  distribution  on  the  splitter 
have  only  minor  effects  on  the  performance  of  the  overall  configuration.  Very  significant 
benefits  result,  however,  if  the  splitter  is  moved  aft  and  just  under  the  trailing  edge  of  the 
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main  blade.  This  tandem  arrangement  gives  large  increases  in  turning  as  well  as  loss  re¬ 
duction.  Geometric  constraints  may  prohibit  this  configuration,  depending  on  the  specific 
application.  The  computational  method  would  serve  very  well  as  a  means  to  efficiently 
evaluate  such  novel  blading  concepts  for  any  new  turbomachine  blading  design. 
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A  new  computational  method,  MISES,  is  developed  for  turbomachinery  design  and 
analysis  applications.  The  method  is  based  on  the  fully  coupled  viscous/inviscid  method, 
ISES,  and  is  applicable  to  blade-to-blade  analysis  of  axial  fan  and  compressor  stator  or 
rotors  with  optional  splitter  vanes.  Quasi-three  dimensional  effects  for  stream  surface 
radius,  streamtube  thickness  and  wheel  rotation  may  be  included.  The  flow  is  modeled 
with  the  steady  Euler  equations  and  the  integral  boundary  layer  equations.  A  robust 
Newton-Raphson  method  is  used  to  solve  the  coupled  non-linear  system  of  equations, 
requiring  only  several  minutes  for  solution  on  a  typical  workstation.  Design  options  are 
implemented  for  either  single  surface  or  camber  redesign. 

The  method  is  exercised  by  comparison  with  transonic  cascade  tests  to  validate  the 
quasi-three-dimensional  formulation.  The  results  show  excellent  correlation  to  measured 
pressure  distributions  and  loss  levels.  The  multiple  blade  capability  is  demonstrated  by 
comparison  to  test  data  for  a  supersonic  cascade  with  splitter  vane.  New  splitter  vane 
configurations  for  improving  the  performance  of  the  supersonic  cascade  are  explored, 
resulting  in  large  increases  in  turning  and  reduced  loss. 
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Chapter  1 

Introduction 


This  thesis  will  present  a  new  computational  method,  MISES1,  which  extends  the  fully 
coupled  viscous/inviscid  method,  ISES,  for  application  to  quasi- three- dimensional  cas¬ 
cade  design.  This  method  is  intended  for  use  as  a  tool  for  compressor  cascade  design 
and  analysis  and  for  blade- to- blade  design  of  axial  stators  and  rotors. 

The  viscous  flow  is  modeled  using  the  steady  Euler  equations  and  integral  boundary 
layer  equations  solved  with  a  Newton- Raphson  technique,  the  approach  used  in  the 
ISES  code  of  Giles  and  Drela  [lj.  This  has  several  advantages  for  design  applications  — 
accuracy,  speed  and  a  robust  inverse  design  capability.  Primarily,  it  is  the  speed  that 
makes  this  approach  attractive  for  design  —  a  typical  case  is  solved  in  3-10  Newton 
cycles,  a  matter  of  minutes  on  a  fast  workstation.  The  accuracy  also  has  a  strong 
influence  on  the  design  process  —  the  drag  predictions  from  ISES  are  sufficiently  good 
that  they  can  be  reliably  used  by  a  designer  to  optimize  an  airfoil  at  both  on-design 
and  off-design  conditions,  using  wind  tunnel  testing  only  to  verify  the  computational 
predictions.  Although  blade-to-blade  design  for  turbomachinery  is  similar  to  the  isolated 
airfoil  design  problem,  there  are  additional  difficulties  in  solving  the  viscous/inviscid 
governing  equations,  which  are  strongly  coupled  for  transonic  cascades.  A  design  method 
for  turbomachinery  requires  accurate  loss  prediction  for  strong  shocks  and  for  separated 
flows. 

This  work  also  focuses  on  the  analysis  and  design  of  cascades  with  multiple  blades, 
specifically  splitter  vanes.  These  are  introduced  to  reduce  loss  and  increase  turning 
by  changing  the  distribution  of  blade  loading.  They  also  increase  blockage  and  can 
increase  viscous  losses.  This  thesis  addresses  some  of  the  design  issues  for  multiple 
'To  answer  the  inevitable  question  -  Multiple  blade  Interacting  Streamtube  Euler  Solver. 
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blade  cascades,  focusing  on  the  optimization  of  the  splitter  vane  in  a  highly  loaded 
supersonic  cascade. 


1.1  Motivation  and  Background 

Improved  levels  of  performance  from  gas  turbine  engines  have  motivated  designers  to 
place  increasing  reliance  on  computational  techniques  in  the  design  process.  This  is 
especially  true  as  blade  loadings  in  compressors  and  turbines  are  increased  and  higher 
levels  of  refinement  in  design  are  attained.  Unlike  typical  external  aerodynamic  design 
problems,  where  the  design  is  driven  to  virtually  eliminate  strong  non-linear  flow  effects, 
the  interacting  compressible  and  viscous  effects  in  cascades  lead  to  complex  flows  that 
frequently  involve  significant  three-dimensionality,  flow  separation  and  shock  loss. 

The  current  state  of  the  art  in  computational  methods  for  cascade  analysis  is  three- 
dimensional  Navier-Stokes  simulations,  normally  using  the  Reynolds-averaged  form  of 
the  equations  with  a  turbulence  model.  These  methods  have  the  advantage  of  modeling 
all  of  the  relevant  viscous  blockage  effects  and  may  be  applied  to  analyze  virtually  any 
cascade  flow.  Navier-Stokes  methods  are  presently  too  expensive  and  cumbersome  to 
use  for  the  bulk  of  the  design  process,  typically  requiring  300  to  6000  time  steps  to 
approximate  steady  state,  involving  hours  of  time  on  fast  computers. 

Although  the  flow  in  axial  turbomachinery  is  three-dimensional,  a  useful  and  often 
necessary  simplification  for  design  purposes  is  to  approximate  the  flow  through  a  stage 
as  a  set  of  two-dimensional  blade-to-blade  problems  defined  on  axisymmetric  stream 
surfaces.  Axisymmetric  through-flow  codes  are  used  early  in  the  preliminary  design 
process  to  define  circumferentially  averaged  conditions  in  one  or  more  stages  of  the 
machine  based  on  initial  estimates  of  work  and  loss.  These  calculations  define  the 
flow,  in  terms  of  axisymmetric  stream  surfaces.  At  the  next  level  of  design  refinement 
the  stream  surface  radius  and  spacing  can  be  used  to  define  a  set  of  quasi-3D  blade- 
to-blade  design  problems  for  each  stage.  These  allow  the  designer  to  select  or  design 
blade  profiles  at  several  radial  stations  to  define  the  complete  three-dimensional  rotor 
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or  stator  blade.  The  blade-to-blade  technique  works  surprisingly  well  for  most  design 
applications,  limited  in  effectiveness  largely  by  the  estimates  for  boundary  layer  effects 
on  the  inner  and  outer  casing  walls  and  by  three-dimensional  effects  not  accounted  for 
with  the  axisymmetric  assumptions.  In  the  context  of  simpler,  faster  and  reasonably 
accurate  preliminary  design  tools,  the  more  complex  three-dimensional  methods  can  be 
most  effectively  used  as  ‘numerical  wind  tunnels’  to  verify  the  design. 


1.2  Cascades  with  Splitter  Vanes 

One  of  the  goals  of  this  work  is  to  develop  capability  for  analysis  and  design  of  cas¬ 
cades  with  intermediate  blades,  or  splitter  vanes,  that  are  added  to  the  basic  blading  to 
increase  turning  or  reduce  flow  separation  and  loss.  The  case  used  as  a  guide  for  this 
development  was  a  high  turning  supersonic  compressor  stage  with  splitter  vane  devel¬ 
oped  and  tested  by  the  Aerospace  Research  Lab  (ARL)  in  the  early  to  mid  1970’s.  The 
high  inlet  Mach  number  and  strong  sidewall  contraction  (2:1)  made  this  an  extremely 
challenging  test  case.  The  blade  cross-sections  for  the  ARL  two-dimensional  cascade 
are  shown  in  Figure  1.1.  The  splitter  vane  for  this  case  was  intended  to  operate  in  the 
subsonic  portion  of  the  passage  downstream  of  the  shock.  The  basic  effect  of  the  splitter 
vane  on  blade  loading  is  shown  in  Figure  1.2,  where  the  main  blade  loading  is  relieved 
by  the  overlapping  vane. 

The  ARL  cascade  has  been  the  subject  of  several  previous  studies,  both  experimental 
and  analytical.  This  cascade  started  out  as  part  of  a  high  pressure-ratio,  high  diffusion 
axial  rotor  developed  for  the  USAF.  Tests  had  indicated  excessive  deviation  angles  that 
were  reduced  by  the  addition  of  a  partial  flow  splitter.  The  two-dimensional  cascade 
was  tested  by  Holtman  et  al  as  a  single  blade  [2j  and  later  with  a  splitter  vane  [3].  A 
rotor  using  this  concept  was  later  designed  and  tested  by  Wennerstrom  and  Frost  [4], 
encountering  test  problems  from  stator  choking.  A  later  experimental  study  [5]  in  1977 
examined  optimization  of  splitter  vane  location. 

An  analytical  study  by  Dodge  [6]  in  1975,  using  a  2D  inviscid  potential  method, 
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Figure  1.1:  Blade  sections  for  ARL  cascade  with  splitter  vane. 

examined  the  ARL  cascade  with  and  without  the  splitter  vane,  attempting  to  do  a 
design  study  with  limited  success.  A  second  effort  [7]  in  1978,  again  by  Dodge,  brought  a 
much  more  capable  inviscid  3D  potential  method  to  bear  on  the  ARL  cascade  and  rotor, 
although  the  lack  of  viscous  blockage  effects  made  comparisons  to  test  data  somewhat 
inconclusive. 

The  large  loss  levels  and  deviation  angles  from  the  experiment  indicate  that  the  ARL 
cascade  has  extensive  separated  flow,  which  is  likely  the  source  of  problems  in  earlier 
analytical  work.  The  current  method,  with  its  full  inviscid/ viscous  coupling  scheme,  is 
capable  of  modeling  these  blockage  effects  and  should  give  a  more  complete  picture  of 
the  splitter  vane  interaction  with  the  main  blade. 
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Figure  1.2:  Effect  of  splitter  vane  on  pressure  distributions  for  ARL  cascade. 

1.3  Overview  of  Thesis 

The  thesis  is  structured  so  the  introductory  material  briefly  discusses  the  development 
of  the  method,  while  the  bulky  derivations  and  detailed  material  is  presented  in  the 
appendices.  A  brief  sketch  of  the  contents  of  the  thesis  is  given  below: 

Chapter  2  introduces  the  coordinate  systems,  governing  equations  and  boundary  con¬ 
ditions  appropriate  to  the  quasi-3D  blade  to  blade  problem. 

Chapter  3  presents  the  discrete  form  for  the  equations  and  boundary  conditions  on 
the  intrinsic  streamline  grid.  The  Newton- Raphson  method  used  to  solve  the  non¬ 
linear  system  is  discussed.  Design  options  and  the  modifications  for  multiple  blade 
analysis  and  design  are  introduced.  The  offset-periodic  grid  and  solver  used  for 
supersonic  cascades  is  discussed. 


24 


Chapter  4  exercises  the  method  on  several  analysis  test  cases  at  subsonic,  transonic 
and  supersonic  inlet  Mach  numbers  to  validate  the  quasi-3D  capabilities.  Solutions 
for  the  supersonic  ARL  cascade  with  splitter  vane  are  compared  to  test  data. 

Chapter  5  examines  the  ARL  cascade  with  splitter  vane  in  greater  detail.  The  effects 
the  splitter  vane  position  are  analyzed  and  compared  to  test  data.  The  splitter 
and  main  blade  are  redesigned  to  optimize  loss  and  turning.  Additional  avenues 
for  cascade  modification  are  explored  using  the  method. 

Appendix  A  presents  the  derivation  of  the  discrete  form  of  the  quasi-3D  Euler  equa¬ 
tions  on  the  streamline  grid. 

Appendix  B  outlines  the  modifications  to  the  integral  boundary  layer  method  to  in¬ 
clude  quasi-3D  effects. 

Appendix  C  presents  the  loss  model  and  the  loss  sensitivity  calculation  used  in  design 
optimization. 

Appendix  D  contains  the  details  of  the  Jacobian  matrices  and  matrix  solution. 
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Chapter  2 

Theory 

2.1  Cascade  Coordinate  Systems 

The  geometric  quantities  defining  a  two-dimensional  cascade  are  illustrated  in  Fig¬ 
ure  2.1.  The  cascade  blades  are  periodic  in  y  with  pitch  ypu,h •  The  blade  chordlines 


Figure  2.1:  Cascade  coordinate  system 

are  inclined  to  x  direction  by  the  stagger  angle  The  inlet  flow  direction  is  d\  and  the 
outlet  flow  angle  is  0j. 

The  three-dimensional  flow  in  an  axial  stator  or  rotor  is  treated  as  an  axisymmetric 
meridional  flow  and  a  quasi-30  blade-to-blade  flow  defined  on  the  axisymmetric  stream 
surfaces,  as  shown  in  Figure  2.2.  In  this  figure  blade  elements  for  the  blade-to-blade  flow 
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are  shown  defined  by  a  stream  surface  of  revolution.  The  axisymmetric  stream  surfaces 
are  defined  in  terms  of  the  radius  r(z)  and  streaxntube  thickness  b(z),  where  z  is  the 
axial  coordinate,  and  are  shown  as  streamlines  in  the  meridional  plane  in  Figure  2.3. 
Meridional  streamlines  can  be  determined  using  a  through-flow  method  to  solve  the 


Figure  2.2:  Blade-to-blade  flow  on  a  stream  surface  of  revolution 


Figure  2.3:  Meridional  (r-z)  plane  for  defining  stream  surfaces 

axial  and  radial  equilibrium  relations  for  mass,  momentum  and  energy  given  the  basic 
stage  inflow  and  outflow  parameters  and  the  hub  and  tip  casing  shapes.  The  flow  over 
a  complete  blade  can  be  approximated  by  a  set  of  quasi-3D  blade-to-blade  problems 
for  stream  surfaces  at  several  radial  stations,  a  technique  that  forms  the  basis  of  many 
preliminary  design  systems  for  compressor  design. 
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The  cascade  definitions  in  Figure  2.1  are  also  used  for  axial  stator  blades  with  x  and  y 
corresponding  to  the  axial  coordinate  z  and  circumferential  coordinate  6.  Alternatively, 
t  and  y  will  correspond  to  a  local  m\  9  coordinate  system  defined  on  the  stream  surface. 
This  choice  is  advantageous  because  it  is  a  conformal  mapping  of  the  surface  of  revolution 
to  a  plane,  preserving  angles  and  shapes.  The  new  coordinate  m!  is  defined  as  the 
normalized  arc- length  in  the  r  -  z  plane. 


In  the  current  method  the  quasi-30  blade-to-blade  flow  is  solved  on  a  stream  surface  in 
an  m',  9  coordinate  system  with  specified  radius  and  stream  tube  thickness. 

These  definitions  also  apply  for  a  rotor,  except  that  the  flow  problem  is  defined  in 
a  blade-relative  coordinate  system  rotating  with  respect  to  an  inertial  frame.  Referring 
to  Figure  2.1,  the  blades  on  an  axial  rotor  are  defined  to  rotate  with  singular  velocity  fi 
in  the  -y  direction. 

2.2  Viscous/Inviscid  Flow  Solution 

Unlike  most  airfoil  flows,  cascade  flows  normally  exhibit  strong  coupling  between  the 
inviscid  flow  and  the  viscous  boundary  layers,  particularly  at  transonic  and  supersonic 
Mach  numbers.  This  coupling  can  be  sufficiently  strong,  involving  extensive  flow  separa¬ 
tion,  that  the  flow  cannot  be  solved  without  considering  viscous  effects  on  blockage,  e.g. 
the  ARL  cascade  in  Chapter  4.  Viscous  effects  must  also  be  included  for  the  prediction 
of  cascade  losses. 

Instead  of  solving  the  viscous  flow  directly,  the  zonal  approach  of  ISES  is  used.  An 
equivalent  inviscid  flow,  EEF  (8),  is  postulated  using  a  displacement  surface  to  represent 
the  viscous  layer,  as  illustrated  in  Figure  2.4.  The  inner  boundary  of  the  EIF  is  displaced 
outward  from  the  wall  by  the  boundary  layer  displacement  thickness  S'.  The  outer 
inviscid  flow  is  solved  using  the  steady  state  Euler  equations  while  the  boundary  layer 
flow  is  solved  using  an  integral  boundary  layer  method.  The  outer  flow  is  coupled 


Figure  2.4:  Equivalent  inviscid  flow  defined  by  BL  displacement  thickness 

to  the  boundary  layer  through  the  edge  velocity  and  density  that  drive  the  integral 
boundary  layer  solution  for  the  displacement  thickness.  These  two  interdependent  flow 
domains  are  fully  coupled  by  a  simultaneous  Newton-Raphson  solution  procedure  for 
the  non-linear  equations. 


2.3  Steady-State  Euler  Equations 

The  inviscid  flow  is  solved  using  the  steady-state  Euler  equations  in  a  blade-relative 
coordinate  system.  The  integral  form  for  the  Euler  equations  is  derived  using  the 
conservation  laws  for  mass,  momentum  and  energy  applied  to  a  control  volume  V  with 
boundary  dV  and  unit  normal  n. 

Mm 

Momentum 

Energy 

These  equations  apply  to  a  reference  frame  rotating  with  angular  velocity  H  with  respect 
to  the  inertial  frame,  adding  rotational  terms  to  the  momentum  and  energy  equations. 


&  p( q  •  n)d4  =  0 

(2.2) 
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The  source  term  in  the  momentum  equation,  f  represents  the  Coriolis  and  centrifugal 
force  terms  in  the  blade  relative  system. 

f  =  ft  X  (il  x  f)  +  2fi  x  q  (2.5) 


The  conserved  quantity  in  energy  equation  is  the  rothalpy ,  I,  which  is  invariant  in  the 
rotating  system. 


I 


(nrf 

2 


(2.6) 


The  inviscid  formulation  is  completed  with  the  assumption  of  the  ideal  gas  law. 


h  = 


7  P 
7-1  P 


with  7  =  cp/c„ 


(2.7) 


The  strategy  for  solving  the  inviscid  equations  is  to  discretize  them  on  a  streamline 
grid,  eliminating  convection  terms  across  the  streamline  faces.  This  leads  to  a  particu¬ 
larly  simple  form  for  the  Euler  equations,  where  the  continuity  and  energy  equations  are 
replaced  by  statements  of  constant  mass  flow  and  rothalpy  in  each  streamtube,  leaving 
only  the  equations  for  streamwise  and  normal  momentum  to  be  solved. 


2.4  Boundary  Layer  Solution 


The  viscous  flow  is  solved  using  the  integral  boundary  layer  method  from  ISES  [9], 
modified  to  include  quasi-3D  effects,  see  Appendix  B.  This  is  a  two-equation  method 
employing  the  venerable  von  Kdrman  integral  momentum  equation  and  the  kinetic 
energy  shape  parameter  equation.  In  turbulent  flow  regions  a  dissipation  lag  equation, 
similar  to  Green’s  lag  equation  [10],  is  added  to  model  upstream  history  effects  on  the 
Reynolds  stress.  These  equations  have  the  form 

^i(M',ur)  (2.8) 

F2(6,6\uK,Cr)  (2.9) 

JFj  (M*,u„Cr)  (2.10) 


d9_ 

da 

d[H-(9,6',uK)  ] 
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Relations  derived  using  Swafford’s  turbulent  flow  profiles  provide  turbulent  closure  for 
the  shape  factors,  skin  friction  and  dissipation. 

In  laminar  flow  regions  the  same  basic  scheme  is  used,  with  the  dissipation  lag 
equation  (2.10)  replaced  by  a  transition  equation.  An  envelope  method  for  maximum 
disturbance  amplification  is  used,  based  on  spatial  amplification  of  Tollmein-Schlichting 
waves  in  the  Orr-Sommerfeld  equation.  The  maximum  disturbance  amplitude  variable 
N  replaces  Cr  in  the  BL  equation  set.  The  amplification  equation  has  the  form 

dN 

—  =  (2.11) 

The  transition  point  is  determined  by  the  variable  N  exceeding  a  specified  critical 
value  Nr The  result  is  roughly  equivalent  to  the  classical  el)  transition  prediction 
method.  The  laminar  flow  equations  are  completed  using  closure  relations  derived  from 
correlations  with  Falkner-Skan  profiles. 

The  boundary  layer  solution  is  influenced  by  the  edge  velocity  u,. ,  which  is  also 
the  surface  velocity  in  the  EIF.  Correspondingly,  the  displacement  thickness  from  the 
boundary  layer  influence'  .he  inviscid  flow  by  offsetting  the  edge  of  the  inviscid  region 
from  the  surface.  Because  the  inviscid  and  viscous  solutions  are  fully  coupled  in  the 
Newton  system,  flows  involving  separation  may  be  calculated. 

2-5  Boundary  Conditions 

The  principal  requirement  for  the  boundary  conditions  is  that  they  must  completely 
specify  the  flow  problem  and  be  physically  realizable.  This  is  a  more  serious  problem 
for  cascades  than  for  airfoils  due  to  the  finite  domain  and  transonic  or  supersonic  con¬ 
ditions  into  or  within  the  passage,  presenting  far  more  opportunities  to  specify  ill-posed 
boundary  conditions. 

The  boundary  conditions  must  specify  the  inflow  thermodynamic  quantities  —  stag¬ 
nation  density  and  stagnation  enthalpy.  Inflow  direction  must  also  be  specified,  although 
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Mnch  number 


Figure  2.5:  Variation  of  mass  flow  with  Mach  number 

the  form  this  takes  depends  on  the  Mach  number.  One  additional  quantity  must  be  spec¬ 
ified,  either  the  mass  flow  or  the  exit  pressure.  Additional  boundary  conditions  that 
constrain  added  variables  (or  DOF’s)  are  required  on  a  one-to-one  basis.  These  are 
often  used  to  drive  the  system  to  a  specified  flow  condition  or  for  design  applications. 

Some  combinations  of  boundary  conditions  are  ill-posed  -  consider  the  mass  flow  in 
a  passage  of  area  A,  shown  as  a  function  of  Mach  number  in  Figure  2.5.  Up  to  the  point 
at  which  the  passage  chokes,  the  specification  of  mass  flow  is  well-posed.  Beyond  the 
choke  point,  further  increases  of  mass  flow  are  ill-posed.  In  fact,  the  discrete  equations 
with  specified  mass  flow  become  increasingly  ill-conditioned  as  the  choking  point  is 
approached,  leading  to  numerical  problems.  These  can  be  alleviated  by  specifying  exit 
pressure  and  allowing  the  mass  flow  be  determined  by  the  solution,  which  is  well-posed 
because  mass  flow  is  uniquely  specified  by  exit  pressure,  while  the  reverse  is  not  true. 

The  Newton  method  used  to  solve  the  non-linear  equations  requires  an  initial  flow 
condition  from  which  to  begin  iterating  to  a  solution.  For  a  given  mass  flow  in  Figure  2.5 
there  exists  two  possible  Mach  numbers,  one  subsonic  and  one  supersonic.  It  is  essential 
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that  the  initial  state  correspond  to  the  desired  branch. 


2.5.1  Subsonic  BC’s 

For  subsonic  inflow  the  boundary  conditions  that  must  be  specified  depend  on  whether 
the  flow  in  the  passage  is  choked.  These  take  the  form 


Subsonic  BC 


inflow  BC 

outflow  BC 

determined 

Pi h 

Peril 

or 

Pih  ),01 

Peril 

m 

Either  of  these  combinations  may  be  used  for  unchoked  flow.  Exit  pressure  need  not 
be  specified  since  pexn  is  uniquely  determined  for  unchoked  flow  by  the  inlet  conditions 
and  cascade  geometry.  Alternatively,  pext,  can  be  specified,  provided  one  of  the  inflow 
boundary  conditions  is  dropped  to  avoid  over-specifying  the  problem.  The  inlet  flow 
angle  i  must  be  specified,  although  it  should  be  imposed  sufficiently  far  forward  of  the 
cascade  that  this  specification  is  relatively  unaffected  by  its  flow  disturbance.  The  inlet 
mass  flow  specifies  the  inlet  Mach  number,  given  the  other  inlet  quantities. 

As  the  flow  becomes  choked  there  is  no  longer  a  unique  exit  pressure  for  a  given  set 
of  inflow  conditions.  The  second  boundary  condition  must  be  used,  with  pr:T,i  specified 
and  one  of  the  inflow  boundary  conditions  dropped,  usually  mass  flow  since  it  is  always 
uniquely  determined  by  exit  pressure. 

2.5.2  Supersonic  BC’s 

As  a  result  of  the  high  stagger  angles  used  on  supersonic  cascades  the  inflow  is  normally 
also  axially  subsonic,  i.e.  M\  cos0\  <  1.0.  The  incoming  flow  is  influenced  by  the 
fixed  pattern  of  expansion  and/or  compression  waves  that  run  axially  forward  from  the 
supersonic  upper  surface  of  the  blade,  see  Figure  2.6.  These  waves  turn  the  incoming 
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inlet 


Figure  2.6:  Supersonic/axially-subsonic  inflow  for  cascade,  showing  inlet  wave  system 

flow  to  approximately  align  with  the  supersonic  portion  of  the  blade  upper  surface.  At 
the  point  where  the  cascade  is  choked  the  upstream  flow  into  the  infinite  cascade  is  at 
the  so-called  “unique  incidence”  condition  for  the  cascade.  In  practice,  whenever  the 
inflow  is  supersonic  the  range  of  inlet  flow  angles  into  the  cascade  is  quite  close  to  the 
unique  incidence  point. 
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The  inflow  angle  cannot  be  arbitrarily  specified  for  supersonic  flow  due  to  the  in¬ 
fluence  of  the  cascade  wave  system.  The  correct  inlet  condition  for  supersonic/axially- 
subsonic  inflow  is  to  specify  an  inlet  characteristic  condition  so  that  the  inlet  boundary 
acts  as  a  free  boundary  to  expansion  or  compression  waves  in  the  flow  domain,  see  [11, 
page  291].  This  inflow  condition  is  defined  by 

cm  =  9  ±  i/(M)  =  cinitpee  (2.12) 
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where  9  is  the  local  inlet  flow  angle  and  v{M)  is  the  Prandtl-Meyer  function, 


v  =  — -arctan  \  (M2  -  1)  - — \  -  arctan  V M2  -  1  (2.13) 

7+1  V  7+1 

With  inlet  characteristic  conditions  the  inflow  angle  and  Mach  number  are  no  longer 
independent.  As  the  inflow  angle  is  changed  by  expansion  or  compression  waves  from 
the  cascade,  the  local  Mach  number  changes  along  the  inlet  characteristic.  The  Mach 
and  local  flow  angle  dependence  for  a  family  of  characteristics  is  shown  in  Figure  2.7. 
Since  the  slope  d9/dM  vanishes  at  the  sonic  point,  the  Mach  number  becomes  very 
sensitive  to  flow  angle  changes  near  M  =  1.0. 


Figure  2.7:  Characteristic  families  for  local  angle  9  and  Mach 


2.5.3  Periodic  BC 

Periodic  conditions  for  geometric  continuity  must  be  imposed  on  the  stagnation  stream¬ 
lines  upstream  and  downstream  of  the  blade,  in  addition  the  streamlines  must  be 
force-free,  set  by  imposing  equal  pressures  across  the  periodic  boundary. 
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2.5.4  Solid  Wall  BC 


Since  a  streamline-based  solution  technique  is  used  for  the  inviscid  flow,  solid  wall 
conditions  are  easily  imposed  since  no  mass  flow  crosses  the  blade  surface  streamline, 
by  definition.  The  solid  wall  boundary  condition  fixes  the  streamline  position  to  the 
blade  surface,  i.e.  Sn  =  0,  where  n  is  the  coordinate  normal  to  the  surface  streamline. 
For  viscous  cases,  where  the  EIF  is  displaced  by  the  boundary  layer,  this  becomes 
Sn  =  S  (£"),  where  S'  is  the  displacement  thickness  from  the  BL  solution. 

2.5.5  Design  BC 

One  of  the  useful  options  for  design  applications  is  a  form  of  inverse  boundary  condition 
where,  instead  of  specifying  the  surface  position  and  having  the  pressure  set  by  the  flow 
field,  the  surface  pressure  is  specified  directly  and  the  streamline  geometry  is  determined 
by  the  flow  solution.  This  is  easily  done  within  the  framework  of  the  streamline  solver 
by  imposing  a  wall  pressure  condition  in  place  of  the  solid  wall  condition.  An  arbitrary 
specified  pressure  distribution  may  not  give  a  physically  realizable  blade,  so  specified 
pressure  distributions  are  subject  to  several  compatibility  conditions,  as  discussed  in  [9]. 

Inverse  design  boundary  conditions  cannot  be  applied  over  a  separated  region  be¬ 
cause  there  exists  only  a  very  weak  local  constraint  for  surface  displacement  in  the 
separated  BL.  A  modal  design  scheme  has  been  implemented,  using  a  set  of  geometric 
shape  functions  (or  bump  functions),  to  circumvent  this  limitation.  Modal  design  is 
discussed  in  greater  detail  in  Chapter  3.  Camber  changes  to  achieve  a  specified  pres¬ 
sure  loading  for  fixed  thickness  may  also  be  determined  using  a  variant  of  the  periodic 
boundary  condition. 
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2.6  Loss  Calculation 


Losses  through  a  cascade  result  from  the  entropy  rise  from  inlet  to  outlet  produced  by 
inviscid  shock  losses  and  viscous  dissipation  in  the  blade  boundary  layers.  The  increase 
in  entropy  causes  a  reduction  of  the  stagnation  pressure  from  the  equivalent  isentropic 
value.  The  loss  coefficient  w  is  defined  as  the  reduction  in  exit  stagnation  pressure 
non-dimensionalized  by  the  inlet  dynamic  pressure. 

.  _  “  PO, 

U )  —  - 

p<h  -  Pi 

where  pn7  is  the  mass-averaged  exit  stagnation  pressure  and  p»iaen  is  the  isentropic  (i.e. 
lossless)  stagnation  pressure  at  the  exit.  Two  approaches  are  used  for  calculating  loss 
—  a  mass-averaged  stagnation  pressure  defect  at  the  exit  plane,  and  an  extrapolation 
to  a  fully  mixed  condition.  Further  details  on  the  loss  calculation  are  presented  in 
Appendix  C. 


(2.14) 


2.7  Non-Dimensionalization  and  Reference  Quantities 


Although  the  basic  equations  can  be  used  with  quantities  having  any  consistent  set  of 
units,  it  is  advantageous  to  impose  a  systematic  non-dimensionalization  of  the  quantities 
defining  the  geometry  and  flow.  The  physical  quantities,  such  as  i,  etc.,  are  non- 
dimensionalized  using  reference  quantities  c,  p,,  and  a,\. 

x  y 

x  =  T—  y  =  - - 

Caxiat  Canal 


(2.15) 
ao 

This  means  that,  in  the  non-dimensional  variables,  inlet  reference  values  for  flow  vari¬ 
ables  are  defined  at  r  =  0  as 


POr«/  -  1  P<  'rrl 
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“Or «/  -  1 
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(2.16) 
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Chapter  3 

Numerical  Formulation 


This  chapter  describes  the  numerical  implementation  of  the  Newton  system  and  bound¬ 
ary  conditions  for  the  streamtube-based  Euler  solution  with  integral  boundary  layers. 
A  design  capability  for  viscous  cascade  redesign  using  a  set  of  geometric  perturbation 
modes  is  introduced.  The  basic  scheme  is  generalized  to  treat  multiple  blade  cascades 
and  several  new  modes  for  the  design  of  splitter  vanes  are  introduced.  Problems  with 
regular  topology  grids  for  supersonic  cascades  are  discussed  along  with  the  formula¬ 
tion  for  the  new  offset-periodic  grid  topology.  Additional  details  on  the  derivation  and 
implementation  are  presented  in  the  Appendices. 


3.1  Overview  of  Numerical  Scheme 

The  viscous  cascade  flow  is  solved  using  a  zonal  approach,  dividing  the  flow  into  an 
inviscid  outer  flow  modeled  using  the  steady-state  Euler  equations,  and  viscous  layers 
over  blades  and  in  wakes  modeled  with  an  integral  boundary  layer  method.  The  two 
zones  are  interacted  by  driving  the  boundary  layer  with  the  edge  flow  quantities  from 
the  inviscid  flow,  whose  edge  is  defined  by  the  boundary  layer  displacement  thickness. 
The  non-linear  equations  for  the  two  zones  are  solved  simultaneously  using  a  Newton- 
Raphson  solver  making  it  possible  to  treat  fully  separated  flows.  The  Newton  solver 
uses  a  direct  solution  of  the  linearized  system,  making  the  method  robust,  and  allowing 
additional  well-conditioned  degrees-of-freedom  (DOF’s)  and  constraint  equations  to  be 
freely  added  to  apply  useful  boundary  conditions  or  design  functions. 
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3.1.1  Inviscid  Flow 


The  inviscid  flow  is  solved  with  a  conservative  formulation  of  the  steady-state  Euler 
equations  in  the  blade-relative  coordinate  system.  The  equations  for  mass,  momentum 
and  energy  are  solved  on  an  I  x  J  streamline  grid  of  nodes  that  define  the  problem 
domain,  see  Figure  3.1.  The  grid  is  defined  by  families  of  streamlines  and  quasi-normal 
lines  indexed  by  j  by  i,  respectively. 


Figure  3.1:  Streamline  grid  system  for  discrete  Euler  equations. 


The  discrete  form  of  the  equations  is  derived  by  applying  the  integral  equations  2.2 
through  2.4  to  the  basic  conservation  cell  at  i,  j,  defined  on  a  streamtube  between  three 
quasi-normal  lines,  as  shown  in  Figure  3.2.  The  use  of  the  streamline  grid  considerably 


Figure  3.2:  Conservation  cell  and  variable  locations. 
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simplifies  the  discrete  form  of  the  equations  because  it  eliminates  convection  terms 
across  the  streamline  grid  faces.  The  velocity  and  density  need  only  be  defined  on  the 
two  quasi-normal  faces  of  the  cell.  Pressures  are  defined  on  quasi-normal  faces  as  p, 
and  on  streamwise  faces  as  11*.  The  mass  equation  reduces  to  an  implicit  statement  of 
constant  mass  flow  in  each  streamtube. 

m.\  —  m2  =  •  •  •  =  rh, 

Similarly  the  energy  equation  becomes  an  statement  of  constant  rothalpy 

Ii  =  I2  =  •  •  •  =  I, 

The  two  equations  remaining  explicitly  are  the  components  of  momentum.  The  momen¬ 
tum  forces  are  resolved  into  a  streamwise,  or  S-momentum,  equation  and  a  quasi-normal, 
or  N-momentum,  equation. 

S  ‘  ^2  fi.j  =  FSi.j{p,P,q,x,y,  r,b)  =  0  (3.1) 

N  •  F,,j  =  Fxt  j  (p,  p,  q,  n~,  n+,  x,  y,  r,  b)  =0  (3.2) 

where  the  forces  consist  of  pressure,  mass  flux  and  rotational  forces.  Rather  than  use 
the  N-momentum  equation  directly,  the  II*  streamline  pressure  variables  are  eliminated 
by  differencing  the  N-momentum  equations  on  adjacent  streamtubes.  The  resulting 
equation,  the  reduced  N-momentum  equation,  gives  the  net  force  across  the  streamline 
directly. 

?xu{piP,<l,z,y,r,b)f'  =0  (3.3) 

where  the  dependence  is  now  on  the  variables  in  two  streamtubes  j  and  j  +  1.  Note 
that  the  solution  of  the  Euler  equations  must  result  in  the  balancing  of  the  streamwise 
forces  in  each  streamtube  and  the  normal  forces  across  each  streamline. 

One  of  the  advantages  of  this  formulation  is  that  the  streamline  equations  require 
no  added  dissipation  for  subsonic  flow.  Dissipation  is  added  in  supersonic  regions  in 
the  form  of  an  upwinded  velocity  term,  similar  to  a  bulk  viscosity.  This  term  becomes 
significant  in  strong  velocity  gradients  and  stabilizes  the  scheme  for  capturing  shocks. 
The  first-order  dissipation  scheme  from  reference  (9)  is  normally  used,  although  a  second- 
order  dissipation  scheme  has  been  added  that  reduces  spurious  dissipation  loss.  Stability 
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problems  at  supersonic  Mach  numbers  have  limited  the  use  of  the  second-order  scheme 
to  subsonic  and  transonic  Hows. 

Corresponding  to  the  two  momentum  equations  for  each  conservation  cell  there  must 
be  two  unknowns.  These  have  been  chosen  as  the  changes  in  density  bp  on  each  quasi¬ 
normal  cell  face,  and  the  normal  displacement  6n  of  each  streamline  node,  as  shown  in 
Figure  3.2.  The  grid  node  positions  are  therefore  determined  as  part  of  the  solution, 
making  it  simple  to  extend  the  basic  method  to  design  problems  where  the  streamline 
shape  corresponding  to  a  specified  pressure  distribution  is  determined.  The  discrete 
form  for  the  interior  equations  is  fully  derived  in  Appendix  A. 

3.1.2  Boundary  Conditions 

The  system  of  interior  equations  must  be  closed  by  applying  boundary  conditions  at 
the  edges  of  the  domain,  see  Figure  3.3  These  take  the  form  of  periodic  and  solid  wall 


Figure  3.3:  Boundary  conditions  on  domain  boundaries. 

boundary  conditions  on  the  j  =  1  and  j  =  J  fares  of  the  domain  and  inlet  and  exit 
slope  conditions  on  the  t  =  1  and  i  =*  I  faces.  Although  they  are  not  unknowns  in  the 
equations,  quantities  for  rothalpy  and  mass  flow  are  also  implicitly  or  explicitly  set  at 
i  =  1.  Each  streamtube  also  has  a  specified  stagnation  density  condition  set  at  the  inlet. 
For  choked  flow  the  static  pressure  is  imposed  at  the  exit. 
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The  discrete  form  for  the  edge  boundary  conditions  is  summarized  in  Table  3.1.  Pe¬ 
riodic  conditions  are  applied  across  stagnation  streamlines  upstream  and  downstream  of 
the  blade,  specifying  An,,  the  gap  between  streamlines  j  =  1  and  j  =  J,  accounting  for 
the  pitch  offset.  Solid  wall  conditions  are  applied  on  each  of  the  stagnation  streamlines 
between  the  blade  leading  and  trailing  edge  by  imposing  An,  the  gap  between  stream¬ 
line  and  the  wall.  These  definitions  for  An  are  modified  for  viscous  flows  to  account  for 
blade  and  wake  boundary  layer  displacement  thickness. 


Table  3.1:  Discrete  boundary  conditions  on  domain  edges. 


Boundary  type 

Discrete  boundary  condition 

solid  surface 

o  o 

ll  II 

-i 

c  c 

<1  <1 

periodic 

An,  =  0 

nr,  -  n,t/  =  o 

subsonic  inflow 

V7.J  ~V\.j 

—  stnl 

*  l  -j 

P'h.j  ~  P«r ,/  =  0 

supersonic  inflow 

arctan  .j)  =  c»>> 

P'h.j  ~  PUrrl  =  0 

outflow 

vij  -  yi-u  . 

—  aout 

* fj  ~  */—!•> 

The  inlet  flow  direction  is  imposed  by  requiring  that  all  streamline  segments  at 
i  =  1  have  a  specified  slope  =  tan/?].  For  supersonic  inlet  conditions  the  inlet  slope 
condition  is  replaced  by  an  inlet  characteristic  condition  specifying  the  local  angle  for 
streamlines  1  <  j  <  J  -  l.  The  angle  for  the  remaining  j  =  J  streamline  is  set  by 
geometric  periodicity. 
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The  outlet  flow  direction  is  treated  analogously  to  the  inlet  direction,  and  is  imposed 
by  requiring  that  all  segments  at  i  =  I  have  slope  sout  =  tan  fa.  In  practice,  to  avoid 
a  large  pressure  jump  at  the  blade  trailing  edge,  the  exit  angle  should  be  treated  as  an 
additional  unknown  which  is  set  by  the  trailing  edge  Kutta  condition. 

The  inlet  stagnation  density  condition  is  imposed  using 

Ptj  (y)  '  '  "  P»re/  =  0  where  h,,  =  I  -  ^  (3.4) 

The  streamtube  approach  encourages  several  of  the  boundary  conditions  to  be  applied 
implicitly,  written  directly  into  the  momentum  equations.  The  rothalpy  is  constant 
in  the  domain  and  is  not  treated  as  an  unknown.  The  mass  flow  is  also  constant  for 
each  streamtube  and  is  treated  implicitly.  Mass  flow  is  set  by  inflow  Mach  number 
and  streamtube  area  at  the  inlet.  The  implicit  treatment  of  mass  flow  requires  special 
treatment  for  choked  flow,  where  the  mass  flow  is  effectively  determined  by  an  imposed 
exit  pressure.  The  details  of  the  choked  flow  solution  are  discussed  in  Appendix  D. 


3.1.3  Newton-Raphson  Solution 


The  complete  set  of  flow  equations  and  boundary  conditions  specify  a  large  non-linear 
system  which  is  solved  using  a  full  Newton-Raphson  method.  An  initial  streamline  grid, 
roughly  corresponding  to  an  incompressible  flow  solution,  is  generated  about  the  blades 
using  an  elliptic  grid  generator  and  supplies  the  initial  guess  for  the  Newton  iteration. 
The  flow  field  is  initialized  with  a  uniform  density  for  subsonic  flow,  or  a  smooth  gradient 
of  density  to  a  specified  exit  condition  at  supersonic  or  choked  conditions. 


The  basic  Newton  method  solves  the  general  non-linear  system 


R(Q)  =  0 


(3.5) 


where  Q  is  the  vector  of  variables  and  R  is  the  vector  of  equations.  At  iteration  level 
m,  the  Newton  solution  procedure  is 


0R 


0QJ 


6  Q"‘  =  -R'1 


(3.6) 
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qm-M  _  Q«  +  $Q">  (3,7) 

Care  must  be  taken  with  the  last  step  when  updating  the  variables  to  limit  changes 
to  ensure  that  the  variables  remain  within  reasonable  limits.  The  Newton  changes  are 
clamped  so  that  density  and  velocity  may  only  change  by  30%  per  iteration. 


For  the  inviscid  flow  R  consists  of  the  S- momentum  equations  on  each  streamtube 
and  the  N-momentum  equations  on  each  streamline  plus  the  boundary  conditions  at 
the  domain  edges.  The  vector  of  unknowns  Q  is  made  up  of  unknowns  for  the  interior 
equations  consisting  of  perturbations  to  the  streamtube  cell  density  8p  defined  on  quasi¬ 
normal  faces  and  perturbations  to  the  streamline  normal  node  displacement  8n.  A  few 
additional  global  unknowns  are  normally  introduced  to  control  boundary  conditions 
such  as  inlet  or  exit  angle.  One  additional  global  constraint  equation  must  be  added  to 
the  system  for  each  global  unknown. 

The  Jacobian  matrix  must  be  generated  and  inverted  for  each  Newton  step. 

This  system  represents  the  linearized  change  in  the  residual  with  respect  to  the  un¬ 
knowns  and  is  constructed  by  taking  variations  of  the  equations  with  respect  to  each 
of  flow  and  geometric  variables.  These  variations  are  resolved  by  repeated  application 
of  the  chain  rule  for  partial  derivatives  to  be  minimally  expressed  only  in  terms  of  the 
Newton  system  unknowns  8p  and  8n.  For  example,  taking  the  S-momentum  equation 
3.1  at  i,j 


d?s  c  ,  dTs  c  ,  d?s  c  ,  d?s  c  ,  d?s  c  ,  9FS  c  ,  dFs  tk 
tfsi.i  =  -jf  6p  +  -gf  sp  +  +  — Sx  +  —  6y +  — Sr  +  — 6b 


dx 


dy 


dr 


db 


The  variations  Sp,  8q,  etc.  are  resolved  in  terms  of  their  dependencies  on  6p  and  8n 
using  repeated  applications  of  the  chain  rule,  as  an  example  the  pressure  variation 
p(p,  q,  z,  y,  r ,  6)  is  resolved  using 

_  (dp  dp  dq\  (dp  dr  dx  dp  ( dq  dx  dqdy  dqdrdx\\  . 

^  \dp  dqdp)  \dr  dx  dn  dq  \5xdn  dy  dn  drdxdn)) 

where  the  intermediate  derivatives  such  as  are  generated  and  stored  at  each  iteration 
level.  Although  it  looks  awkward  on  paper,  the  chain  rule  process  is  reasonably  efficient 
at  the  Fortran  level,  at  least. 
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3.1.4  Matrix  System 


The  Jacobian  matrix  is  a  large  linear  system  of  J  x  (2 J  -  1)  equations,  made  up  of 
(I  -  2)  x  (/  -  1)  S-momentum  equations  and  (J  -  2)  x  (J  -  2)  reduced  N-momentum 
equations.  Boundary  conditions  supply  equations  on  the  j  =  1  and  j  =  J  streamlines 
and  the  i  =  1  and  i  =  I  quasi-normal  lines. 

The  domain  of  dependence  for  the  S  and  N-momentum  equations  is  relatively  com¬ 
pact,  as  shown  in  Figure  3.4.  The  basic  equations  span  three  streamwise  nodes,  corre- 


Figure  3.4:  Domain  of  dependence  for  S-momentum  and  N-momentum  equations. 

sponding  to  blocks  B,,  4;  and  C,.  The  S-momentum  equation  spans  one  streamtube  (two 
streamlines)  and  the  N-momentum  equation  spans  two  streamtubes  (three  streamlines). 
Additional  upstream  dependence,  due  to  artificial  dissipation,  increases  the  domain  of 
dependence  to  include  one  or  two  additional  streamwise  nodes.  The  small  ‘molecule’ 
of  dependence  results  in  a  sparse  Jacobian  matrix  whose  equations  are  most  efficiently 
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arranged  as  I  block  rows  of  (2J-1)  equations  in  each  block  and  have  the  form 


'  A, 

C, 

.  H,  ' 

r  i 

-Ri 

b2 

A2 

c2 

i  h2 

$2 
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'  g3. 

69  \ 

L-R-J 

The  unknown  vector  for  each  block  is  =  (Snjj—j,  corresponding  to  the 

density  changes  and  displacements  on  the  t  quasi-normal  line.  The  R,  terms  are  the 

Newton  residuals  for  the  S-mom.  and  N-mom.  equations.  Additional  global  unknowns 

6,,  and  constraint  equations,  made  up  of  G,  entries  for  the  interior  flow  unknowns, 
= 

and  Gr;  entries  for  the  global  unknowns,  do  not  fit  within  the  compact  block  structure. 
The  H,  blocks  represent  the  effect  of  the  global  unknowns  on  the  S-mom.  and  N- 
mom.  equations.  Further  details  on  the  block  structure  of  the  matrix  and  solution  are 
presented  in  Appendix  D. 


3.1.5  Global  Variables  and  DOF’s 

The  large,  sparse  matrix  in  3.8  is  solved  using  a  specialized  direct  solver  that  divides 
the  full  matrix  into  a  banded  matrix  and  a  small  matrix  for  the  global  unknowns.  The 
global  influence  columns  H  are  simply  treated  as  additional  right-hand  sides,  so  that 
the  addition  of  global  unknowns,  or  DOF’s  to  the  system,  increases  the  solution  cost 
by  a  negligible  3%  per  DOF.  This  treatment  of  global  variables  as  separate  unknowns 
with  global  constraint  equations  adds  considerable  flexibility  in  the  implementation  of 
boundary  conditions  or  design  conditions.  The  equation  added  to  constrain  a  DOF 
need  only  be  well-posed,  it  does  not  even  need  to  directly  involve  the  global  variable. 
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For  example,  the  exit  angle  sout  can  be  explicitly  set,  typically  resulting  in  a  pressure 
mismatch  at  the  blade  trailing  edge.  The  unique  (but  unknown)  exit  slope  that  satisfies 
the  Kutta  condition  can  be  found  setting  sou<  as  a  DOF  and  adding  a  pressure  equality 
constraint  equation  at  the  trailing  edge 

nrTE.i  -  Kr.j  = 0  (3-9) 

Most  of  the  quantities  that  affect  the  boundary  conditions  can  be  used  as  DOF’s  and 
have  corresponding  direct  constraint  equations,  see  Table  3.2. 

Table  3.2:  Direct  constraint  equations  for  global  DOF’s. 


DOF 

Global  Unknown 

Direct  Constraint  Equation 

S„i  1 

Ss.ni 

S.ill  ~  <S|n/.prr  =  0 

&onl 

Stout 

&ntit  3out,ptc  =:  0 

Ctnl 

6c„,i 

C,,i/  —  Cln/  lp0(.  =  0 

m 

6m 

m  -  m,p,c  =  0 

Prrtt 

6pr.nl 

Pf  .nl  ~  Pmt  =  0 

The  power  of  the  global  variable/global  constraint  arrangement  can  be  used  to  pro¬ 
vide  some  powerful  capabilities  for  cascade  work.  For  example,  the  inlet  characteristic 
global  variable  ctnj  can  be  used  to  set  the  Mach  number  into  the  infinite  cascade  when 
it  is  paired  with  a  constraint  equation  for  the  mass-averaged  Mach  number  at  the  inlet 
station  i  =  2. 

T  M,.J  -  A V,  =  0  (3.10) 

ft  m 

The  Mach  constraint  equation  can  be  generalized  and  applied  across  a  line  at  any 
specified  position  in  the  cascade.  This  is  particularly  valuable  for  determining  the  exit 
pressure  for  a  cascade  corresponding  to  a  given  flow  condition.  The  exit  pressure  prT,i 
is  set  as  a  DOF,  constrained  by  the  mass  averaged  Mach  number  constraint  applied  at 


a  specified  position  in  the  cascade.  If  the  Mach  number  is  specified  as  Mspec  =  1.0  this 
constraint  can  be  used  to  set  the  position  of  a  normal  shock,  although  this  condition 
must  be  used  with  some  care  so  that  the  desired  shock  position  and  Mach  number  is 
physically  realizable.  The  exit  pressure  to  attain  a  given  exit  Mach  number  can  also  be 
determined  using  this  technique. 

Blades  with  round  leading  edges  add  a  DOF  for  tangential  movement  of  the  leading 
edge  stagnation  point.  This  is  handled  by  sliding  the  streamline  nodes  along  the  blade 
surface  accordian-fashion  to  match  the  stagnation  node  movement.  The  leading  edge 
movement  is  implemented  as  a  a  global  grid  displacement  DOF  that  also  moves  nearby 
nodes  in  the  interior  mesh  to  minimize  local  grid  skewing.  The  displacement  of  each 
grid  node  depends  on  normal  streamline  displacements  and  on  the  global  DOF  s/,£-. 

6  {x,j  i  +  y,jj  )  =  8n,.j  hlv,  +  S3LE(xLEi.j  *  +  yLE.j  j)  (3-11) 

The  position  of  the  stagnation  point  aur  is  constrained  by  a  leading  edge  Kutta  condi¬ 
tion  similar  to  the  trailing  edge  Kutta  condition  3.9.  This  is  covered  in  more  detail  in 
reference  [9]. 

3.2  Boundary  Layer  Coupling 


Viscous  effects  are  included  by  modifying  the  wall  boundary  condition  on  blade  and 
wake  streamlines  using  the  displacement  surface  concept.  The  condition  that  fixes  the 
wall  streamline  to  the  blade  surface  from  Table  3.1  is  replaced  by 

Srii  -  S(6 ’)  =  S’  -  An,  (3.12) 

where  S’  is  the  BL  displacement  thickness,  and  An  is  the  streamline  offset  distance, 
as  shown  in  Figure  3.5.  The  same  condition  is  used  on  the  wake,  but  the  thickness  is 
applied  between  the  two  streamlines  bounding  the  wake. 

The  theory  and  numerical  implementation  of  the  boundary  layer  solution  is  given  in 
references  [9,  12,  13].  Details  of  the  quasi-3D  modifications  and  matrix  implementation 
are  contained  in  the  Appendices. 
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Figure  3.5:  Boundary  conditions  for  viscous  coupling. 


3.3  Design  Capabilities 


The  streamline- based  method  makes  possible  a  variety  of  useful  design  capabilities.  The 
classical  inverse  problem,  where  the  geometry  corresponding  to  a  prescribed  surface 
pressure  is  determined,  is  easily  implemented  with  the  current  approach  by  replacing 
the  boundary  condition  imposing  the  displacement  of  surface  nodes  with  an  imposed 
surface  pressure  condition.  Full  inverse  design  is  rarely  used  in  practice,  however,  and 
has  been  removed  from  the  current  method,  in  favor  of  mixed-inverse  design.  The 
mixed-inverse  problem,  where  one  portion  of  the  geometry  is  prescribed  and  the  other 
is  to  be  designed,  is  better  suited  to  the  iterative  nature  of  design  where  problem  areas 
are  identified  and  fixed  on  a  local  basis.  For  cases  with  areas  of  strongly  constrained 
geometry,  such  as  the  supersonic  zones  on  a  transonic  airfoil  or  the  supersonic  surface 
of  a  supersonic  cascade,  it  is  essential  that  only  localized  geometric  changes  be  made. 
Both  of  these  methods  are  fully  described  in  reference  [9].  Note  that  these  techniques 
are  only  marginally  useful  for  separated  flow  cases,  where  the  local  geometry  is  poorly 
constrained  by  the  pressure. 

For  external  aerodynamics,  where  viscous  effects  are  usually  Second-order,  design 
changes  are  made  inviscidly  to  maintain  profile  smoothness,  and  a  s,eries  of  viscous 
check  cases  run  to  assess  performance.  This  inviscid  redesign  strategy  is  likely  to  fail 
for  cascades  due  to  the  strong  viscous  coupling,  involving  some  degree  of  separated  flow 


at  the  operating  point1 .  An  alternative  strategy  is  to  deal  with  a  fixed  set  of  geometric 
modes  to  perturb  the  airfoil  contour,  assuring  that  any  design  changes  remain  smooth 
without  locally  wild  behavior.  For  a  modal  scheme  the  surface  node  displacement  is 

given  by  a  sum  of  geometric  perturbations 

L 

Sn i  =  6fi{  gt(8t)  (3.13) 

r=\ 

where  the  L  geometric  modes  gi(s)  are  defined  in  terms  of  surface  arc  length  and  have 
mode  amplitudes  hi.  Seven  geometric  modes  are  implemented  in  the  current  method  and 
are  used  in  for  a  redesign  case  in  Chapter  5.  The  mode  shapes  are  shown  in  Figure  3.6, 
although  any  reasonably  orthogonal  set  of  perturbations  could  be  used.  The  modes  can 
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* 


Camber-line 


perturbation 


< 


rSn  --  £8n  gls) 


5 ingle -surface  perturbation 


Figure  3.6:  Mode  shapes  and  design  modes  for  geometric  perturbations. 

be  used  to  perturb  all  or  any  specified  region  of  a  single  blade  surface,  or  the  same 
perturbation  could  be  applied  to  opposing  sides  on  both  blade  surfaces.  The  latter 
possibility  is  equivalent  to  a  camber  redesign  mode,  preserving  blade  thickness.  The 
mode  amplitudes  hi  can  be  set  as  global  DOF’s  with  global  unknowns  Shf.  The  modes 
can  be  specified  directly  using  the  corresponding  direct  modal  constraint  equations 

=  0  (3.14) 

'The  possible  exception  being  for  blading  designed  to  avoid  separation,  e.g.  DCA  sections. 
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For  modal-design  cases,  the  surface  displacement  boundary  condition,  3.12,  is  mod¬ 
ified  with  an  extra  term  corresponding  to  the  geometric  modes. 

/, 

Snt  -  S(S-)  =  6‘  -  An,  +  ^  Snf  g,(s,)  (3.15) 

l=  i 

Note  that  the  streamline  is  still  offset  from  the  airfoil  by  the  boundary  layer  thickness, 
but  the  blade  surface  can  be  displaced  by  the  geometry  modes. 

The  inverse  modal  formulation  minimizes  the  integrated  mismatch  between  a  spec¬ 
ified  and  resulting  surface  pressure  distribution.  The  discrete  form  for  this  is 

1  =  ~  Ps, .«■«■( a,')‘  (3-16) 

I 

summed  only  over  the  inverse  segment.  Minimization  constraints  are  obtained  in  the 
standard  way  by  setting  the  variations  with  respect  to  the  mode  amplitudes  to  zero. 


3.4  Multiple  Blades 


One  of  the  goals  for  the  current  work  is  the  development  of  methodology  for  analysis  and 
design  of  cascades  with  splitter  vanes.  The  problem  was  approached  from  a  general  point 
of  view  so  that  anv  number  of  intermediate  blades  rnav  be  used,  although  the  discussion 
and  examples  focus  on  one  splitter  vane.  The  basic  configuration  for  a  cascade  with 
splitter  vane  is  shown  in  Figure  3.7.  The  intermediate  blade  introduces  an  additional 
stagnation  line  that  splits  at  the  leading  edge  to  define  the  upper  and  lower  surface 
streamlines  with  associated  boundary  layers.  This  stagnation  line  also  acts  as  a  dividing 
streamline  between  the  upper  and  lower  passages,  each  with  its  own  mass  flow.  The 
relative  passage  mass  flows  are  set  by  the  position  of  the  stagnation  line  at  the  inflow 
boundary,  which  is  in  turn  set  by  the  trailing  edge  Kutta  condition,  All  =  0  on  the 
intermediate  blade. 

The  streamline  grid  structure  is  modified  for  a  splitter  vane  as  shown  in  Figure  3.8, 
where  the  stagnation  streamline  for  the  intermediate  blade  defines  two  grid  stream¬ 
lines  and  a  null  streamtube  through  which  there  is  no  flow.  The  conditions  on  the  new 
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Figure  3.7:  Cascade  with  splitter  blade,  showing  stagnation  lines. 


Str#amtub*s 


Figure  3.8:  Grid  arrangement  for  multiple  blades. 
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stagnation  streamlines  are  handled  by  generalizing  the  periodic  boundary  conditions  in 
Table  3.1  for  the  j  =  1  and  j  —  J  domain  boundaries  to  apply  across  adjacent  stream¬ 
lines.  Further  details  on  the  treatment  of  multiple  blades  are  covered  in  Appendix  D. 

3.4.1  Global  Variables  and  DOF’s  for  Multiple  Blades 

The  stagnation  line  for  the  intermediate  blade  divides  the  cascade  into  upper  and  lower 
passages  with  mass  flows,  m/„.r  and  m„pr.  This  effectively  introduces  a  new  global 
unknown  for  mass  flow  into  the  system.  It  is  convenient  to  treat  the  two  unknown 
mass  flows  as  a  total  mass  flow  m  =  mi,,.,  ~  m,()„  and  a  differential  mass  flow  Am  - 
m„,„  -  rh/„.r.  This  allows  the  total  mass  flow  to  be  imposed  from  upstream  conditions 
and  the  differential  mass  flow  determined  by  a  single  added  DOF.  The  global  unknown 
for  differential  mass  flow  is  normally  constrained  by  a  global  equation  for  the  trailing 
edge  Kutta  condition  on  the  intermediate  blade.  It  may  also  be  directlv  specified  by  a 
constraint  equation  for  the  mass-fraction-ratio  nit,  —  rh/,,,.. 

^  rn  (, ..  .  n  (3.17, 

Till.,  I 


Each  additional  blade  with  a  blunt  leading  edge  also  requires  a  DOF  and  constraint 
equation  for  LE  stagnation  point  movement,  similar  to  the  main  blade.  The  global 
grid  movement  mode  that  accompanies  leading  edge  movement  extends  only  over  the 
passages  immediately  above  and  below  the  blade. 

Two  new  design  modes  have  been  added  for  intermediate  blade  redesign  -  a  blade 
displacement  mode  Dmm.  and  a  blade  shearing  mode  Drnl.  The  movement  mode  moves 
the  blade  in  the  circumferential  direction  and  is  used  to  change  the  mass  flow  fraction 
between  the  flow  passages.  The  shearing  mode  maintains  axial  chord  and  trailing  edge 
position  and  is  used  primarily  for  modal  design  and  modal  inverse  redesign.  Figure  3.9 
shows  the  associated  mode  shapes  at  the  splitter  blade.  These  modes  differ  from  the 
modal  design  modes  in  Figure  3.6  because  they  are  not  limited  to  modal  displacements 
along  the  streamline  normal  direction  and  they  are  global  in  the  sense  that  they  specify 
displacements  for  every  grid  node.  The  global  modes  are  treated  in  a  similar  way  as 
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the  leading  edge  movement  mode.  In  practice  any  reasonably  well-defined  mode  shape 
could  be  used  in  place  of  the  movement  or  shearing  modes.  The  form  for  the  node 
displacement  equation  includes  both  normal  displacements  Sn  and  global  displacements 
from  movement  and  rotation  modes  which  specify  x  and  y  nodal  displacements  for  the 
entire  grid,  for  example 

4"  y,.j  ;  )  =  Sn, ,j  fti.j  4"  $ Dmov{%movjj  *  4”  J/mor;  j  j)  (3.18) 

The  new  modes  can  be  set  as  DOF’s  and  have  direct  constraint  equations  similar  to  the 
direct  modal  constraints  3.14.  The  intermediate  blade  global  modes  may  be  mixed  with 
the  geometric  modes  in  inverse  design  calculations  where  an  inverse  pressure  constraint 
similar  to  3-16  is  used.  Alternatively  the  mass  fraction  constraint  in  3.17  can  be  used 
to  drive  the  blade  position  to  give  a  specified  mass  flow  ratio  across  the  blade. 


Figure  3.9:  Design  modes  for  splitter  movement  and  rotation. 

3.4.2  Grid  Generation 

Since  the  velocity  and  pressure  are  determined  by  streamtube  geometry,  a  streamline 
grid  generator  is  used  to  supply  the  initial  guess  for  the  flow  field.  Grid  generation 
for  a  single  blade  cascade  is  relatively  straightforward  —  once  a  distribution  of  nodes 
is  established  on  the  domain  boundaries  an  elliptic  grid  generator  is  used  to  establish 
the  position  of  interior  nodes,  determining  the  streamlines  and  quasi-normal  lines.  The 
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normal  spacing  of  streamlines  is  set  by  the  fractional  mass  flow  through  each  streamtube. 
The  resulting  grid  essentially  corresponds  to  an  incompressible  flow  solution  for  the 
assumed  positions  of  the  inflow  and  exit  stagnation  lines. 

The  multiple  blade  grid  is  generated  in  similar  fashion,  treating  each  of  the  flow 
passages  as  an  independent  domain,  bordered  by  stagnation  streamlines.  Two  difficulties 
are  introduced  for  multiple  blades. 

1.  The  position  of  the  intermediate  stagnation  lines  is  not  known  a  priori. 

2.  The  initial  node  spacing  on  stagnation  lines  strongly  affects  grid  skew  in  the 
interior  mesh.  Grid  skew  has  little  effect  in  subsonic  flow  but  can  lead  to  serious 
problems  for  supersonic  cascades  (this  will  be  discussed  in  the  following  section). 

The  approach  used  for  cascades  assumes  the  main  and  intermediate  blades  are  not  closelv 
spaced  (as  they  would  be  in  a  multi-element  airfoil)  and  that  a  cascade  ‘meanline'  based 
on  axial  position  can  be  defined  using  the  stagnation  line  and  camberline  of  the  main 
blade.  The  position  of  an  intermediate  stagnation  line  ran  be  determined  from  the 
‘mealtime'  using  the  relative  circumferential  position  of  the  intermediate  blade  at  the 
leading  or  trailing  edge.  i.e. 

,  ,  ,  ,  .  yi  i.  -  y<„.i .(*/  / ) 

y.t.n(r)  =  y„ ,./.(*)  -  ny,:hi,  with  rj  - - 

I, 

The  resulting  stagnation  line  is  smoothed  to  even  out  any  remaining  discontinuities. 

Spacing  on  stagnation  lines  is  handled  using  a  rurvature-based  node  distribution 
on  the  blades  and  a  primitive  blocking  scheme  to  roughly  match  spacing  on  inlet  and 
exit  streamlines.  This  sometimes  requires  user  interaction  to  achieve  a  satisfactory  grid 
but  is  largely  automatic.  As  the  blades  become  more  closely  spaced,  e.g.  the  ARL 
tandem  splitter  in  Chapter  5,  the  spacing  requirements  become  more  draconian  and 
require  greater  user  intervention.  For  such  cases  a  different  approach  may  be  needed, 
subdividing  the  domain  into  grid  blocks  with  an  arbitration'  scheme  for  resolving  spac¬ 
ing  requirements  between  blocks  to  obtain  a  compatible  grid  layout.  A  variant  of  the 
multiple  blade  cascade  code  has  been  applied  to  the  multi-element  airfoil  problem  by 
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Drela  and  uses  a  panel  method  solution  to  determine  the  initial  stagnation  line  positions 
and  node  distribution. 


3.5  Leading  Edge  Problems 


The  leading  edge  radii  for  cascade  blades,  especially  supersonic  blading,  is  typically  much 
smaller  than  for  airfoils  —  0.5%  for  the  DFVLR  blade  and  0.25%  for  the  ARL  supersonic 
cascade  discussed  in  Chapter  4.  Problems  are  sometimes  encountered  in  resolving  the 
flow  near  the  leading  edge  stagnation  point,  particularly  due  to  spurious  velocity  or 
total  pressure  loss  from  dissipation  terms  that  become  significant  in  the  strong  velocity 
gradients  near  the  leading  edge.  The  small  radii  can  produce  gradients  that  are  difficult 
to  distinguish  from  shocks,  although  the  actual  flow  is  nearly  isentropic.  Decelerating 
flow  produces  spurious  loss,  while  accelerating  flow  can  produce  a  total  pressure  increase. 
Typically  only  the  first  few  streamtubes  immediately  adjacent  to  the  blade  are  strongly 
affected.  Unfortunately  loss  or  gain  in  the  inner  stream  tube  has  a  strong  effect  on  the 
boundary  layer  through  changes  in  the  edge  velocity.  Overall  spurious  loss  effects  thus 
depend  on  details  of  the  alteration  in  edge  velocity  and  on  the  inviscid  loss. 

One  solution  to  this  problem  is  to  modify  the  leading  edge  geometry  to  alleviate 
the  strong  gradients  by  ellipticizing  or  sharpening  it,  although  some  of  the  genuine 
physics  may  be  lost,  particularly  for  supersonic  blades  with  bow  shocks.  The  ARL 
case,  for  example,  shows  a  6%  loss  in  the  innermost  streamtubes  due  to  the  bow  shock, 
although  this  represents  a  much  smaller  portion  of  the  mass-averaged  total  loss.  An 
alternative  solution,  implemented  as  an  option,  is  to  treat  a  region  near  the  leading  edge 
isentropically,  completely  eliminating  dissipation  loss  or  gain.  The  special  treatment  is 
limited  to  the  two  streamtubes  immediately  above  and  below  the  blade  and  extends 
from  just  downstream  of  a  bow  shock  to  a  point  behind  the  leading  edge  expansion. 
This  allows  the  loss  in  the  bow  shock  to  be  correctly  modeled  and  minimizes  spurious 
dissipation  effects. 
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3.6  Grid  Skew  Problems  for  Supersonic  Cascades 


While  the  current  approach  for  solving  the  inviscid  flow  works  well  for  subsonic  and 
transonic  flows,  and  works  well  for  supersonic  airfoils,  problems  were  encountered  with 
supersonic  cascades  at  high  stagger  angles.  Grids  with  highly  skewed  cells  were  found  to 
exhibit  spurious  odd-even  instabilities  in  the  flow  field  at  supersonic  conditions.  These 
problems  were  associated  with  strong  passage  shocks  or  the  leading  edge  wave  system. 
Figure  3.10  shows  the  converged  inviscid  result  for  a  sharp  leading  edge  blade  with 
50°  stagger.  The  waves  upstream  of  the  shock  result  from  dispersion  problems  with 
the  numerical  scheme,  similar  to  those  found  by  Shapiro  [14]  for  the  linearized  steady- 
state  Euler  equations.  Shapiro  showed  that  the  dispersive  behavior  of  several  discrete 
operators  were  well  predicted  by  his  Fourier  analysis  but  offered  no  wisdom  on  reducing 
the  effect.  Higher  dissipation  reduces  the  dispersive  behavior  only  for  mild  cases,  but 
results  in  increased  shock  smearing.  The  geometry  of  the  streamlines  and  quasi-normals 


Figure  3.10:  Dispersion  effects  at  shock  for  skewed  grid. 


at  the  instability  in  a  skewed  grid  is  shown  in  Figure  3.11.  The  strongly  oblique  initial 
grid  is  further  distorted  at  alternate  streamline  nodes,  where  the  local  Mach  angle 
roughly  matches  the  mean  cell  skew. 
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Figure  3.11:  Supersonic  odd-even  instability  for  skewed  grids. 

The  problem  becomes  more  serious  as  the  skew  increases,  leading  to  non-convergence 
or  totally  inconsistent  solutions.  The  grid  and  Mach  contours  for  a  round  leading  edge 
blade,  similar  to  the  sharp  edged  blade  above,  is  shown  in  Figure  3.12.  The  Mach 
contours  show  a  ‘lozenge’  of  accelerated  flow  (M  =  1.T)  located  below  and  upstream  of 
the  blade  leading  edge  where  the  skew  is  highest.  Note  that  the  passage  shock  is  also 
badly  smeared  by  the  axial  elongation  of  cells  in  the  skewed  mesh. 
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3.7  Offset-Periodic  Grid  and  Solution 


Using  the  direct  approach  to  problem  solving,  the  skewed  grid  problems  were  solved  by 
un-skewing  the  mesh  in  the  inlet  area.  The  topology  of  the  grid  was  changed  to  make  the 
quasi-normals  roughly  orthogonal  to  the  streamlines  in  the  inflow  region.  This  results 
in  a  ‘stairstep’  grid,  as  shown  in  Figure  3.13,  where  the  j  =  1  and  j  -  J  domain  edges 
overlap  but  are  shifted  by  a  streamwise  offset  i„jj. 


Figure  3.13:  Layout  of  offset-periodic  grid. 

Skipping  the  details  of  the  formulation  for  the  moment,  Figures  3.14  and  3.15  show 
results  for  the  cascade  from  Figure  3.12  modeled  with  an  offset-periodic  grid.  The  dis¬ 
persion  problems  are  absent  and  the  shock  is  well  resolved  as  a  result  of  the  approximate 
alignment  of  the  quasi-normals  with  the  shock. 
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distribution  for  super 


Figure  3.14:  Pressure 


3.7.1  Offset-Periodic  Matrices  and  Solution 


The  offset  grid  is  configured  with  the  offset  portion  of  the  grid  extending  from  the 
inflow  boundary  to  the  trailing  edge,  1  <  i  <  ijE •  The  outflow  conditions  are  thus  not 
affected.  Periodicity  for  the  stairstep  inlet  is  achieved  with  a  modification  of  the  periodic 
conditions  from  Table  3.1.  The  An,  =  0  geometric  matching  condition  is  applied  across 
nodes  i  on  streamline  J  and  i  -  i0jj  on  streamline  1.  Similarly  for  the  pressure  matching 
condition,  where  n,_,o//.i  -  II,../  =  0.  The  offset  inlet  conditions  alter  the  arrangement 
of  terms  in  the  block  matrix  3.8  by  moving  these  periodic  terms  to  off-diagonal  blocks, 
as  shown  in  Figure  3.16. 


Figure  3.16:  Arrangement  of  off-diagonal  periodicity  terms  for  offset  grid. 

The  off-diagonal  terms  complicate  the  matrix  solution  somewhat  and  require  mod¬ 
ifications  to  the  specialized  direct  solver.  Note  that  to  minimize  fill-in  of  the  upper 
triangle  during  the  matrix  solution  the  indicial  order  of  streamlines  in  Figure  3.13  must 
be  used.  Further  details  on  the  offset  matrix  and  solver  are  included  in  Appendix  D. 
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3.7.2  Offset-Periodic  BC 


The  new  grid  topology  requires  that  boundary  conditions  be  imposed  on  the  free  bound¬ 
ary  streamline  at  j  =  1  in  the  first  offset  block.  The  conditions  specified  on  this  stream¬ 
line  segment  depend  on  the  inlet  Mach  number. 

Subsonic  The  local  angle  of  j  =  1  streamline  segments  is  prescribed  to  be  s,„i,  the 
same  condition  used  at  the  i  =  1  face.  The  offset  segments  should  be  located  far 
enough  upstream  that  they  are  relatively  unaffected  by  the  disturbance  from  the 
blade,  implying  that  subsonic  cases  should  use  at  least  three  offset  inlet  blocks. 

Supersonic/Axially  Subsonic  The  local  angles  on  the  j  —  l  streamline  segments 
are  set  by  the  inlet  characteristic  c,„/.  An  area  behind  the  i  =  1  inlet  face  is  in 
a  forbidden-signal  zone  with  respect  to  the  wave  system  from  the  blades,  so  the 
local  streamline  angles  on  the  i  =  1  face  are  set  to  match  the  angle  on  the  i  -  i,,ff 
streamline  segment  on  ;  =  1. 

3.7.3  Offset-Periodic  Grid  Generation 

Initial  grid  generation  for  the  offset  grids  uses  the  concept  of  an  inlet  offset  block,  which 
becomes  the  basic  unit  determining  the  length  of  the  inlet  and  the  number  of  quasi- 
normal  lines.  The  leading  edge  index  on  the  upper  and  lower  surface  streamlines  of  the 
main  blade  are  set  by 

tLEj=  t  =  'LEj.j  +  *o//  with  iLEj=x  ~  ^blocks  i0Jf 

The  offset  index  t„//  is  set  by  the  position  on  the  j  -  J  upper  surface  of  the  quasi¬ 
normal  line  from  the  lower  surface  leading  edge.  The  spacing  of  nodes  on  the  streamline 
edges  of  each  offset  block  is  prescribed  by  its  opposite  upstream  neighbor  to  maintain 
the  shifted  periodicity.  One  of  the  disadvantages  of  the  offset  grid  is  that  it  is  generally 
wasteful  of  streamwise  resolution  over  the  blades.  This  is  chiefly  a  problem  at  inlet 
Mach  numbers  near  unity,  where  inlet  waves  are  nearly  normal  to  the  inflow  and  longer 
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inlets  are  required.  Fortunately,  resolution  in  the  streamwise  direction  is  cheap,  0(1), 
due  to  the  banded  matrix  structure. 
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Chapter  4 

Results  -  Analysis  Cases 


Previous  chapters  have  described  the  development  of  a  method  for  quasi-3D  blade- 
to-blade  analysis  and  design.  Since  this  effort  is  an  extension  of  the  ISES  algorithm 
developed  by  Giles  [15]  and  Drela  [9]  and  duplicates  the  characteristics  of  the  original 
method  for  their  subsonic  and  transonic  test  cases,  only  results  pertaining  to  the  new 
capabilities  will  be  presented.  Effects  of  splitter  position  and  redesign  cases  will  be 
covered  in  the  following  chapter.  The  principal  objectives  for  the  test  cases  presented 
here  are 

•  Investigate  accuracy  of  loss  prediction  for  a  range  of  inlet  Mach  numbers. 

•  Validate  the  quasi-3D  capabilities  for  AVDR  (streamtube  width)  effects. 

•  Demonstrate  the  multiple  blade  capability  on  a  highly  loaded  supersonic  cascade 
with  splitter  blade. 

In  this  chapter  the  method  is  exercised  on  several  test  cases  for  purposes  of  comparison 
and  validation  of  the  basic  analysis  capability.  These  cases  do  not  include  rotational 
effects. 

•  UTRC  Build  I  -  A  subsonic,  attached  flow  cascade  at  a  Mach  number  of  0.113 
for  low  speed  loss  comparison. 

•  DFVLR  L030-4  -  A  transonic  and  low  supersonic  compressor  cascade  at  Mach 
numbers  from  0.82  to  1.10  and  a  range  of  AVDR. 

•  ARL  Splitter  -  A  supersonic,  highly  loaded  compressor  cascade  with  splitter 
vane  at  an  inlet  Mach  number  of  1.46  with  strong  streamtube  contraction. 
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4.1  General  Effect  of  AVDR 


The  axial  velocity- density  ratio,  or  AVDR,  is  related  to  the  streamtube  width  by 

AVDR  =  (4.1) 

p  |?.ri  b2 

The  conservation  of  y-momentum  dictates  that  qy ,  =  is  constant,  defining  a 

relation  between  AVDR  and  exit  angle. 

p-iK 

tan/?2  =  where  K  is  a  constant  set  by  inlet  conditions  (4-2) 

The  effect  of  AVDR  is  to  change  axial  velocity  at  the  cascade  exit,  affecting  both  the 
inviscid  flow  (exit  angle  and  pressure)  and  the  boundary  layer  development.  The  effect 
on  the  viscous  flow  is  especially  pronounced  for  highly-loaded  blades,  where  boundary 
layers  at  or  beyond  the  point  of  separation  show  strong  effects  on  blockage  and/or  loss 
from  a  small  relief  in  adverse  velocity  gradient.  In  a  subsonic,  unchoked  cascade  AVDR 
uniquely  determines  both  exit  pressure  and  exit  angle.  In  a  choked  or  supersonic  cascade 
AVDR  determines  only  the  exit  angle,  since  the  exit  pressure  is  not  uniquely  set  by  inlet 
conditions  and  geometry. 


4.2  UTRC  Low  Speed  Cascade 

The  intent  of  this  comparison  is  to  validate  the  cascade  loss  model  for  a  low  speed  case 
without  separation  or  strong  compressible  effects.  Data  was  obtained  from  a  M  =  0. 1 13 
test  of  a  DCA  blade  (Build-I)  by  UTRC  in  Reference  [16].  A  conventional  132x20  grid 
was  used  for  this  case,  as  shown  in  Figure  4.1.  Linear  streamtube  width  variation  was 
assumed,  with  an  AVDR  =  1.023.  Free  transition  was  used  for  the  MISES  comparison, 
as  experimental  pressure  data  indicated  that  separation  bubbles  were  present  on  both 
upper  and  lower  surfaces,  not  surprising  in  light  of  the  low  Reynolds  number  of  4.78  x  10'. 
The  surface  pressures  compare  well  to  test  data,  with  good  definition  of  both  upper  and 
lower  separation  bubbles,  see  Figure  4.2.  The  loss  comparison  is  excellent,  although  the 
predicted  turning  is  slightly  greater  than  the  test  ( ~2.2C). 
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4.3  DFVLR  Cascade 


This  case  analyzes  the  flow  in  a  two-dimensional  low  turning  compressor  cascade  at 
transonic  and  low  supersonic  inlet  Mach  numbers.  Comparison  data  was  obtained  from 
Schreiber&Starken  [17]  with  additional  detailed  data  for  two  specific  supersonic  oper¬ 
ating  points  from  Reference  [18].  This  data  includes  overall  performance  data  for  Mach 
numbers  from  0.82  to  1.1  at  a  range  of  inlet  angles  and  AVDR.  Surface  pressures  were 
measured  at  several  AVDR  for  the  same  inflow  conditions,  making  this  a  useful  test  case 
for  quasi-3D  effects. 

4.3.1  Section  Characteristics 

The  cascade  corresponds  to  the  45%  height  section  from  a  rotor  designed  with  MCA 
profiles  for  a  1.51  pressure  ratio  with  a  tip  Mach  number  of  1.38.  The  blade  section, 
shown  in  Figure  4.3,  has  a  design  Mach  number  of  1.09  and  stagger  angle  £  =  48.51c. 
The  blade  definition  from  Reference  [18  was  used,  with  a  circular  leading  edge  and 
thick  trailing  edge  (0.695%). 

4.3.2  Experimental  Flow  Conditions 

The  testing  configuration  used  for  the  DFVLR  cascade  presents  several  problems  for 
comparison  of  experimental  and  computational  results. 

•  The  exact  inlet  conditions  at  the  cascade  were  not  measured  but  were  assumed  to 
correspond  to  far-upstream  conditions. 

•  For  supersonic  flow,  no  upstream  wedge  or  similar  device  was  used  to  establish 
periodic  flow.  The  first  blade  set  the  inlet  flow  field  for  the  remaining  blades. 

•  Sidewall  and  endwall  suction  was  used  to  minimize  sidewall  boundary  layer  effects 
and  to  change  AVDR  through  the  cascade! 
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0.0 


Figure  4.3:  Blade  section  for  DFVLR  tests 


As  a  result,  the  inlet  flow  conditions  corresponding  to  an  “infinite”  cascade  are  not 
known.  Experimental  exit  pressures  were  not  supplied,  except  for  the  final  two  cases, 
so  supersonic  cases  were  compared  by  adjusting  exit  pressures  to  match  experimental 
pressure  distributions  at  the  blade  trailing  edge. 

The  Reynolds  number  for  these  tests  was  1.45  x  106  and  turbulent  flow  conditions 
were  assumed,  using  transition  at  2%  chord  for  upper  and  lower  surfaces. 
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4.3.3  Comparison  of  AVDR  Effects 


Since  Mach  numbers  for  this  comparison  were  close  to  M  =  1.0,  a  relatively  long  inlet 
length  was  used  to  define  the  grid  for  the  MISES  computations.  A  grid  density  of 
220x20  was  used  for  this  case  with  4  offset  blocks  in  the  inlet  region,  as  shown  in 
Figure  4.4.  The  streamtube  width  between  leading  edge  and  trailing  edge  was  assumed 
to  be  linear.  Isentropic  treatment  was  used  near  the  blade  leading  edges  to  minimize 
spurious  dissipation  losses  near  the  stagnation  point. 


Figure  4.4:  Offset-periodic  220x20  grid  for  DFVLR  cases 


Surface  pressure  data  at  inlet  Mach  numbers  of  0.82,  0.92,  1.03  and  1.10  are  com¬ 
pared  in  Figures  4.5  through  4.8  In  each  figure  experimental  and  computed  pressures  are 
compared  for  two  AVDR  values.  The  upstream  inlet  angle  was  constant  =  58.5°  for 
all  Mach  numbers.  The  pressure  coefficient  C was  used  for  this  comparison,  defined 
by 


Pot  ~  P 
Po,  -  Pi 


(4.3) 
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Figure  4.5:  AYDR  Effect  on  DFYLR  blade:  M  =  0.82,  =  58.5°. 


Transonic  Cases 


Figures  4.5  and  4.6  compare  surface  pressures  for  transonic  inlet  Mach  numbers  of 
M  —  0.82  and  M  =  0.92.  For  these  cases  the  experimental  inlet  angle  and  Mach 
number  were  used  as  boundary  conditions  for  MISES  results.  Surface  pressures  for  the 
M  =  0.82  case  compare  closely  at  the  tested  AVDR.  At  Mach  0.92  the  lower  surface 
and  trailing  edge  pressures  were  best  matched  to  test  data  by  using  a  slightly  higher 
AVDR  (1.04).  The  general  agreement  for  the  pressures  is  good,  although  the  test  data 
for  M  =  0.92,  AVDR  =  1.17  shows  a  longer  supersonic  region  on  the  blade  upper 
surface.  The  effective  inlet  Mach  number  may  be  higher  than  0.92,  possibly  resulting 
from  suction  effects  on  the  sidewall  boundary  layer  changing  the  upstream  streamtube 
width. 
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Figure  4.6:  AYDR  Effect  on  DFYFR  blade:  M  =  0.92.  3\  =  58. 5:. 
Supersonic  Cases 

Figures  4.7  and  4.8  compare  MISES  results  with  test  data  for  the  supersonic  M  -  1.03 
or  M  -  1.10  cases.  Since  exit  pressures  were  not  specified  in  Reference  17  .  they  were 
adjusted  to  match  surface  pressures  at  the  blade  trailing  edge.  A  global  DOF  using 
c,„i  to  drive  mass  flow  was  used  to  match  mass  flow  from  the  upstream  test  conditions, 
resulting  in  a  different  Mach  number  and  inlet  angle  at  the  infinite  cascade  due  to  the 
blade  wave  system.  MISES  pressures  were  re-normalized  to  the  test  condition  dynamic 
pressures  to  permit  comparison  of  Cp„ . 

The  pressures  compare  well  except  for  M  =  1.10,  AVDR  =  1.15  where  the  lower 
surface  shock  near  the  leading  edge  indicates  that  the  specified  mass  flow  may  be  ex¬ 
cessive.  Alternate  boundary  conditions,  using  a  specified  cm/  =  [cin/]t(p,(rratn  inflow- 
condition,  made  very  little  change  in  the  results,  although  the  loss  levels  were  slightly 
lower,  ilui  -  -0.002. 
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Figure  4.7:  AYDR  Effect  on  DFVLR  blade:  M  =  1.03,  0,  =  58.5° 


Figure  4.8:  AVDR  Effect  on  DFVLR  blade:  M  =  1.10,  0{  =  58.5°. 
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4.3.4  Loss  Comparison 


The  loss  data  for  the  flow  conditions  in  Figures  4.5-  4.8  are  compared  with  experimental 
loss  data  in  Figure  4.9.  The  differences  in  loss  for  low  AVDR  =  1.0  -  1.04  and  high 


1  p-'  f  •  .Itti  \  l.t*  it 

Figure  4.9:  Comparison  of  predicted  and  measured  loss  for  DFYLR  cascade 

AVDR  =  1.10  -  1.17  are  relatively  small  and  the  predicted  losses  are  in  generally  good 
agreement  with  tested  losses.  The  lower  loss  levels  predicted  for  M  -  0.92  are  consistent 
with  the  differences  in  the  upper  surface  supersonic  region  observed  in  the  pressure 
comparisons. 

4.3.5  Additional  Supersonic  Test  Cases 

Data  for  two  additional  supersonic  operating  points  from  Reference  [18]  are  compared 
to  MISES  results  in  Figure  4.10  and  4.12  for  upstream  Mach  numbers  of  1.023  and 
1.086,  respectively.  For  these  cases  the  exit  pressures  were  specified,  although  again 
only  far-upstream  inlet  conditions  were  given.  Like  the  previous  supersonic  cases,  a 
DOF  was  used  to  match  the  far-upstream  mass  flow. 

The  M  =  1.023  case  was  run  at  a  low  inlet  angle  d|  =  56.8°,  near  the  spill  point 
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for  the  cascade,  with  an  AVDR  =  1.092.  The  pressure  ratio  was  set  from  the  test 
conditions.  The  pressure  data  comparison  is  shown  in  Figure  4.10,  with  Mach  contours 
in  Figure  4.11.  Considering  the  match  in  the  pressure  data,  the  difference  between  the 
MISES  and  experimental  loss  levels  is  surprising. 

The  M  =  1.086  case  was  run  at  inlet  angle  f3\  =  58.5°,  much  closer  to  the  design 
angle  for  the  blade,  and  at  a  high  AVDR  =  1.184.  The  pressure  ratio  was  adjusted 
upward  +0.04  from  the  test  conditions  to  match  the  trailing  edge  and  lower  surface 
pressures.  The  pressures  are  compared  in  Figure  4.12  and  the  Mach  contours  are  shown 
in  Figure  4.13.  MISES  again  under-predicted  the  loss,  this  time  by  0.0252. 

Considering  the  excellent  match  in  surface  pressures,  and  the  previous  favorable 
match  in  losses  for  the  AVDR  comparison,  the  magnitude  of  the  mismatch  in  loss  levels 
is  surprising.  It  is  especially  puzzling  considering  the  typical  mismatch  between  ISES 
results  and  tested  drag  levels  is  less  than  +0.00101  for  transonic  airfoils  with  similar 
upper  surface  Mach  numbers  to  this  cascade.  The  differences  in  loss  may  result  from 
the  simplistic  modeling  of  of  the  sidewall  contraction.  The  sidewall  boundary  layers  are 
actually  subject  to  different  flow  conditions  on  each  side  of  the  blade.  Although  a  simple 
blockage  correction  using  the  total  sidewall  BL  thickness  to  set  AVDR  is  adequate  for 
subsonic  flow,  it  cannot  represent  the  physics  as  the  flow  becomes  sonic  and  the  local 
sensitivity  of  the  flow  to  blockage  effects  increases.  At  this  point  separate  blockage 
effects  are  required  for  upper  and  lower  surfaces. 

The  sidewall  boundary  layers  also  strongly  influence  the  apparent  thickness  of  the 
mid-span  boundary  layer  in  a  two-dimensional  test.  The  mechanism  for  this  is  simply 
continuity,  changing  the  boundary  layer  thickness  as  a  result  of  the  local  effective  side- 
wall  contraction  downstream  of  the  blade.  Since  the  drag  (or  loss)  is  typically  measured 
by  a  survey  of  the  wake  boundary  layer  at  the  mid-span  of  the  tested  blade,  the  lo¬ 
cally  increased  thickness  can  indicate  falsely  high  losses.  Overall,  the  best  approach  for 
transonic  or  supersonic  testing  is  to  eliminate  the  sidewall  boundary  layer  completely 
through  sidewall  suction  over  the  blade. 

1  For  carefully  conducted  tunnel  tests,  such  as  the  AGARD  RAE  2822  series. 
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Figure  4.12:  DFVLR  blade:  M  =  L.0B6,  /3,  =  58.5°,  AVDR  =  1.184 


Figure  4.13:  Mach  contours  for  DFVLR  blade,  M  =  1.086,  contour  interval  0.05 
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4.4  ARL  Supersonic  Cascade 


This  case  analyzes  the  flow  in  a  two-dimensional  high  turning  supersonic  compressor 
cascade  with  splitter  vane  at  an  inlet  Mach  number  of  1.46.  This  cascade  was  designed 
to  achieve  high  turning  by  high  streamtube  contraction  from  inlet  to  exit  (a  metal-metal 
AVDR  =  2.053!).  The  blade  and  splitter  vane  sections  are  shown  in  Figure  4.14  with 
the  splitter  at  the  nominal  50>o  circumferential  position. 


Figure  4.14:  Blade  and  splitter  vane  for  ARL  cascade 

Comparison  data  was  obtained  from  References  [2,  3,  5j.  This  data  included  turning 
and  loss  data  at  A/  =  1.46  for  a  range  of  pressure  ratios,  both  with  and  without  the 
splitter  vane.  Surface  pressure  data  was  measured  for  the  main  blade  and  splitter*. 


' Unfortunately  only  5  taps  were  used  to  measure  splitter  pressures. 


The  blade  definition  was  obtained  from  Reference  [5],  using  circular  leading  edges  on 
both  main  blade  and  splitter.  The  main  blade  and  splitter  thickness/chord  ratios  were 
3.61%  and  4.5%,  with  thick  trailing  edges,  1.964%  and  1.021%  respectively.  As  a  result 
of  problems  early  in  the  MISES  development,  a  slightly  elliptical  leading  edge  was  used 
for  the  main  blade  to  reduce  expansions  at  the  leading  edge.  These  problems  were  later 
corrected  so  that  circular  leading  edges  may  be  used,  but  the  comparisons  that  follow 
use  the  modified  blade  profile. 

4.4.1  Test  Flow  Conditions 

Flow  conditions  for  the  ARL  tests  were  carefully  set  up  and  thoroughly  documented, 
simplifying  the  comparison  of  computational  results  to  experiment.  An  upstream  wedge 
was  used  to  establish  periodic  flow  in  the  cascade  and  inflow  conditions  were  measured 
in  the  periodic  flow  field.  A  cross-section  view  of  the  test  section,  showing  the  blade  row 
between  the  strongly  contracted  side  walls  is  shown  in  Figure  4.15.  Suction  was  used  to 
remove  sidewall  boundary  layers  immediately  upstream  of  the  contraction.  The  blade 
aspect  ratio  was  approximately  1:1,  resulting  from  structural  limitations  on  the  blading. 
The  two-dimensionality  of  the  flow  over  such  a  low  aspect  ratio  blade  is  questionable. 

The  inlet  flow  conditions  were  an  inlet  Mach  number  of  1.46,  =  66.84°,  and  an 

inlet  Reynolds  number  of  1.45  x  10(’.  Turbulent  flow  conditions  were  assumed,  with 
transition  at  2%  chord  for  all  surfaces. 

Measured  data  from  the  test  was  obtained  by  flow  survey  at  27%  of  axial  chord 
downstream  of  the  trailing  edge.  Mass-averaged  near-field 1  exit  data  for  exit  pressure, 
angle,  loss  and  other  flow  properties  were  generated  from  this  data.  A  mixing  loss 
analysis  of  the  mass-averaged  data  was  used  to  generate  far-field  exit  data. 

'My  terminology,  near- field  refers  to  mass-averaged  exit  data  from  the  flow  survey,  far-field  refers  to 
mixed-out  exit  data 
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Figure  4.15:  Top  view  of  ARL  test  section  showing  sidewall  contraction 
4.4.2  Computational  Model 

Due  to  the  high  inlet  angle  J,  =  67'  which  accentuates  skew  problems  for  standard 
grids,  an  offset-periodic  grid  was  used  to  model  the  ARL  cascade.  Most  of  the  MISES 
results  were  obtained  using  a  220x20  computational  grid  with  2  offset  blocks  in  the  inlet 
region,  see  Figure  4.16.  This  grid  is  from  a  converged  solution,  showing  large  blockage 
effects  from  separated  boundary  layers  on  the  main  blade  and  splitter  vane.  Grids  with 
densities  from  180x16  to  240x32  were  used  with  essentially  identical  results  -  within 
2-3%  of  loss  and  turning. 

The  metal-metal  streamtube  width  for  this  case  is  shown  in  Figure  4.17,  correspond¬ 
ing  to  AVDR  =  2.053.  Also  shown  is  the  reduced  streamtube  width  corresponding  to 
AVDR  -  2.242,  or  109%  of  the  metal-metal  AVDR,  used  for  much  of  the  comparisons. 
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Figure  4.16:  Offset-periodic  220x20  grid  for  ARL  cases 


The  boundary  conditions  used  for  MISES  specified  exit  pressure  and  applied  an 
extra  DOF  using  c, „/  to  drive  the  inlet  Mach  number  to  1.46.  An  isentropic  region  was 
used  from  behind  the  leading  edge  bow  shock  through  the  leading  edge  expansion  to 
minimize  spurious  losses  due  to  dissipation  in  the  near-surface  streamlines. 


4.4.3  AVDR  Effects 

As  discussed  in  Section  4.1  AVDR  has  a  strong  effect  on  exit  angle  and  loss  for  cascades. 
These  effects  are  shown  in  Figures  4.18  and  4.19  for  the  ARL  cascade  with  50%  splitter 
vane  position  and  a  static  pressure  ratio  ps/pi  =  1.93.  The  intent  of  this  comparison 
was  to  choose  an  AVDR  that  best  matched  exit  angles  in  the  test  data.  Unfortunately 
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Figure  4.17:  Streamtube  width  for  ARL  cascade  at  metal-metal  and  109Tc  AVDR. 

a  single  AVDR  value  did  not  match  the  observed  variation  of  far-field  exit  angles,  in 
the  range  28'  -  312.  Rather  than  chase  the  ‘correct'  AVDR  for  each  pressure  ratio,  a 
compromise  AVDR  of  2.242  (109°c)  was  chosen  as  a  best  fit  at  pressure  ratios  of  1.8- 1.9. 
and  was  used  for  further  comparisons  and  design  studies.  The  sidewall  profile  for  the 
modified  AVDR  was  obtained  by  scaling  the  streamtube  thickness  about  the  inlet  value. 
This  is  not  justifiable  on  the  basis  of  using  the  viscous  blockage,  i.e.  boundary  layer 
thickness,  to  modify  the  AVDR.  The  issue  of  the  'real'  blockage  effects  is  considerablv 
more  complex  and  may  invalidate  the  quasi-3D  flow  assumptions.  This  will  be  discussed 
further  in  the  next  section. 

Pressure  distributions  for  the  two  extremes  are  shown  -  for  low  AVDR,  in  Figure  4.20. 
and  for  high  AVDR,  in  Figure  4.21.  These  show  that  the  splitter  loading  increases  at 
low  AVDR  from  the  locally  increased  incidence  at  the  splitter  due  to  separation  on  the 
main  blade.  Note  that  the  pressure  levels  on  the  lower  surface  of  the  main  blade  are 
strongly  affected  by  the  AVDR. 
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Figure  4.20:  Pressure  distribution  for  AVDR  =  1.95. 


Figure  4.21:  Pressure  distribution  for  AVDR  =  2.295. 
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4.4.4  Pressure  Comparison 


The  difficulty  in  comparing  pressure  data  for  supersonic  cascades  lies  in  determining  the 
features  that  should  be  matched  using  the  available  parameters.  Since  MISES  results 
show  strong  curvature  of  exit  streamlines  immediately  behind  the  blade,  see  Figure  4.16, 
exit  angles  from  MISES  were  compared  to  the  experimental  far-field  exit  angle  rather 
than  the  near-held  angle.  Exit  angles  in  MISES  were  measured  at  the  exit  boundary 
one  axial  chord  downstream  and  represent  almost  completely  mixed-out  values.  For  the 
ARL  case,  AVDR  was  adjusted  to  109%  of  the  metal-metal  AVDR  to  roughly  match 
exit  angle  with  the  far-field  exit  angle  from  the  test.  This  corresponds  to  reducing  the 
streamtube  width  by  approximately  one-half  of  the  boundary  layer  thickness  at  the 
trailing  edge.  Since  the  sidewall  boundary  layers  and  blade  boundary  layers  must  be 
roughly  comparable  in  thickness,  this  correction  is  certainly  not  excessive.  In  addition, 
the  assumed  profile  for  streamtube  width  does  not  take  into  account  the  actual  axial 
variation  of  the  sidewall  boundary  layer  thickness. 

I  he  static  pressure  ratio  for  MISES  that  best  matched  surface  pressures  at  the 
trailing  edge  and  splitter  lower  surface  turned  out  to  be  the  same  as  the  far-field  static 
pressure  ratio  from  the  test.  This  indicates  that  the  mixed-out  flow  quantities  are 
more  appropriate  candidates  for  comparison  than  the  near-field  data.  Results  for  static 
pressure  ratio  p>jp\  =  2.034  are  shown  in  Figure  4.22.  displaying  excellent  agreement  to 
test  data  from  Reference  (5 j .  The  near-field  pressure  ratio  in  the  test  v/a.sp>fpt  =  2.011. 
Mach  contours  for  this  flow  condition  are  shown  in  the  companion  Figure  4.23. 

Results  for  static  pressure  ratio  pilp\  =  1.741  are  shown  in  Figure  4.24.  The  agree¬ 
ment  on  the  splitter  is  good,  showing  a  strengthening  of  the  upper  passage  shock.  The 
results  for  the  main  blade  show  considerable  mismatch  on  the  lower  surface  from  mid¬ 
blade  to  the  trailing  edge.  The  measured  near-field  static  pressure  ratio  for  this  case 
was  pijp\  =  1.592.  Mach  contours,  in  Figure  4.25,  show  passage  shocks  above  and 
below  the  splitter  vane. 
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The  effect  of  static  pressure  ratio  on  exit  angle  is  shown  in  Figure  4.26.  Computa¬ 
tional  results  for  exit  angle  are  more  nearly  matched  to  far-field  exit  angles  (no  surprise 
here)  and  show  a  similar  trend  with  pressure  ratio.  The  effect  of  static  pressure  ra¬ 
tio  on  loss  is  summarized  in  Figure  4.27  showing  that  the  computed  loss  levels  were 
25-40%  lower  than  measured  far-field  loss.  The  computed  and  test  results  also  show 
different  trends  with  exit  pressure,  with  the  test  data  indicating  decreasing  losses  at 
lower  pressure  ratios,  a  somewhat  unexpected  result. 
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Figure  4.26:  Effect  of  pressure  ratio  on  exit  angle  for  ARL  cascade  109%  AVDR. 

The  comparison  of  MISES  results  and  test  data  shows  that  there  may  be  addi¬ 
tional  factors  at  work  not  accounted  for  in  the  analytical  modeling  or  the  experimental 
measurements.  The  predicted  and  measured  loss  levels  show  an  uncomfortably  large 
mismatch  and  they  show  different  trends  with  pressure  ratio.  The  surface  pressure 
comparisons  showed  an  excellent  match  at  p-i/p ]  =  2.034,  yet  the  comparison  at  lower 
pressure  ratio  p-j/pi  =  1.741  showed  significant  differences  in  the  pressure  distributions 
on  the  main  blade  lower  surface  near  the  trailing  edge.  Considering  the  boundary  layer 
thickness  variations  and  the  sensitivity  of  loss  and  exit  angle  to  AVDR  it  is  not  surpris¬ 
ing  that  there  are  differences  between  these  results,  especially  because  of  the  simplistic 
modeling  of  cascade  blockage  with  AVDR.  For  the  analytical  model  AVDR  was  adjusted 
solely  to  match  far-field  exit  angle  which  is  largely  insensitive  to  details  in  local  blockage. 
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Figure  4.27:  Effect  of  pressure  ratio  on  loss  for  ARL  cascade  109%  AVDR. 

Since  the  blade  surface  boundary  layers  have  drastically  different  thicknesses  and  the 
sidewall  BL  thicknesses  must  be  comparable  to  those  on  the  blade,  the  actual  blockage 
effects  are  very  different  for  the  upper  and  lower  passages.  The  differences  observed  in 
surface  pressures  are  likely  to  result  from  the  local  differences  in  sidewall  blockage. 

The  losses  measured  in  the  tests  may  also  be  artificially  high  as  a  result  of  the 
nearly  stagnant  sidewall  boundary  layer  fluid  being  driven  by  the  blade-to-blade  pressure 
gradient  -V/>  into  the  suction  surface  boundary  layer,  thickening  it  and  giving  an 
incorrectly  high  loss.  This  flow  mechanism  is  pictured  in  Figure  4.28,  showing  a  view  of 
the  flow  passage  with  paths  for  the  inviscid  and  wall  streamlines.  One  cannot  assume 
that  the  upstream  boundary  layer  bleed  will  eliminate  sidewall  boundary  layer  problems, 
since  the  separated  blade  and  sidewall  boundary  layer  thicknesses  can  reach  10%  of  the 
passage  height. 
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4.5  Convergence  of  Newton  Solution 


The  convergence  of  the  Newton  system  is  largely  determined  by  the  extent  of  strong 
compressible  effects  and  viscous  coupling,  both  of  which  become  progressively  more 
significant  from  the  UTRC  subsonic  case  to  the  supersonic  ARL  cascade  with  its  strong 
viscous  blockage  effects.  Convergence  histories  for  the  analysis  test  cases  are  shown  in 
Figures  4.29  and  4.30.  These  figures  show  the  maximum  density  change  and  loss  level 
for  each  Newton  cycle,  and  demonstrate  the  strong  influence  of  Mach  number  on  the 
convergence  rate.  It  is  characteristic  of  the  Newton  solution  that  the  terminal  rate  of 
convergence  to  machine  accuracy  is  roughly  the  same  for  all  cases,  but  the  non-linear 
cases  spend  many  cycles  moving  non-linear  residual  errors  around  before  getting  to  that 
point. 

The  maximum  density  change  is  normally  used  to  indicate  convergence  for  MISES, 
with  a  10-J  drop  as  the  threshold  for  engineering  accuracy.  Further  convergence  is  rec¬ 
ommended  for  parameter  sensitivity  analysis  (e.g.  loss  sensitivity,  see  C.4),  as  functional 
gradients  are  inherently  much  noisier  than  the  functions  themselves. 

The  actual  speed  of  the  method  depends  on  the  time  required  for  each  Newton  cycle, 
which  is  a  function  of  the  grid  topology  and  density.  For  example,  on  a  DEC  DS5000 
workstation  (approx.  25MIPS,  4MFLOPS),  the  UTRC  test  case  took  9  sec.  per  Newton 
cycle  on  a  132  x  20  standard  grid,  or  1  minute  for  a  converged  viscous  solution.  The 
ARL  case  took  roughly  40  sec.  per  Newton  cycle  on  a  220  x  20  offset-periodic  grid,  or 
roughly  20  minutes  for  a  converged  solution.  By  way  of  contrast  with  other  methods, 
the  ARL  case  was  also  run  with  UNSFLO  [19],  a  time-marching  Euler  code  with  an 
embedded  Navier-Stokes  viscous  layer.  The  case  never  fully  converged,  but  required 
over  5000  iterations  and  many  hours  of  computer  time  before  the  flow  transients  began 
to  settle  out. 

Convergence  from  one  flow  condition  to  another,  ‘near-by’  flow  condition  depends 
on  the  non-linear  phenomena  in  the  flow  field.  Typically,  subsonic  cases  require  only 
3-4  Newton  cycles  to  reconverge  for  near-by  conditions,  while  6-8  cycles  were  required 
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for  the  supersonic  AEL  case.  Shock  movement,  or  changes  in  the  point  of  transition 
or  separation  slow  convergence  until  the  final  position  is  attained.  Note  that  small 
variations  in  specified  flow  conditions  may  not  produce  near-by  solutions.  Sometimes  a 
particularly  difficult  flow  condition  must  be  approached  through  a  succession  of  near-by 
solutions,  or  in  extreme  cases  it  may  be  necessary  to  reconverge  a  case  from  scratch. 

The  relatively  large  number  of  cycles  to  convergence  for  the  ARL  cases  can  be  re¬ 
duced  by  providing  a  better  initial  condition  for  the  flow  field.  The  ‘standard’  initial 
condition  for  cascades  uses  a  smooth  gradient  of  density  along  the  passage  specified  by 
inlet  conditions  and  exit  pressure.  Boundary  layers  are  initialized  to  zero  S'.  A  signifi¬ 
cantly  better  initial  condition  can  be  obtained  by  interpolating  density  and  displacement 
thickness  from  a  previous  solution  onto  a  new  grid1.  This  reduces  the  number  of  cycles 
to  convergence  from  ~  32  to  10,  as  shown  in  Figure  4.31. 

The  best  convergence  for  the  ARL  case  was  obtained  by  specifying  an  exit  pressure 
that  produced  a  started  condition,  with  the  shock  in  the  passage,  but  not  so  low  as  to 
strongly  choke  the  splitter  passages,  typically  pijp \  =  2.03  was  used.  Normally,  MISES 
will  converge  regardless,  except  when  there  may  be  no  steady  solution.  Unstarted 
solutions  for  the  ARL  splitter  at  AVDR.  =  2.242  and  p>!p\  >  2.06  never  converged, 
suggesting  that  only  unsteady  solutions  may  exist  for  that  state'. 


'The  initialled  displacement  thickness  probably  has  the  more  significant  effect  for  the  ARL  case 
with  its  extensive  separation. 

'Although  the  nominal  ARL  cascade  had  no  steady  unstarted  solution,  several  of  the  redesigned 
blades  were  found  to  have  such  solutions. 
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Figure  4.31:  Convergence  of  Max  density  change  for  ARL  Splitter  with  initialized  flow 

field 
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Chapter  5 

Results  -  Design  Applications 


The  previous  chapter  compared  analytical  results  with  experimental  data  for  several 
test  cases,  including  the  supersonic  ARL  cascade  with  splitter  vane  at  its  nominal  50% 
position.  This  chapter  will  focus  on  optimization  of  the  ARL  cascade  with  splitter  vane. 
The  principal  objectives  are 

•  Compare  predictions  for  optimal  splitter  circumferential  position  with  test  data. 

•  Redesign  main  blade  and  splitter  for  improved  performance. 

•  Examine  design  space  for  profitable  directions  for  design  changes. 


5.1  ARL  Cascade  w/o  Splitter 


Before  examining  the  splitter  vane  in  detail,  the  performance  of  the  basic  ARL  cascade 
without  splitter  vane  was  established.  Essentially  the  same  computational  model  de¬ 
veloped  for  the  ARL  cascade  in  the  previous  chapter  was  used,  with  the  splitter  blade 
removed.  The  same  inflow  and  boundary  conditions  were  used. 

Results  at  static  pressure  ratios  of  1.93  and  2.03  are  plotted  with  data  from  Refer¬ 
ence  [2]  in  Figures  5.1  and  5.2.  The  pi/p\  =  1-93  results  are  compared  to  test  data  at 
a  near-field  pressure  ratio  of  1.82.  Corresponding  pressures  and  forces  for  the  cascade 
with  splitter  vane  can  be  compared  between  Figure  5.11  and  Figure  4.22  in  the  previous 
chapter.  Analytical  results  indicate  that  the  addition  of  the  splitter  increased  turning 
by  roughly  3°  for  a  30%  increase  in  loss. 
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Figure  5.1;  Pressure  distribution  for  main  blade  at  p2/p,  =  1.93,  109%  AVDR 


Figure  5.2:  Pressure  distribution  for  main  blade  at  p2/pi  =  2.03,  109%  AVDR 
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5.2  Splitter  Vane  -  Background 


The  original  splitter  vane  was  added  to  the  ARL  cascade  in  an  attempt  to  reduce 
deviation.  It  was  designed  to  extend  between  50%  and  100%  axial  chord  and  used 
the  mean  camber  line  of  the  main  blade  in  the  overlap  region.  The  vane  thickness 
distribution  corresponded  to  that  of  a  circular-arc  blade  with  the  same  inlet  and  exit 
angles  and  same  maximum  thickness  as  the  main  blade.  The  splitter  thickness/chord 
ratio  was  4.5%.  The  nominal  position  for  the  vane  was  midway  between  the  main  blades. 


5.3  Splitter  Vane  Circumferential  Position 


The  effects  of  splitter  circumferential  position  were  investigated  using  the  computational 
model  for  the  ARL  cascade  developed  in  the  previous  chapter.  Data  from  an  ARL 
experimental  study  [5]  for  splitter  vane  position  was  used  for  comparison  of  surface 
pressures,  turning  and  loss.  The  splitter  was  displaced  using  the  blade  movement  DOF 
D„„„ ..  and  specified  mode  constraint  discussed  in  Chapter  3.  Other  boundary  conditions 
were  identical  to  those  used  in  section  4.4.2. 

The  effect  of  splitter  circumferential  position  on  loss  is  summarized  in  Figure  5.3  for 
splitter  trailing  edge  positions  from  44%  to  56%.  The  computed  loss  levels  were  20-30% 
lower  than  the  measured  near-field  loss,  and  show  no  significant  variation  of  loss  levels 
with  splitter  position1. 

Pressure  distributions  for  the  44%  and  56%  splitter  positions  are  shown  in  Fig¬ 
ure  5.4  and  5.5,  compared  to  test  data  at  43%  and  55%  positions,  respectively.  The 
effect  of  upward  movement  of  the  vane  was  to  increase  splitter  loading.  The  effect  on 
the  lower  surface  of  the  main  blade  was  similar  to  the  low  pressure  ratio  comparison 
( pifp\  -  1.741)  in  the  previous  chapter,  where  the  MISES  result  indicates  a  strong 
acceleration  forward  of  the  trailing  edge.  Analytical  results  indicate  that  the  load  dis- 

1  This  conclusion  was  also  reflected  in  calculations  of  loss  sensitivity  to  circumferential  movement 
(see  section  C.4)  where  typical  values  of  dw/dDmov  =  ±0.001/%  were  obtained. 
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Figure  5.3:  Effect  of  splitter  vane  circumferential  position  on  loss 

tribution  between  main  blade  and  splitter  is  changed  by  splitter  circumferential  position, 
but  the  net  loading  and  turning  of  the  combination  is  unaffected.  Given  the  strong  side- 
wall  influence,  it  is  likely  that  the  variation  of  u>  in  the  test  is  due  to  sidewall  boundary 
layer  changes  not  modeled  with  the  fixed  AVDR  assumption.  The  precise  mechanism 
for  the  loss  variation  is  not  identifiable. 
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Figure  5.4:  Pressure  distribution  for  ARL  cascade  with  44%  splitter  position 


Figure  5.5:  Pressure  distribution  for  ARL  cascade  with  56%  splitter  position 
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5.4  Overview  of  ARL  Cascade  Flow  Field 


Before  proceeding  with  changes  to  the  blade  geometry,  we  first  review  the  problems  we 
intend  to  cure  by  redesign.  Analytical  results  have  shown 

•  The  flowfield  is  characterized  by  relatively  mild  forward  passage  shocks.  The 
inviscid  loss  accounts  for  30-60%  of  the  total  loss,  depending  on  pressure  ratio. 

•  The  forward  portion  of  the  splitter  is  not  loaded.  The  splitter  leading  edge  near 
the  sonic  point,  so  that  shocks  are  possible  in  either  or  both  passages. 

•  The  main  blade  pressure  surface  and  the  splitter  upper  surface  form  a  supersonic 
nozzle  terminated  by  a  strong  upper  passage  shock,  usually  with  an  associated 
separation  region  on  the  splitter. 

•  Strong  separation  is  present  on  the  main  blade  with  subsequent  reduced  turning. 

The  goals  for  any  optimization  or  redesign  are  generally  to  reduce  losses  and/or  increase 
turning.  The  trade-off  between  loss  and  turning  is  not  well  defined  here,  and  usually 
depends  on  details  of  the  cascade  application.  Profitable  directions  for  redesign  of  the 
ARL  cascade,  based  on  the  flow  features,  appear  to  be 

•  Reduction  of  the  strong  acceleration  in  the  upper  passage  to  reduce  shock  strength 
and  separation  on  the  splitter. 

•  Shifting  some  of  the  splitter  loading  forward  towards  the  leading  edge. 

•  Reducing  the  strength  of  the  main  blade  separation. 
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5.5  Modal  Re-design  of  ARL  Cascade 


The  splitter  vane  and  main  blade  were  redesigned  using  the  modal  design  option  dis¬ 
cussed  in  Chapter  3.  Modal  design  was  selected  because  of  the  presence  of  separated 
flow  on  both  main  blade  and  splitter.  Camber  redesign  was  used  to  avoid  thinning  the 
already  very  thin  blades. 

The  starting  point  was  the  ARL  cascade  with  the  original  splitter  at  p-,/p ,  ~  1.93 
and  109%  AVDR.  A  series  of  modal  design  steps  was  used,  explicitly  specifying  modal 
changes  while  monitoring  changes  in  loss,  turning  and  loss  sensitivity.  Mode  shapes 
called  out  in  Table  5.5  correspond  to  those  shown  in  Figure  3.6. 


Iteration 

_  _ 

Modes 

Blade 

Surface 

Design  region 

Purpose 

1 

1,  3,  5 

splitter 

both 

0-100% 

decamber 

2 

Brnov 

splitter 

both 

move  down 

3 

Drol 

splitter 

both  j 

increase  incidence 

4 

7 

splitter 

both 

70-100% 

add  rear  camber 

5 

i 

main 

lower  j  50-100% 

fill  out  lower  surface  ! 

6 

7 

main 

both  j  75-100% 

add  rear  camber 

Table  5.1:  Stepwise  modal  redesign  of  ARL  cascade  with  splitter  vane. 

The  intent  in  iteration  2  was  to  reduce  the  ‘nozzle’  effect  in  the  upper  passage.  Step 
5  pushed  out  the  main  blade  lower  surface  to  send  compression  waves  onto  the  splitter 
suction  surface  to  reduce  the  Mach  numbei.  Finally,  iteration  6  added  rear  camber  to 
the  main  blade  to  restore  the  turning  to  the  original  level.  These  changes  could  also 
have  been  accomplished  using  the  inverse  design  mode,  but  the  technique  of  explicit 
stepwise  reshaping  combined  with  loss  sensitivity  analysis  gave  more  insight  into  the 
effects  of  each  perturbation.  The  strong  effect  of  turning  on  loss  presents  a  problem  for 
design  guided  by  sensitivity  analysis  because  a  steepest  descent  down  the  loss  sensitivity 
gradient  in  modal  design  space  also  tends  to  reduce  turning.  A  more  automated  design 
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scheme  should  incorporate  DOF’s  and  constraints  to  evaluate  sensitivity  with  turning 
held  constant. 

The  blade  shapes  for  the  modified,  MOD1C,  cascade  are  compared  to  the  original 
blades  (dashed  lines)  in  Figure  5.6.  The  pressure  distribution  for  the  modified  cascade  is 


Figure  5.6:  Redesigned  blade  section  compared  to  original  ARL  blades. 

shown  in  Figure  5.7.  For  comparison,  the  pressure  distribution  for  the  original  blading  is 
shown  in  Figure  5.10.  The  splitter  loading  has  been  moved  forward,  but  loading  on  the 
main  blade  has  decreased.  The  design  changes  resulted  in  0.32°  lower  turning  for  a  2.7% 
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reduction  in  loss,  a  disappointingly  small  effect.  The  modified  cascade  was  run  for  the 
range  of  pressure  ratios  1.74-2.03  and  demonstrated  essentially  identical  performance 
to  the  original  cascade  as  shown  in  Figure  4.27.  This  insensitivity  to  design  changes  is 
similar  to  that  observed  for  turbulent  flow  airfoils  where,  so  long  as  the  flow  remains 
attached  and  there  are  no  strong  compressible  effects,  the  airfoil  drag  is  basically  fixed 
by  the  Reynolds  number2.  Since  the  transition  was  fixed  for  both  blades  in  the  ARL 
cascade,  a  similar  effect  may  be  operating  for  cascades.  Still,  considering  the  shock  on 
the  upper  surface  of  the  splitter,  this  degree  of  insensitivity  is  unexpected. 


Figure  5.7:  Surface  pressures  for  redesigned  ARL  cascade  at  p^/pi  =  1-93 


‘The  possibilities  for  optimisation  are  expanded  considerably  for  laminar  flow  airfoils,  however. 
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5.6  Alternate  Design  Directions 


The  modal  redesign  resulted  in  little  improvement  in  cascade  performance.  This,  com¬ 
bined  with  the  small  effect  of  splitter  circumferential  position,  implies  that  the  design 
space  should  be  explored  in  other  directions.  This  section  examines  the  effects  of  splitter 
rotation,  splitter  sizing  and  axial  position  on  cascade  performance. 

5.6.1  Splitter  Vane  Rotation 

The  splitter  position  study  showed  that  the  circumferential  position  changed  the  relative 
loading  of  main  and  splitter  blades,  but  had  almost  no  effect  on  net  turning  and  loss.  The 
effect  of  splitter  incidence  on  loading  was  investigated  for  a  fixed  splitter  trailing  edge 
position  (50%).  The  splitter  vane  was  rotated  from  -3°  to  +3°  using  the  blade  rotation 
DOF  and  fixed  mode  constraint.  Computed  pressure  distributions  for  p>/p\  =  1.93  and 
109%  AVDR  are  compared  in  Figure  5.8.  The  loss  and  exit  angles  for  splitter  rotation 
are  summarized  in  Table  5.2,  showing  an  increase  in  both  turning  and  loss  for  positive 
incidence  on  the  splitter.  Note  that  the  -3°  rotation  caused  the  lower  passage  to  form 
a  diverging  nozzle  with  a  mild  shock  on  the  splitter  lower  surface. 


Incidence 

w 

02 

+  3.0 

0.1422 

29.38° 

0.0 

0.1274 

30.26° 

-3.0 

_ 

0.1300 

_ 

30.58° 

Table  5.2:  Effect  of  splitter  rotation,  pz/pi  =  1.93,  109%  AVDR. 
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Figure  5.8:  Surface  pressures  for  ARL  cascade  at  ±3°  splitter  incidence 
5.6.2  Short  Chord  Splitter  Vane 

The  results  for  splitter  position  and  rotation  imply  a  strong  interaction  between  main 
blade  and  splitter  loading.  Since  the  effect  of  the  splitter  is  to  cancel  main  blade  loading 
in  the  overlap  region,  a  shortened  splitter  was  investigated  as  a  means  to  increase  loading 
on  the  main  blade.  The  splitter  was  shortened  to  70%  of  its  chord  with  the  same 
relative  camber  and  trailing  edge  at  the  nominal  (50%)  position.  The  relative  thickness 
was  increased  to  maintain  the  same  blockage  as  the  original  splitter.  The  pressure 
distributions  for  the  original  splitter  and  shortened  splitter  are  compared  in  Figures  5.9 
and  5.10.  The  shortened  splitter  gave  basically  the  same  performance  as  the  original 
splitter  —  a  slight  decrease  in  turning  for  a  slight  decrease  in  loss  (probably  due  to  the 
reduced  skin  friction). 
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Figure  5.9:  Surface  pressures  for  shortened  chord  splitter  at  p2/p t  = 


Figure  5.10:  Surface  pressures  for  original  splitter  at  P2/P1 


5.6.3  Splitter  Vane  Axial  Position 

Another  means  to  reduce  the  interaction  between  the  main  blade  and  splitter  is  to 
stagger  them  axially  by  moving  the  whole  splitter  blade  downstream.  The  original 
splitter  vane  was  used,  moved  and  rotated  along  a  slight  arc  to  maintain  a  position 
midway  between  the  extended  camber  lines  of  the  main  blades.  Since  a  strong  effect  on 
turning  was  found  for  the  downstream  positions,  a  higher  static  pressure  ratio  p?/Pi  = 
2.03  was  used.  Table  5.3  summarizes  the  effects  of  splitter  axial  movement,  showing  a 
large  increase  in  turning  (+4°)  for  a  slight  increase  in  loss  levels.  The  higher  turning 


Axial  TE  Position 

IjJ 

02 

vane  incidence 

1.00  (nominal) 

0.1294 

29.19° 

O 

O 

o 

1.10 

0.1274 

28.66° 

3.5° 

1.15 

0.1321 

27.44° 

5.2° 

1.20 

0.1347 

25.18° 

7.0° 

Table  5.3:  Effect  of  splitter  axial  position,  p_>/Pi  =  2.03,  109%  AVDR. 

available  at  the  most-downstream  position  increased  the  pressure  ratio  obtainable  before 
unstarting  from  pi>/pi  =  2.06  for  the  nominal  position  to  Pu'P\  =  2.23.  The  increased 
turning  results  from  higher  main  blade  loading,  as  shown  in  Figure  5.12.  For  comparison, 
Figure  5.11  shows  the  pressure  distribution  on  the  original  cascade  at  p-_» / p i  =  2.03. 
Attempts  to  increase  the  turning  further  by  redesigning  the  splitter  for  higher  loading 
were  unsuccessful  -  increased  forward  loading  on  the  splitter  simply  reduced  loading  on 
the  main  blade  at  the  expense  of  loss. 
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Figure  5.11:  Surface  pressures  for  AEL  cascade  at  p2/pi  =  2.03 


Figure  5.12:  Surface  pressures  for  ARL  cascade  with  splitter  moved  +0.20  downstream 
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5.6.4  Tandem  Splitter  Vane 


The  high  turning  from  the  aft-most  splitter  position  indicates  that  a  different  approach 
to  splitter  design  may  be  required.  The  appropriate  model  for  splitter  design,  if  axial 
displacement  is  possible,  may  be  to  treat  the  splitter  as  the  aft  element  of  a  multi¬ 
element  airfoil.  Multi-element  airfoils  achieve  high  lift  without  flow  separation  on  the 
main  element  by  reducing  the  trailing  edge  pressure  with  the  suction  from  a  smaller  aft 
element.  The  second  element  can  then  be  used  to  recover  the  flow  to  the  trailing  edge 
pressure  because  its  boundary  layer  is  thin. 

A  tandem  configuration  was  defined  using  the  70%  chord  splitter  moved  aft  by  0.35 
and  up  by  0.35  toward  the  main  blade  pressure  surface,  and  rotated  by  13°  to  increase 
incidence  to  match  the  main  blade  camber.  The  tandem  blading  is  compared  to  the 
original  ARL  splitter  (dashed  lines)  in  Figure  5.13.  Some  care  was  required  to  generate 
the  grid  for  this  configuration  as  a  result  of  the  proximity  and  overlap  of  the  blades. 
Like  the  previous  case,  the  tandem  splitter  was  run  at  a  higher  static  pressure  ratio 
p>!p\  =  2.23  as  a  result  of  its  higher  turning. 

The  effect  of  the  tandem  configuration  was  surprisingly  large  -  resulting  in  both 
dramatically  increased  turning  (  +  12°)  and  a  15%  reduction  in  loss.  The  improvement 
in  turning  is  due  to  a  higher  main  blade  loading,  as  shown  in  Figure  5.14.  Mach 
contours  for  this  case  are  shown  in  Figure  5.15.  The  separated  wake  from  the  main 
blade  ‘bursts’  as  it  encounters  the  pressure  rise  downstream  of  the  splitter,  apparently 
without  any  significant  loss  penalty.  The  pressure  ratio  at  the  unstart  point  was 
raised  to  P2/P1  =  2.35.  Redesign  changes  made  to  the  splitter  to  minimize  the  shock  at 
the  trailing  edge  resulted  in  surprisingly  small,  hard- won  improvements  to  the  loss  and 
turning.  Again,  increased  forward  loading  on  the  splitter  acted  to  decrease  aft  loading 
on  the  main  blade  and  reduce  turning. 
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Figure  5.14:  Surface  pressures  for  70%  chord  tandem  splitter,  109%  AVDR. 


Figure  5.15:  Mach  contours  for  70%  chord  tandem  splitter 


Chapter  6 

Conclusions  and  Recommendations 


This  thesis  has  described  the  development  of  MISES,  a  new  computational  method 
for  application  to  cascade  analysis  and  design.  This  work  extends  the  fully  coupled 
viscous/inviscid  ISES  method  for  turbomachinery  application.  Quasi-three-dimensional 
effects  for  streamtube  width  and  radius,  as  well  as  for  rotational  effects  have  been 
incorporated  so  that  the  method  is  applicable  to  blade- to- blade  design  of  axial  stator 
and  rotor  stages. 

MISES  uses  a  zonal  approach  to  solve  the  coupled  Euler  and  integral  boundary  layer 
equations.  An  intrinsic  streamline  grid  is  used,  resulting  in  considerable  simplification 
of  the  inviscid  equations.  The  large  non-linear  system  of  equations  is  solved  with  a  full 
Newton-Raphson  method  using  a  specialized  direct  solver.  This  allows  the  addition  of 
extra  unkno%vns  and  constraints  that  enforce  specified  design  conditions.  The  modal 
design  capability  makes  possible  well-conditioned  redesign  of  blades  with  areas  of  flow 
separation. 

The  flow  solver  was  generalized  to  treat  cascades  with  splitter  vanes.  A  new  grid 
topology  using  an  offset-periodic  inlet  grid  was  introduced  for  supersonic  cascades  to 
circumvent  problems  with  grid  skew.  The  new  grid  also  has  the  benefit  of  improving 
resolution  of  shock  waves.  Loss  can  be  calculated  from  exit  conditions  or  from  a  uniform 
mixed-out  state.  Loss  sensitivities  to  cascade  parameters  such  as  Mach  number,  exit 
angle  or  exit  pressure  may  also  be  calculated  as  a  ‘free’  byproduct  of  the  Newton  solver. 
These  sensitivities  can  be  used  to  guide  design  optimization.  The  method  is  attractive 
for  design  applications  due  to  its  speed  and  robustness.  Solutions  are  obtained  in  6-30 
Newton  cycles  —  a  matter  of  1-20  minutes  on  a  fast  workstation. 
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6.1  Test  Cases 


Results  were  presented  to  validate  the  new  capabilities  of  the  method  on  subsonic, 
transonic  and  supersonic  cascades.  The  results  for  the  UTRC  subsonic  cascade  showed 
excellent  correlation  of  loss  with  experimental  data.  The  DFVLR  test  case  verified 
the  quasi-3D  capabilities  with  comparisons  of  surface  presssures  and  loss  for  a  blade 
at  several  AVDR  and  for  transonic  and  supersonic  Mach  numbers.  This  case  suffered 
from  poorly  specified  conditions  for  supersonic  inlet  conditions,  making  comparisons 
difficult.  The  results  showed  that  the  streamtube  width  effects  were  accurately  modeled 
and  losses  generally  compared  well  with  test  data. 

The  ARL  case  proved  to  be  the  most  difficult  test  of  the  quasi-3D  capability  as  the 
result  of  its  sensitivity  to  streamtube  contraction  and  viscous  blockage.  Discrepancies 
between  computational  results  and  experimental  data  are  undoubtedly  due  to  three- 
dimensional  effects  from  sidewall  boundary  layers.  It  is  naive  to  assume  that,  in  presence 
of  such  strongly  separated  flow,  the  upstream  boundary  layer  bleed  eliminated  sidewall 
boundary  layer  effects.  The  sidewall  boundary  layers  must  be  comparable  in  thickness 
to  the  blade  boundary  layers,  hence  the  upper  and  lower  passages  will  show  radically 
different  degrees  of  blockage.  This  violates  the  basic  assumptions  for  the  quasi-3D  flow, 
making  detailed  comparison  to  test  data  questionable. 


6.2  Splitter  Vane  Optimization 

The  analytical  results  for  the  ARL  cascade  indicate  that  the  splitter  circumferential 
position  has  almost  no  effect  on  loss  or  turning.  This  surprising  result  was  reinforced 
when  major  redesign  changes  were  made  to  the  splitter  vane  and  main  blade  to  reduce 
loss.  After  significant  geometric  changes  were  made  to  both  blades  to  optimize  the 
surface  pressure  distributions,  the  performance  was  almost  unchanged  from  the  original 
cascade. 

The  general  effect  of  the  splitter  is  to  reduce  the  loading  of  the  main  blade  in  the 
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overlap  region.  With  the  splitter  vane  located  within  the  main  passage  the  net  loading 
and  loss  appear  to  be  independent  of  the  distribution  of  loading  between  the  vane  and 
the  main  blade.  When  the  vane  is  moved  aft,  significant  re-loading  of  the  main  blade 
occurs  and  turning  is  increased.  The  best  position  for  the  splitter  vane  appears  to  be 
in  a  tandem  arrangement  similar  to  the  aft  element  in  a  multi-element  airfoil.  This 
configuration  resulted  in  +12°  of  turning  and  a  15%  reduction  in  loss. 


6.3  Recommendations  for  Further  Work 


There  are  several  extensions  that  should  be  considered  for  this  work: 

•  Validation  of  MISES  rotational  capability  using  analytical  test  cases  or  high  qual¬ 
ity  experimental  data. 

•  Redesign  to  reduce  loss  using  sensitivity  information  also  reduced  turning.  A 
more  automated  design  strategy  should  allow  the  designer  to  fix  important  design 
parameters  such  as  exit  angle.  Such  constraints  can  be  easily  implemented  with 
the  framework  for  the  current  method. 

•  The  entrainment  of  the  bow  shock  wake  into  the  boundary  layers  is  not  fully  rep¬ 
resented  in  the  present  viscous/inviscid  method.  Although  it  is  currently  not  clear 
how  significant  this  effect  is,  a  model  for  mixing  within  the  near-wall  streamtubes 
could  be  implemented. 

•  The  linearization  of  the  Jacobian  entries  using  the  chain  rule  is  laborious  and 
error-prone.  The  possibility  of  developing  Newton  me  hods  using  finite  differences 
to  calculate  Jacobian  entries  should  not  be  overlooked  [20].  These  approximate 
schemes  could  be  effectively  used  to  experiment  with  different  approaches  or  to 
check  analytic  derivatives  for  errors. 
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6.4  Characteristics  of  Ill-posed  Boundary  Conditions 


As  discussed  briefly  in  Chapter  2,  it  is  possible  to  specify  combinations  of  boundary 
conditions  that  are  inconsistent  with  respect  to  the  flow  physics.  With  experience  it 
is  possible  to  recognize  such  ill-posed  conditions  by  the  characteristic  behavior  of  the 
Newton  solver.  Figure  6.1  is  presented  as  a  guide  to  aid  in  recognizing  these  conditions. 
Good  Luck! 


Wbal-  Could  fee  dangerous  olsouV 
rurv\'mq  a  New4on  SoWtr  on  a 
^erioXc.  cascc.de  \  c.  cpp 

"'an . 


Figure  6.1:  Characteristic  behavior  for  ill-posed  boundary  conditions. 
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Appendix  A 

Development  of  Quasi-3D  Euler  Method 


The  development  of  the  quasi-3D  inviscid  equations  is  an  extension  of  the  work  by 
Giles  [15]  and  Drela  [9]  on  the  ISES  code.  This  chapter  derives  the  quasi-3D  formulation 
for  the  inviscid  solver,  and  draws  on  the  material  in  these  references. 


A.l  Coordinate  System 


The  three-dimensional  flow  problem  is  approximated  as  a  quasi-3D  blade- to-blade  flow 
defined  on  axisymmetric  stream  surfaces.  The  basic  geometry  of  the  problem  is  illus¬ 
trated  in  Figure  2.2.  Axisymmetric  stream  surfaces  are  defined  in  terms  of  the  radius 
r(z)  and  streamtube  thickness  6(z),  where  z  is  the  axial  coordinate,  see  Figure  2.3.  The 
axisymmetric  stream  surface  radius  and  thickness  distribution  can  be  obtained  by  a 
meridional  through-flow  method. 

The  blade- to-blade  flow  is  solved  in  a  local  m',  9  stream  surface  coordinate  system 

with  stream  surface  radius  and  thickness.  The  new  coordinate  m'  is  defined  as  the 

arc-length  in  the  r  -  z  plane,  normalized  by  r. 

dm  _  f  \f{dr)2  +  (dz)2 
r  J  r 

Note  that  m',  9  is  a  conformal  mapping  of  the  surface  of  revolution  to  a  plane,  preserving 
angles  and  shapes,  and  that  both  m'  and  9  are  dimensionless  and  defined  in  radians. 
Since  m'  is  an  integral  quantity  it  can  be  referenced  to  any  convenient  point,  such  as 
the  leading  edge.  The  relationship  between  m'  and  z  is  unique  so  that  r  and  b  may  be 
defined  r  =  r(m/)  and  b  =  b(m'). 

In  the  following  treatment  m\  9  will  be  replaced  by  x,  y  for  convenience,  except 


115 


where  noted.  The  coordinates  of  a  point  (z,  y)  in  the  dimensionless  system  transforms 
to  ( rx ,  ry)  in  the  physical  system.  The  distance  between  two  points  A  and  B  in  the 
physical  system  is  approximated  for  jr. v  -  fa)  <  r.i  +  ru  by 

\AB\  ss  r--~B  |(zB  -  *,»)*  +  {y b  -  y,i)il 

where  i  and  j  are  the  unit  vectors  in  the  m!  and  the  9  directions. 


A. 2  Euler  Equations 


As  discussed  in  Chapter  2  the  inviscid  flow  is  solved  using  the  steady-stace  Euler  equa¬ 
tions  in  a  blade-relative  coordinate  system  rotating  with  angular  velocity  ft.  The  finite 
volume  form  was  derived  for  conservation  of  mass,  momentum  and  energy  (equations  2.2 
through  2.4). 

The  discrete  equations  approximate  the  integral  form  of  the  Euler  equations  on 
an  intrinsic,  or  streamline,  grid  which  evolves  with  the  solution.  The  computational 
domain  is  made  up  of  conservation  cells  defined  on  a  set  of  streamtubes  that  make  up 
the  problem  domain.  The  continuity  and  energy  equations  become  simple  algebraic 
statements  of  constant  mass  flow 

6m  =  0 

and  constant  rothalpy  for  each  streamtube. 

61  =  0 


The  momentum  equation 

j>  p(q  •  n)q  +  ph  dA  =  -  p(  dV  with  f  =  ft  x  ^ft  x  rj  +  2ft  X  q 

r'M'  f 

is  also  simplified  since  the  conservation  cells  are  defined  such  that  only  one  pair  of 
opposing  faces  has  mass  flux  across  them,  the  other  faces  are  streamline  faces  and 
contribute  only  pressure  terms.  This  has  the  effect  that  velocity  and  density  are  only 
required  on  the  two  streamwise  faces,  simplifying  the  system  greatly. 
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Figure  A.l:  Conservation  cell  on  nodal  grid 

A. 3  Conservation  Cell 

The  flow  equations  are  solved  on  an  intrinsic  grid  that  aligns  in  one  of  its  directions 
with  the  streamline  direction.  Figure  A.l  shows  a  view  in  the  x,y  ( m',9 )  plane  the 
basic  arrangement  of  a  streamtube  conservation  cell  superimposed  on  its  defining  six 
grid  nodes.  The  lines  of  upper  +  and  lower  -  grid  nodes  define  streamwise  faces  of 
grid  which  have  no  mass  flux  across  them,  while  the  lines  connecting  the  corresponding 
pairs  of  upper  and  lower  nodes  define  the  quasi-normals  to  the  grid.  Another  view 
of  conservation  cell  is  shown  in  Figure  A. 2  illustrating  the  directed  face  area  vectors, 
consisting  of  streamwise  vectors  B+  and  B~  defined  between  midpoints  of  contiguous 
streamline  segments  and  pseudo-normal  area  vectors  A  \  and  .4;  defined  midway  between 
pairs  of  quasi-normal  grid  lines. 

The  faces  of  the  quasi-3D  conservation  cell  are  illustrated  in  Figure  A. 3,  where  the 
edge  vectors  A  and  B  define  face  areas  of  the  six-sided  conservation  volume.  The  edge 
faces  of  the  conservation  cell  are  denoted  by  1  and  2  for  the  streamtube  inlet  and  outlet 
faces  and  +  and  -  for  the  upper  and  lower  streamline  faces.  The  quasi-3D  cell  adds 
new  faces  on  the  inner  and  outer  stream  surfaces  of  revolution,  denoted  by  ()j  on  the 
inner  face  and  ()0  on  the  outer  face.  Note  that  the  areas  of  faces  At,  A2,  B+  and  B~ 
are  functions  of  b(x)  and  these  faces  lie  between  the  inner  and  outer  stream  surfaces 
whose  dimensional  coordinates  are  functions  of  r(z). 


117 


The  face  areas  are  defined  as  vectors  in  the  x  -  y  plane 


Ar,  =  r{b  | 


xj  +  x+  -X--  x2 


Ay.  =  ribi 


yt  +  vt  ~  V\  ~  V2 


B*  =  r+i*  *» -T  *t 
x  2 

s;  =  r-6-  ELiiL 

x  2 


^yj  =  f2&2 


y2+  +  ^3+  -  Vi  -  vt 


n+  =  r+6+  HLJL L 


-SlzIl 


where  the  face  radii  and  streamtube  thicknesses  are  given  by 


X  +  +  X+  +  X]  +x2 
r,  _r( - - - ) 


bi  =  b{^L^L) 


(A.l) 


*+  +  *J  +  x2-  +I3 

^  - 2 - ] 


x2+  +  x+  +  *2  +  xt  x 
62  ~6( - 2 - } 

*+  =  &(?*  +_£?) 
b-  =  H^) 


It  is  useful  to  define  a  streamwise  area  vector  5  and  a  quasi-normal  area  vector  jV 
from  the  averages  of  the  d  and  B  vectors 


B~  +  B+ 


.4 1  -I-  .4> 


The  grid  nodes  are  defined  to  lie  on  streamlines,  so  that  the  velocity  vectors  at  the 
center  of  the  1  and  2  faces  are  defined  lie  on  a  line  midway  between  the  upper  and  lower 
streamlines.  The  velocity  on  the  midpoint  of  the  left  cell  face  (at  point  1)  is  g*i  =  qi i|, 
with  the  density  and  pressure  at  that  point  p\  and  p\.  The  ii  vector  is  defined 


«  _  £ i  +  -  xt  -  *r 

2 


=  vjl  + 32 


yt  4-  y.2  -  yj  -  yf 
2 


i,  =  ^-i  +  ^7 
3 1 


(A.4) 


Similarly  the  flow  quantities  on  the  right  cell  face  2)  are  q-2  =  ?2i2,  Pi  and  p2,  where  the 
3-2  vector  is  defined 


_  x+  +  X-  -  x\  -  x£ 
r?  "  2 


_  yt  +  y3  -  yt  -  y2~ 
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J2  =  \/322  +  3l?  j2  =  (A.5) 

v  32  J2 

The  streamtube  cross-section  area  is  the  area  normal  to  the  the  a  face  velocity  vector. 
The  normal  areas  in  and  out  of  the  conservation  cell  are 

4n,  =  |i,  x  i’ll  An,  =  |i2  x  v42|  (A. 6) 

To  obtain  the  area  contributed  by  the  inner  and  outer  stream  surface  faces  consider 
the  normal  area  vector  of  the  inner  face  i  as  aj  =  ajnj,  similarly  for  the  outer  face  a0  and 
the  four  edge  faces.  The  normal  area  vectors  of  the  faces  of  the  the  closed  conservation 
cell  must  sum  to  zero,  by  simple  geometric  identity.  This  defines  the  net  normal  area 
vector  of  the  inner  and  outer  faces  to  be 


+  clo  —  — (di  -f  fflt  +  -fa) 


(A.7) 


The  normal  area  vectors  of  the  faces  are  expressed  in  terms  of  the  face  area  vectors  as 


d\  —  (  —  Ayit,  ATlj)  a-i  —  (Ay.,t,  Ar,j) 

a+  =  (-B^i,  Bti)  a~  =  (B~i,  - B~) ) 

The  volume  of  the  conservation  cell  is 

AV  = 

where  b  is  the  cell  average  streamtube  thickness. 


(A.8) 


(A.9) 


A. 4  Continuity  Equation 


The  mass  equation  is  simplified  by  use  of  the  streamline  grid  to  be  a  statement  of 
constant  mass  flow  everywhere  along  a  stream  tube.  For  each  cell,  therefore 

rh  =  P\  ?i  Ani  =  p2  ?2  An?  (A.10) 

The  mass  flow  in  each  streamtube  is  set  from  the  inlet  conditions  and  inlet  node  positions 
defining  the  streamtube  normal  area. 
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A. 5  Momentum  Equation 


The  momentum  equation  for  the  conservation  cell  is  composed  of  pressure  terms,  mass 
flux  terms  and  source  terms  due  to  rotation  (Coriolis  and  centrifugal  forces). 


F Jlux  Fprca*  —  F zonal is  4  Fc.r.nt 


(A.ll) 


The  discrete  equations  for  x  and  y  momentum  are  resolved  into  a  local  streamwise,  or 
S-momentum  equation  and  a  quasi-normal,  or  N-momentum  equation  by 


S- momentum 
N-momentum 


s'-Ef 


A. 5.1  Pressure  Terms 


The  pressures  on  the  inner  and  outer  faces  are  assumed  to  both  take  the  average  value 
of  the  pressures  on  the  four  edge  faces 

Pi  4  p>  -f-  n+  +  n 
Pi  =  P°  =  - 4 - 

In  this  case  the  sum  of  the  pressure  forces  acting  on  the  cell  faces  is 

-F,)r, =  P|d,  +  pidi  -)-  n+o+  +  n-a~ 

—  -  (pi  4  p?  4  n+  +  II  )  (ci|  4  (i'2  4-  5+  4a) 

4 


After  tedious  manipulation  this  becomes 

(Pi  ~  P2) 


Fprea<i  — 


_  (n+-n-).^ 

(a,  _  a2)  +  1 - - - L  (a+  -  a  ) 


4 


2  '  '  2 

(p,  4  Pi  -  n+  -  n-) . 


(ai  4  o2  -  a+  -  a  ) 


(A.  12) 


The  combined  pressure  term  has  special  significance  and  is  denoted  by 

„  n+  4  n-  -  Pi  -  P’i 

+  t r.nrr  — 


(A.  13) 
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which  will  approach  zero  in  the  limit  as  the  grid  is  refined,  where  the  averages  of 
the  pressures  on  opposite  face  pairs  will  be  equal.  Unfortunately,  Equation  A. 13  is 
transparent  to  an  odd-even  grid  mode.  To  prevent  grid  sawtooth  instabilities  another 
form  for  Pcorr  is  used,  based  on  the  isentropic  pressure-area  ratio 


Pro 


k  7PM2  j s'  x  s2~j  -  [sf  x  4| 


where  p  = 


\Sx  N\ 

and  Ml  = 


(A.14) 


P1+P2  ,  ,2  _  m*  +  Mi 


The  constant  k  is  used  to  set  the  relative  magnitude  of  the  Pcorr  term  and  typically  has 
a  value  0.1. 


The  x  and  y  components  of  the  pressure  forces,  using  the  area  normal  vectors  in  A. 8 
and  the  definitions  A. 3,  are 


-F, 


Xprr.%* 


=  -(pi-P2)Ny  -  -n~)sy 

_  p 

■»  rn  1 


•^vi  +  &>j 


-f,.,„  =  (Pl-p>)Nr  +  (n+-n-)5r 


+  Pro 


-ATi  +  ,4r,  +  B+  -  B~ 


(A.15) 


The  S  and  N  components  of  the  pressure  terms  are  formed  by  taking  the  dot  product 
of  A.15  with  the  5  and  N  vectors. 


~  S  ■  Ppresa  =  -  (Pt  “  P2)  1$  X  N] 

n  |5x(-A, +  I2  +  B+-5-)| 

_  rcorr  - 

-  n  ■  Fprr.„  =  -  (n-  -  n+)  15  x  n  1 

„  |iVx  (-A,  +i2  +  5+-F-)| 

—  Pcorr  „ 


(A.16) 
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Figure  A. 4:  Flux  terms  from  circumferential  displacement  of  velocity  vectors 
A. 5. 2  Flux  Terms 

The  mass  flux  terms  are  given  by  the  net  momentum  flux  through  cel!  faces  1  and 
2.  Additional  terms  are  required  for  the  quasi-3D  cell  because  the  x,y  coordinates 
are  defined  on  a  surface  of  revolution  and  the  unit  vectors  of  the  local  system  change 
direction  with  circumferential  coordinate  y  (or  9).  For  the  following  discussion  the  m'  -9 
terminology  will  be  used  in  place  of  x,y  for  clarity  with  regard  to  angles  on  the  curved 
stream  surface. 

The  basic  statement  for  the  mass  flux  forces  is 

-  Fjinx  =  m(A(gm.em<)  +  A(qgeg))  where  A()  =  ()2  -  ()i  (A. 17) 

where  m  =  pqA  is  the  mass  flux  in  the  streamtube  and  em<  and  eg  are  unit  vectors  on 
the  stream  surface.  Additional  terms  arise  for  the  A()  expressions  due  to  the  change 
in  directions  of  the  unit  vectors  with  a  displacement  AO  in  the  circumferential  direc¬ 
tion.  Figure  A. 4  illustrates  the  basic  geometry,  where  a  vector  {qm>,qg)  defined  on  a 
stream  surface  tilted  by  dr/(rdm')  from  the  axial  direction  is  displaced  by  A 6  in  the 
circumferential  direction. 
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The  Aqm>...  term  in  A.17  expands  to 

=  Aqm>emi  +  —  ) 

dr 

=  Agm'em>  +  qm'A9—^eo  (A.18) 

where  the  e#  term  is  due  to  the  radial  component  of  qmi  on  the  tilted  stream  surface, 
q„,>dr/(rdm'),  which  is  rotated  by  a  change  of  direction  A9,  producing  a  component 
normal  to  the  etni  vector  at  the  original  position. 


The  A qo—  term  in  A.17  expands  to 

A(qoeo)  =  Aqoeo  +  qo{*67  -  ) 


=  A qgtg  -  qgA9 


dr 

rdm‘ 


€ rn 1 


(A.  19) 


where  the  em/  term  is  due  to  rotation  of  qg  by  a  change  of  direction  A 9,  producing  a 
radial  component  -qoA9  that  has  a  component  along  the  tilted  stream  surface  in  the 
e,,,i  direction. 


The  force  due  to  the  mass  flux  is  then 

-  Ffinr  =  ™  ^(A(?r.i'  -  qo A9 —  )e,„.  +  (A(?fl  +  q,„> AS (A. 20) 
Continuity  for  the  cell  requires  that  qm'A9  -  qoAm'. 


The  eg  term  becomes 


&qo  +  qm'&9 


dr 

rdm' 


-  Aqg  +  qoAm 

&{rqo) 


A  r 
rAm' 


(A.21) 


The  em>  term  gives 

&qm'  +  qe&0 


dr 

rdm1 


=  A  qm<  +  q0  Am 


i  Ar 


rAm' 

&rq 


—  qgA9 


dr 

rdm' 


A (rqm>)  .  Ar  ,  dr 

- - - (^W?n,'Am  +qoA677^> 


rdm' 


~  ~~  -  {qm'Arn'  +  qeA9)^l 
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A. 5. 3  Rotational  Forces 

The  use  of  a  blade-relative  coordinate  system  brings  additional  acceleration  terms  into 
the  momentum  equation.  The  coordinate  system  is  rotating  at  angular  velocity  SI 
leading  to  Coriolis  forces  and  centrifugal  forces  on  the  fluid  in  the  cell. 
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Coriolis  Forces 


The  Coriolis  force  per  unit  volume  in  a  rotating  fluid  is 

%ii  =  -2,  (fi  X  f) 

Since  the  velocity  is  defined  in  the  S  direction  the  force  on  a  computational  cell  purely 
— * 

in  the  N  direction. 

N  ■  Fconoh 3  =  -2/3  (n  x  q'j  ■  N  AV 

=  -2pqQ,  ■  (s  x  N)  AV 

using  A„  =  Ji  x  iV|  and  m  =  pqA„ 
dr 

=  -2  mfi—  AV' 
rdx 

This  can  be  simplified  using  approximate  identities  for  dr /(rdx)  in  A. 25  and  the  defi¬ 
nition  of  Alr  to  be 


N  ■  Fcor.oi,,  =  -2mfl  {Nu{r  >  -  r,)  -  St/(r+  -  r  ))  (A.26) 


Centrifugal  Forces 


The  centrifugal  force  per  unit  volume  in  a  rotating  fluid  is  pCl2r  in  the  radial  direction. 
This  has  a  component  in  the  x  direction  due  to  the  inclination  of  the  stream  surface  by 
dr /(rdx)  from  the  axial  direction.  The  x  component  of  the  centrifugal  force  on  a  cell  is 

Fi“n*  =  pQ2rvh^v  (A,27) 

where  p  =  j(/>i  +  p 2),  the  average  cell  density. 


The  streamwise  and  normal  centrifugal  forces  on  a  cell  are 


S  •  Fteni  =  PCl2rSs-^-AV 
rdx 

>na 


=  (5  x  N\  using  A. 25 
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N-Fcenl  =  pCl2rNx-^AV 


using  A. 25 


(A. 28) 


A. 5. 4  S  and  N  Momentum  Equations 


In  the  preceeding  sections,  the  x  and  y  components  of  the  cell  forces  were  used  to 
define  momentum  components  in  both  a  “streamwise”  or  S-momentum  equation  and  a 
“normal”  or  N-momentum  equation.  Collecting  terms  from  Equations  A. 16,  A. 24,  A. 26 
and  A. 28  and  dividing  through  by  (S  x  ,/V|  we  obtain  the  S-momentum  equation 


-  (Pi  ~  Pi)  ~ 


where  the  new  terms  are  defined 


m  (-gi/i  +  q-ih) 

Prnrr  |S  X  (  — /T|  +  At  +  B+  —  B  )| 


|S  X  N\ 


|S  x  /Vj 


/,-■£• 5 


(A. 29) 


|5  x  iV| 


(A. 30) 


The  N-momentum  equation  is  similarly  defined  as 


j  -  2mft- 


IS  x 


2  m 
ri  +  r2 


r  —  r 

+  r2q2g2 - - - (?i/i  +  q2fi) 


-  (ir  -  n+)  - 


Pcorr  \Nx(-At  +A2  +  B+-B-)\ 

2  |5  x  1V| 


with  new  terms 


i[  •  N 
9  l  =  7JT 


|S  x  JV| 


i,  •  N 

g>  =  — - — 

IS  x  N | 


(A.31) 


(A. 32) 
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A. 5. 5  Upwinding  Scheme 


The  S  and  N  momentum  equations  are  modified  in  supersonic  zones  by  the  use  of 
artificial  dissipation  in  the  form  of  a  bulk  viscosity  term,  defined  as  in  Drela’s  thesis  [9], 
by  upwinding  the  speeds  and  q-z  in  A. 29  and  A. 31.  This  scheme  has  been  extended  to 
include  an  optional  second-order  dissipation  term  that  can  significantly  reduce  spurious 
total  pressure  losses  due  to  added  dissipation.  The  basis  for  this  idea  comes  from 
Giles’  second-order  dissipation  scheme  in  [15]  but  has  been  applied  to  the  bulk  viscosity 
dissipation  model.  The  form  of  the  upwinded  speed  is 

q,  =  q,  -  -M?.  -  ?,-i)  +  /*,(?.- 1  -  q,- 2)  (A. 33) 

The  upwinded  speed  thus  depends  on  the  speed  at  cell  face  i  and  one  or  two  upstream 
speeds,  depending  on  whether  second-order  dissipation  is  used. 

The  definition  for  the  first-order  upwinding  coefficient,  A  is 

A ,  =  fin)„  max  ^0,  — - — ^  where  M1  =  max(A/,J,  M,2_ , )  (A. 34) 

where  fir„„  is  a  scaling  factor  for  dissipation  (normally  1.0)  and  jl/,.r,/  is  a  threshold 
slightly  less  than  unity  (typically  0.9). 

The  second  order  upwinding  coefficient  is  set  such  that  fi  =  A,  producing  second- 
order  dissipation  when  the  grid  nodes  are  equally  spaced.  This  does  not  produce  enough 
dissipation  in  strong  gradients  (shocks)  and  p  is  reduced  by  a  factor  proportional  to 
velocity  gradient,  preferentially  for  decelerating  flow  to  minimize  spurious  loss  in  ac¬ 
celerations.  The  use  of  second  order  dissipation  complicates  the  numerical  form  of  the 
equations  by  introducing  dependence  on  the  two  previous  cell  face  speeds,  rather  than 
simply  the  previous  speed.  This  will  add  terms  to  the  matrix,  increasing  the  solution 
cost  slightly.  In  practice,  the  use  of  the  second  order  dissipation  has  been  limited  to 
transonic  Mach  numbers.  For  supersonic  inlet  Mach  numbers  the  second  order  scheme 
appears  to  have  some  unsolved  stability  problems. 


streamtube 


j-1  streamtube 


Figure  A. 5:  Pressures  on  neighboring  streamtubes 
A. 5. 6  Reduced  N-Momentum  Equation 

The  S-momentum  equation  A. 29,  using  the  upwinded  speeds  q |  and  q>,  depends  only 
on  the  geometry,  upwinded  speed,  density,  pressure  at  faces  1  and  2,  and  P,„rr.  Unfor¬ 
tunately  the  N-momentum  equation,  A. 31,  contains  additional  pressure  terms  Il+  and 
IT”  on  the  streamline  faces.  These  can  be  eliminated  to  reduce  the  number  of  unknowns 
and  form  the  reduced  N-momentum  equation  that  depends  only  on  P,,}rr  and  the 
geometry,  upwinded  speed,  density  and  pressure  at  faces  1  and  2.  This  is  done  in  the 
same  manner  as  in  Drela’s  thesis,  where  the  N-momentum  equation  A. 31  can  be  written 
in  the  form 

n+  -  IT  =  (•••) 

The  P,„rr  equation,  A.  13,  can  be  written  in  the  form 

n+  +  n  =  Pl  +  P’2  +  2  Pcorr 

The  sum  and  difference  of  these  two  equations  define  equations  for  11+  and  for  II”  on 
streamtube  j. 

Since  the  II  pressures  on  matching  edges  of  neighboring  streamtubes  correspond,  as 
shown  in  Figure  A. 5,  the  II  variables  can  be  eliminated  by  taking  the  difference  of  the  II + 
equation  for  cell  i  on  streamline  j  —  l  and  the  II  ~  equation  for  cell  i  on  streamline  j.  The 
resulting  equation  effectively  couples  each  pair  of  streamtubes  through  the  corresponding 
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pressures  on  the  streamline  faces,  which  never  explicitly  appear  as  unknowns  but  whose 
difference  will  approach  zero  as  the  solution  converges.  It  is  interesting  to  observe 
that  the  Newton  system  solves  the  steady  state  equations  by  balancing  the  streamwise 
momentum  equation  and  by  driving  the  streamtube  coupling  pressure  difference,  IIJ  - 
njl  i ,  to  zero. 


A. 6  Energy  Equation 


The  energy  equation  in  the  blade-relative  coordinate  system  rotating  with  angular  ve¬ 
locity  (l  is  a  statement  of  constant  rothalpy  along  the  streamtube, 

I  =  I,  =  I2  where  I  =  h  +  £  -  ^  (A.35) 

with  q  as  the  speed  in  the  blade  relative  frame  and  h  as  the  static  enthalpy  of  the  fluid. 
This  statement  implies  that  the  reference  rothalpy  is  the  stagnation  enthalpy  at  r  =  0. 
The  following  definitions  are  used  for  quantities  in  the  rotating  system. 


The  static  enthalpy  on  the  stream  surface  at  radius  r  is 

*-i-£  +  !2£ 

2  2 


The  stagnation  enthalpy  at  radius  r  is 


ho  —  I  + 


(n  r? 


The  pressure  is  derived  from  the  definition  of  enthalpy  h  =  ~\^  so  that 


P  = 


7-1 


P  I-*T  + 


(nr)2 


The  Mach  number  is  defined  using  a2  =  2E,  which  gives 

2  _  i!  _ 


JVr  = 


(7  -  i)  (I  -  $  +  !2£) 

Finally,  the  temperature  ratio  is  defined  as 


%  _  I  I 

T  h  i  _  al  + 

J.  2  '  2 


(A. 36) 

(A.37) 


(A. 38) 


(A. 39) 


(A. 40) 
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The  usual  definitions  for  stagnation  quantities  can  be  obtained  from  A.40,  but  the 
radius  must  be  considered  in  defining  reference  quantities.  For  example,  the  stagnation 
density  at  r  =  0  in  terms  of  the  density  at  r  with  velocity  q  is 


i 


Po 


I  -  21  + 

V  2  ' 


(»«•)* 


(A.41) 


A. 7  Linearization  of  Equations 


Since  the  flow  equations  and  boundary  conditions  are  non-linear  a  Newton-Raphson 
procedure  is  used  to  solve  the  resulting  numerical  system.  The  S-momentum  equation 
and  reduced  N-momentum  equation  are  linearized  in  terms  of  perturbations  to  the  cell 
geometry  x,  and  cell  face  velocity  q,  density  p  and  pressure  p  in  the  same  manner 
as  in  the  two  ISES  theses  (15,  9).  The  principal  differences  from  their  procedure  stem 
from  additional  terms  due  to  dr/dx  and  db/dx.  This  set  of  unknowns  are  systematically 
reduced  using  the  chain  rule  until  the  momentum  equations  are  expressed  solely  in  terms 
of  the  cell  face  density  changes  6p  and  normal  movement  of  the  streamline  nodes  6n. 
The  details  of  this  process  are  exhaustive  and,  in  the  interests  of  brevity,  are  best  left 
for  a  perusal  of  the  Fortran  coding. 

Several  additional  global  variables  are  also  required  for  perturbations  to  mass  flow, 
inlet  and  exit  angles,  inlet  characteristic,  leading  edge  position,  exit  pressure  and  design 
variables.  These  are  used  to  implement  boundary  conditions  for  analysis  and  design 
applications. 
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Appendix  B 

Development  of  Quasi-3D  Boundary  Layer 


The  integral  boundary  layer  method  developed  for  ISES  is  modified  to  include  effects  for 
streamtube  thickness,  radius,  and  rotation.  These  effects  have  been  included  without 
drastically  altering  the  framework  of  the  basic  method,  as  reviewed  below.  The  details 
of  Drela’s  derivation  of  the  closure  relations  are  not  repeated  as  they  are  unchanged  and 
can  be  best  obtained  from  his  thesis  [9]  and  subsequent  papers  [13]. 


B.l  Background  -  2D  Boundary  Layer 


This  section  presents  the  basic  relations  from  Dreia’s  2D  integral  boundary  layer  method. 
The  boundary  layer  development  is  defined  by  two  principal  equations,  the  classic  von 
Kdrmdn  integral  momentum  equation 


IdS 
9  da 


1 Cj_ 

9  2 


-  (H  +  2  -  M;)  — 

V  'UP 


l  du,: 
da 


(B.l) 


and  the  integral  kinetic  energy  equation,  obtained  by  integrating  the  differential  mo¬ 
mentum  equation  multiplied  by  u. 


J_  d9^_ 

9m  da 


1  duc 
ue  da 


(B.2) 


where  a  is  the  arc-length  along  the  boundary,  uf  is  the  edge  velocity  and  Me  is  the  edge 
Mach  number.  The  shape  parameters  ff,  Hm  and  H"  are  defined  as 


(B.3) 


expressed  using  the  integral  thicknesses 


(B.4) 
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= r  f1  -  -)  -dy  r = r  f1  -  (-)l 

Jo  \  pj  Ue  Jo  \  \U'J  }  peue 

where  q  is  the  wall  normal  coordinate  and  pr  and  u€  are  flow  quantities  at  the  edge  of 

the  EIF.  The  skin  friction  coefficient  C j  and  dissipation  coefficient  Cq  are  defined  by 


Cf  —  _  o  rtt 


-  _L_  /**  — 

peul  Jo  dy 


The  kinetic  energy  equation  B.2  is  not  directly  used,  instead  it  is  combined  with  the 
momentum  equation  to  give  the  so-called  shape  parameter  equation 


j_djr  1  (2 Co  _Cf\_  (2 H2  .  _  H\  1 
H'  da  6\H-  2  J  \  H-  J  uc 


1  due 
uc  da 


In  the  boundary  layer  solution  the  momentum  equation  and  the  shape  parameter 
equation  are  augmented  by  a  third  equation,  depending  on  boundary  layer  type. 

Laminar  A  disturbance  amplification  equation  is  used  to  determine  transition  location. 
The  amplification  equation  is  an  envelope  method,  based  on  correlations  from 
maximum  spatial  amplification  growth  rates  from  the  Orr-Sommerfeld  equation. 

Turbulent  A  dissipation  lag  equation  is  used  to  model  upstream  history  effects.  The 
lag  equation  is  similar  to  Green’s  lag  equation  [10]. 


Laminar  closure  relations,  derived  using  Falkner-Skan  profiles,  provide  correlations  be¬ 
tween  shape  factors,  skin  friction  and  dissipation.  Similarly  for  turbulent  flow,  Swaf¬ 
ford’s  profiles  were  used  to  derive  relations  for  turbulent  closure. 

The  three  equations  characterize  the  development  of  the  boundary  at  any  streamwise 
station  and  are  discretized  along  with  the  Euler  outer  flow  and  solved  by  Newton’s 
method.  Since  the  inviscid  outer  flow  and  the  boundary  layer  are  solved  simultaneously, 
separated  flows  and  transitional  separation  bubbles  can  be  modeled. 
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B.2  Quasi-3D  Boundary  Layer 


The  quasi-3D  equations  are  derived  in  the  same  manner  as  the  2D  equations,  integrating 
the  differential  form  of  the  momentum  equation  but  including  the  effects  of  streamtube 
thickness  b,  radius  r  and  rotation  Cl.  The  momentum  equation  becomes 


1  d9_  _  1  Cj_ 
9da  ~  9  2 


1  due 
uc  da 


Idb 
b  da 


1  dr 
r  da 


(B.8) 


and  the  shape  parameter  equation  is 


1  dH’ 
7F  da 


1  (2  Cp 

6  V  h- 


C/\  /2  H" 

2  )  V  H- 


+  1 


1  dur 
uc  da 


2 H"’\  /Dr\2  Idr 
H~  )  \uc  )  r  da 


(B.9) 


The  integral  thicknesses  used  in  the  shape  parameters  for  B.9  and  B.10  are  defined  as  in 
B.4  and  B.5  but  rdy  is  used  in  place  of  dy  as  the  integration  variable  increment  normal 
to  the  wall.  Similarly  for  the  definition  of  the  dissipation  coefficient  Ci).  The  new  shape 
parameter  that  appears  is  Hp  =  bp/9y  where  the  density  thickness  8P  is  given  by 

8p"  I  rdy  (B-10) 

The  shape  parameter  Hp  is  defined  by  Green  (iOj  as  a  function  of  kinematic  shape  factor 
and  Mach  number 

Hp  ~  (0.185/f*  +  0.15) M*  (B.ll) 


where  Hk  is  the  kinematic  shape  factor  correlated  to  H  and  Mg  by  Whitfield  [21]  as 


L  V1  =  g  -  0.290 Mg 

i+o-u3M« 

y0  \  UgJ  Ur 


(B.  12) 


and  used  by  Drela  in  his  most  of  his  closure  relations  to  account  for  compressible  effects. 


The  additional  effects  of  stream  surface  radius,  streamtube  convergence  and  rotation 
are  included  in  the  right-hand  sides  of  equations  B.9  and  B.10.  The  streamtube  conver¬ 
gence  has  a  relatively  strong  effect  on  the  boundary  layer  as  a  result  of  the  new  term 
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db/da  in  the  integral  momentum  equation.  Rotational  effects  have  a  relatively  small 
influence,  but  affect  both  the  momentum  equation  and  the  shape  parameter  equation. 

The  amplification  equation  and  the  dissipation  lag  equation  are  fortunately  unaf¬ 
fected  by  r,  b  and  Cl  terms.  Additionally,  the  laminar  and  turbulent  closure  relations 
for  shape  parameters,  skin  friction  and  dissipation  were  used  unchanged  from  Drela’s 
derivation  [9],  although  the  current  implementation  of  the  boundary  layer  differs  some¬ 
what  from  that  in  his  thesis  as  a  result  of  later  work  [13]  that  addressed  the  effects 
of  thick  blade  trailing  edges.  Wall  quantities  for  the  boundary  layer,  such  as  ,  A/,., 
T,./T, i  are  defined  using  rothalpy  I  rather  than  enthalpy. 

B.2.1  Inviscid/ Viscous  Coupling 

The  coupling  between  the  boundary  layer  and  the  displacement  surface  is  altered  slightly 
for  the  quasi-3D  treatment.  The  boundary  layer  is  calculated  in  dimensional  quantities 
for  the  integral  thicknesses  using  the  dimensional  surface  length  and  velocity.  The 
presence  of  the  boundary  layer  modifies  the  solid  wall  boundary  condition  by  displacing 
the  stagnation  streamline  outward  from  the  blade  by  the  boundary  layer  thickness.  The 
surface  streamline  is  displaced  An,  in  the  non-dimensional  m'  -  9  coordinates  of  the 
blade  by 

An,  =  (B.13) 

r, 
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Appendix  C 

Cascade  Loss  Calculation 


As  discussed  in  Chapter  2,  losses  through  a  cascade  result  from  inviscid  shock  losses 
and  viscous  dissipation  in  the  blade  boundary  layers.  This  is  manifested  as  a  reduction 
of  the  stagnation  pressure  from  the  equivalent  isentropic  value  at  the  outlet,  assuming 
adiabatic  flow.  The  loss  coefficient  w,  is  defined  as  the  loss  in  stagnation  pressure  in  the 
exit  flow  non-dimensionalized  by  the  inlet  dynamic  pressure. 


Pll;  —  pn-, 

*  i .«  r  n _ 

Pci,  -  Pi 


(C.l) 


where  p, i.,._  is  the  isentropic  (lossless)  exit  stagnation  pressure  and  p„.,  is  the  mass- 
averaged  exit  stagnation  pressure.  The  isentropic  stagnation  pressure  evaluated  at  the 
exit  is 


P",, 


7-1 

7 


hupnr,f 


where  h„  =  I  +  ^  (C.2) 


The  mass-averaged  total  pressure  is  defined  by  integrating  total  pressure  with  respect 
to  mass  flow  at  the  exit. 


u>  = 


Ah  -/' 

J,  I 


dm 


P  cl. 


m 


(C.3) 


P'i,  -  Pi 

The  integration  variable  here  is  dm  =  pubdy ,  the  incremental  mass  flux  at  the  outlet, 
with  b  as  the  streamtube  width.  This  definition  is  used  for  losses  in  rotors  as  well  as 
stators. 


In  MISES  a  zonal  approach  is  used  to  solve  the  viscous  flow  indirectly.  An  equivalent 
inviscid  flow,  EIF,  is  postulated  with  a  displacement  surface  to  represent  the  viscous 
layer.  The  EIF  is  defined  such  that  it  is  locally  irrotational  and  contains  all  the  mass 
flow,  while  the  viscous  effects  are  modeled  by  the  mass  flow  defect.  As  a  result,  the 
inviscid  and  viscous  loss  contributions  must  be  calculated  over  their  respective  domains 
-  inviscid  shock  losses  are  given  by  the  streamtubes  defining  the  EIF,  while  viscous  losses 
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are  calculated  as  a  defect  contribution  from  the  wake  boundary  layer  at  the  cascade  exit 
where  the  flow  has,  at  least  partially,  mixed  out.  This  calculation  takes  the  form  of  C.l 
with  po,  defined  by  the  two  contributions  at  the  exit  plane 


'  dm]  f  dm 

Pu,  =  /  Pu  -T-  -  /  (po«  -  Pu)  -r- 

.J  o  rn  Uo  m 


(C.4) 


C.l  Inviscid  Loss  Calculation 


Inviscid  losses  are  defined  directly  by  the  integral  in  C.3.  The  total  pressure  in  the  EIF 
is  related  to  the  local  density  and  stagnation  enthalpy  by, 


The  inviscid  losses  are  given  by, 


_  7  -  1  , 

Pu  —  - pu  ft u 

7 


r 1  ~ 1  a  /  \drh 

i  fto  (pu,.,,,  -  P"?)  -7- 

u  7 _ 

Pu,  -  Pi 


(C.6) 


where  the  stagnation  density  evaluated  at  the  exit,  using  the  stagnation  enthalpy  h„  = 
I  -f  and  local  density  and  velocity,  is 


(C.7) 


p„,  =  p2  1  - 


The  isentropic  stagnation  density  evaluated  at  the  exit  is 


fh"\~ 

N*. I,  P*' r ,  ,  l  j  I 


(C.8) 


C.2  Viscous  Loss  Calculation 


The  viscous  losses  should  ideally  be  evaluated  far  downstream  where  the  flow  has  mixed 
out  to  a  uniform  condition.  As  an  approximation  to  the  fully  mixed  out  condition  the 
loss  can  be  evaluated  at  the  exit  boundary  where  the  wake  boundary  layer  is  reasonably 
well  developed.  Downstream,  the  wake  shape  factor  H  -*  1,  implying  that  the  velocity 
defect  6u  =  -  u  is  small.  This  will  certainly  be  a  good  approximation  for  unseparated 
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flows  and  will  be  less  valid  as  the  extent  of  flow  separation  increases  to  the  point  where 
the  wake  cannot  fully  develop  in  the  cascade  exit  length. 


The  total  pressure  in  the  boundary  layer  depends  on  p,  u  and  h(). 


(c'9> 

It  is  assumed  that  /i( ,  is  constant  in  the  boundary  layer,  limiting  this  treatment  to  adi¬ 
abatic  flows.  In  the  boundary  layer  p  =  pr,  so  that  C.9  can  be  linearized  for  small 
perturbations  from  the  edge  conditions  (pe,uv)  to  give  an  expression  for  the  total  pres¬ 


sure  defect, 


which  reduces  to, 
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^  6u 


Spn  =  ~p,.U„<5u 
Pr 


(C.10) 


(C.ll) 


The  total  pressure  defect  is  integrated  with  respect  to  mass  flux  across  the  boundary 
layer  to  give, 

.  Pn,  f’tr'i-h 

/  Spudm  =  — -p,.u,b  /  pu(ur-u)rdy  ( C.  1 2 ) 

J  n  p,  Jn 

where  r  and  b  are  the  streamtube  radius  and  width  at  the  exit  plane.  The  last  integral 
is  recognized  as  simply  pru-  9  from  B.4  so  that  the  mass-averaged  loss  in  the  viscous 


yer  is, 
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where  9  is  defined  by  the  integral  BL  and  p, ,  u,  is  defined  by  the  EIF. 


The  viscous  loss  calculation  assumes  near-unity  wake  shape  factors.  The  sensitivity 
of  the  loss  to  the  outlet  length  and  exit  shape  factor  is  shown  in  Table  C.l.  The 
outlet  length  is  measured  from  the  blade  trailing  edge  to  the  exit  plane  and  is  non- 
dimensionalized  by  blade  chord.  The  test  case  was  the  DFVLR  cascade  at  M  =  0.82, 
Rn  =  1.4  x  106  with  no  flow  separation.  The  calculated  loss  levels  are  not  particularly 
sensitive  to  the  exit  length,  and  hence  degree  of  mixing,  even  very  close  to  the  trailing 
edge  for  wake  shape  factors  of  almost  1.7. 
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Table  C.l:  Loss  as  a  function  of  outlet  length 


Outlet  length 

u 

^ in  via 

WUJC 

Hcxit 

0.25 

0.0310 

0.00720 

0.0239 

1.677 

0.50 

0.0309 

0.00715 

0.0238 

1.196 

1.00 

0.0310 

0.00717 

0.0238 

1.055 

2.00 

0.0309 

0.00714 

0.0238 

1.027 

C.3  Loss  at  Uniform  Exit  Condition 


The  combined  inviscid  and  viscous  losses  have  been  calculated  by  mass-averaging  at 
the  exit  plane  of  the  computational  domain.  The  following  conservation  scheme  can  be 
used  to  extrapolate  the  mass-averaged  flow  quantities  to  a  hypothetical  fully  mixed  out 
flow  downstream. 


The  mass  flow,  total  momentum  components  and  enthalpy  are  known  at  the  mass¬ 
averaging  plane,  referred  to  using  ()j  subscripts.  Conditions  at  the  fully  mixed  condition 
are  denoted  by  ()  t.  Conservation  of  mass,  momentum  and  energy  at  the  two  planes  give 

P\q.r,S]b;\  =  piq -  m 


rhqtn  — 


mqT,  +  p.\h  b  \ 


mq, ,  —  P j  bj 


7  P  i  ?  ;i  _  7  P  :  |  } 

7-1 P.)  2  7-1  h  2 


(C- 14) 


where  5  =  r  yptlch  is  the  physical  cascade  pitch  at  r  and  b  is  the  streamtube  width.  The 
fixed  quantities  are  m  and  h»,  which  are  constant  at  all  stations,  and  the  radius  1*3  and 
width  6  )  at  the  mixed-out  plane.  The  mass-averaged  quantities  &  and  p>,  and  the  exit 
angle  0-<  are  known  from  the  flow  calculation.  The  mass-averaged  velocity  is  given  from 
the  momentum  integral  over  the  inviscid  outer  flow  and  a  momentum  defect  integral 
over  the  boundary  layer  at  the  exit. 
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\.Jo  m>  J  K/f  in  Iv  o  1  m j 
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[Jo  ^  .  E/F  771 


(C.15) 


A  simple  mass-averaging  in  the  EIF  suffices  for  defining  p2 .  The  exit  angle  (32  is  given 
by  a  global  unknown,  either  specified  or  part  of  the  solution.  The  x  and  y  components 
of  qt  are  defined  using  the  exit  angle  (32.  This  system  defines  four  equations  in  four 
unknowns,  so  the  solution  is  straightforward. 

?!/!  =  Qy? 


Pi  = 


h  faVj  j 

Pi  a ■>  b2  +  m{qI7  -  q „) 


/?;)  =  arctan 


where  y,.n  is  defined  by  the  solution  to  the  quadratic 


f  7  +  1  1  ,  T  7  +  Pihbi]  q;n 


9r"  [2(1  -7)J  +  qr*  1 7  -  1  m 


The  loss  in  the  mixed-out  flow  is  defined  using  the  mixed  out  total  pressure 


Ph  =  Pi 


hi,  - 


(C- 16) 


+  22--A,,  =  0  (C.17) 


(C.  18) 


in  place  of  p,,,  in  the  loss  equation  C.l.  The  loss  defined  using  this  approach  normally 
differs  only  slightly  from  the  mass-averaged  result. 
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C.4  Loss  Sensitivity  Calculation 


One  of  the  useful  byproducts  of  using  a  Newton  solution  procedure  is  the  Jacobian 
matrix  J  that  is  generated  and  inverted  for  each  Newton  step.  The  Jacobian  defines 
the  sensitivity  of  the  solution  with  respect  to  the  grid,  density  and  global  unknowns. 
This  information  can  also  be  used  to  to  calculate  the  sensitivity  of  any  flow  dependent 
quantity,  such  as  the  loss  coefficient  w  =  /(n,p),  to  the  global  variables.  The  global 
variables  can  include  the  global  parameters  such  as  inlet  angle  or  Mach  number,  exit 
pressure  or  geometric  design  modes  that  perturb  all,  or  part  of,  the  blade  surface.  This 
opens  the  door  to  a  host  of  interesting  possibilities  for  design  optimization. 


The  linear  system  for  the  Newton  procedure  has  the  form 


r  ,  1 

<5n 

’ 

J'  h 

-R' 

1 

< 

bp 

.  =  < 

— 

G,,.,,  |  G, 

/g. 

,-R*. 

(C.19) 


where  the  unknown  vector  contains  the  Eulef  flow  field  unknowns  bn,  bp,  and  the  global 
unknowns  b g.  The  Jacobian  matrix  J  is  composed  of  flow  equation  rows  (J'  +  h)  and 
global  equation  rows  (G„.#)  4-  Gq).  The  vector  R'  is  the  residual  for  the  flow  equations, 
R,;  is  the  residual  for  the  global  constraint  equations. 


The  flow  equations  are  tightly  banded  with  respect  to  the  flow  variables  but  not 
with  respect  to  the  global  variables.  This  linear  system  is  most  efficiently  solved  by 
rearranging  the  columns  corresponding  to  the  global  variables  to  the  right-hand  side  to 
get 


J' 

< 

M 

i  i 

»  =  < 

H 

'  +  6g  i  - 

1  J 

>  »  • 

*  +  Sgm  < 

[-  h„,  ’ 

W 

i  J 

l  j 

l  1 

(C.20) 


where  the  right-hand  side  consists  the  flow  variable  residuals,  R',  and  the  h  global 
variable  columns.  The  global  equations  define  a  second,  much  smaller  system 


Gn./; 


+  GgS  g  -  -Rg 


(C.21) 


141 


The  system  of  flow  equations  C.20  is  solved  by  elimination  to  obtain  a  partially 
determined  solution. 


f Sn ] 

f(fn) 

L 

i 

>  =  <  >  +  Sg\  < 

- 

h, 

'  •  ■  •  +  5gm  < 

h„ 

UpJ 

J 

J 

(C.22) 


where  the  vectors  Sn  and  Sp  are  the  flow  solution  without  the  effect  of  the  global 
variables  and  h*  are  the  vectors  of  sensitivities  of  flow  variables  to  the  global  unknowns. 


ht  = 


On 

%k 

Op 

■  Ofik  ■ 


(C.23) 


The  solution  is  obtained  by  solving  the  small  global  system  C.21  using  the  partially 
determined  solution  Sh  and  tip. 

[G,,„  (h,  ■••hrn)  4-  Gv]  Sg  =  -R,  -  G„.„  |^J  ( C . 2 4 ) 

The  flow  unknowns  Sn  and  Sp  are  determined  from  the  partial  solution  using  the  global 
variables  S g  and  the  h  vectors. 


The  components  of  the  h  vectors  can  be  used  to  calculate  the  sensitivity  of  the  loss 
coefficient  u(n.p)  to  a  global  variable  change.  The  sensitivity  is  obtained  by  linearizing 
the  loss  w.r.t.  the  flow  unknowns,  generating  the  vectors  du >/dp,  and  du  dn,.  The 
sensitivity  to  the  global  unknowns  is  given  by 

du  _  du  dn,  du  dp,  (C  2") 

dgk  ,  dn,  dgk  "  dp,  dgk 

Note  that  the  sensitivity  of  any  flow  dependent  quantity,  such  as  mass  flow  or  turning, 
could  also  be  linearized  in  this  manner. 


Unfortunately  the  global  unknowns  are  not,  in  general,  orthogonal  with  respect  to 
one  another  so  that  the  effect  of  a  perturbation  in  a  single  global  variable  on  the  set  of 
global  unknowns  must  be  accounted  for.  Using  the  small  global  matrix  system  C.24, 
each  global  variable  Sg}  is  in  turn  perturbed  by  replacing  its  regular  constraint  equation 
with  Sgj  =  1.  The  right  hand  sides  of  the  remaining  rows  are  cleared  and  the  global 
system  is  solved  for  the  resulting  vector  of  perturbations  to  the  global  unknowns  Sgk, 


142 


This  is  the  response  of  the  global  unknowns  to  a  single  unit  global  variable  change  6gj. 
The  linearized  loss  sensitivity  to  a  global  variable  change  is  defined  by  the  inner  product 

Sw  8lj  6gk 

&9j  ~  j  d9k  Sg} 

This  procedure  is  repeated  for  each  desired  sensitivity. 

The  linearized  loss  sensitivity  strictly  only  applies  in  the  vicinity  of  the  point  at  which 
it  was  generated,  so  that  some  caution  must  be  used  in  interpreting  this  information  for 
design  changes  using  the  modal  redesign  feature.  Experience  shows  that  the  sensitivity 
can  be  used  as  a  guide  to  a  steepest  descent  optimization,  provided  sufficiently  small 
changes  are  introduced.  The  sensitivity  information  is  only  valid  near  convergence  of 
the  Newton  system  so  that  the  best  strategy  is  to  introduce  a  modest  change,  guided  by 
the  calculated  optimization  direction,  and  run  for  several  Newton  cycles  and  reevaluate. 
One  of  the  difficulties  associated  with  this  design  strategy  is  that  it  tends  to  lead  to 
“point  optimizations”  rather  than  robust  designs. 
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Appendix  D 

Matrices  and  Linear  System  Solution 


The  theory  and  numerical  implementation  of  the  quasi- 3D  flow  solver  has  been  described 
in  previous  chapters  and  appendices.  This  appendix  includes  details  of  the  Jacobian 
matrices  generated  by  the  Newton  system  and  the  approach  used  for  solving  the  linear 
system.  Modifications  to  the  basic  scheme  for  multiple  blades  and  offset-periodic  grids 
are  discussed. 


D.l  Newton  System 


As  introduced  in  chapter  3  the  inviscid  flow  is  solved  with  the  steady-state  Euler  equa¬ 
tions  discretized  on  a  streamline  grid,  see  Figure  D.l.  A  Newton- Raphson  solution 


Figure  D.l:  Newton  variables  for  discrete  Euler  equations. 

technique  is  used  to  iteratively  solve  the  resulting  non-linear  equations.  As  a  result 
of  the  streamline  approach,  only  streamwise  convection  terms  are  present,  effectively 
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reducing  the  system  of  equations  to  a  streamwise,  or  S-momentum  equation  at  conser¬ 
vation  cell  centers  in  each  streamtube  and  a  normal  momentum  (reduced  N-momentum) 
equation  at  nodes  on  each  streamline.  The  Jacobian  system  is  constructed  by  taking 
variations  of  these  equations  with  respect  to  each  of  flow  and  geometric  variables.  These 
variations  are  compactly  expressed  in  terms  of  one  unknown  Newton  variable  for  each 
equation,  a  perturbation  to  the  streamtube  cell  density  6p  and  a  perturbation  to  the 
streamline  normal  node  displacement  6n,  see  Figure  D.l. 

A  streamline  grid  is  defined  with  streamwise  index  i  and  quasi-normal  index  j  and  is 
made  up  of  (I  nodes  on  each  streamline)  x  J  streamlines.  The  application  of  the  discrete 
equations  defines  I  x  (J  -  1)  S-momentum  equations  and  I  x  J  reduced  N-momentum 
equations,  which  are  efficiently  arranged  as  I  block  rows  of  ( 2  J  - 1 )  equations  in  each 
block.  The  i  =  1  and  i  =  I  block  rows  are  used  to  apply  the  inlet  and  exit  boundary 
conditions.  Within  each  block,  the  j  =  1  and  j  =  J  rows  are  used  apply  periodicity 
conditions  to  the  top  and  bottom  streamlines. 

The  domain  of  dependence  ‘molecule’  for  the  equations  and  variables  on  interior 
streamlines  is  shown  in  Figure  D.2.  The  equations  span  four  streamwise  (i)  nodes  and 


Figure  D.2:  Domain  of  dependence  for  S-momentum  and  N-momentum  equations. 

one  streamtube  for  the  S-mom.  equation,  two  streamtubes  for  the  N-mom.  equation. 
The  Z,  block  dependence  from  nodes  at  i  -  2  results  solely  from  the  influence  of  the 
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upstream  Mach  number  in  the  artificial  dissipation  and  the  influence  of  ue  in  the  BL 
equations.  An  additional  V,  block,  similar  to  the  Z,  block  but  defined  at  node  line  i  -  3, 
is  added  if  second-order  dissipation  is  used. 

Additional  unknowns  are  often  required  that  do  not  fit  within  the  compact  block 
structure  of  the  interior  equations.  These  global  variables  and  associated  global  equa¬ 
tions  are  specially  treated  in  the  solution  process. 


D.1.1  Jacobian  Matrix 


At  each  Newton  iteration  a  linear  system  must  be  solved  for  the  changes  in  density  and 
streamline  node  displacement.  The  basic  form  of  the  block  matrix  system  is 


A, 

c, 

B. 

A. 

c. 

Z:, 

B, 

A:) 

c 

v, 

z, 

B, 

A, 

i  H, 
l  H; 
I  H. 

.  H, 


-it, 
-it. 
-ft  i 
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(D.l) 
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-R 

-R, 

-R, 


The  linear  system  is  arranged  in  I  blocks  of  dimension  2J-1,  where  each  block  has  an 
unknown  vector  6,  =  (5nt>i_y,  R,  is  the  Newton  residual  for  the  S-mom. 

and  N-mom.  equations.  The  vector  of  global  unknowns,  S9,  is  associated  with  the  global 
equations,  made  up  of  G,  entries  for  the  6p  and  in  unknowns,  and  Gg  entries  for  the 
global  unknowns.  The  H,  blocks  represent  the  effect  of  the  global  unknowns  on  the 
S-mom.  and  N-mom.  equations. 
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D.1.2  Block  Structure 


The  component  blocks  for  the  interior  equations  in  D.l  are  shown  in  Figures  D.3  and  D.4. 
Periodic  conditions  for  pressure  and  geometry  are  enforced  by  the  first  and  last  rows  in 
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Figure  D.3:  Structure  of  A,  and  C,  blocks. 
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Figure  D.4:  Structure  of  Z,  and  B,  blocks. 


the  upper  sub-block  corresponding  to  Sn  unknowns.  Geometric  continuity  is  imposed 
by  the  identity  condition  in  the  A,  block  only.  Note  that  the  periodic  pressure  condition 
introduces  two  node  displacement  entries  and  one  density  entry  for  the  corresponding 
periodic  streamline  (i.e.  streamline  1  is  coupled  to  streamline  J,  involving  streamtubes 
1  and  J-l).  The  geometric  continuity  condition  introduces  only  one  ‘off  diagonal’  entry. 
For  a  solid  surface  BC,  such  as  over  a  blade  surface,  the  periodic  pressure  and  geometric 
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rows  we  cleared  and  are  replaced  by  1  on  the  diagonal. 


The  blocks  corresponding  to  the  inlet  conditions  are  defined  in  Figure  D.5.  At  the 
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Figure  D.5:  Structure  of  inlet  blocks. 

inlet,  only  the  A*  and  C|  blocks  are  present.  These  implement  linearized  slope  or 
inlet  characteristic  conditions  in  the  upper  sub-blocks  and  linearized  stagnation  density- 
conditions  in  the  lower  sub-blocks. 


The  blocks  for  the  last  line  of  grid  nodes  at  i  -  I  implement  the  exit  conditions 
and  are  defined  in  Figure  D.6.  At  the  exit,  only  blocks  V/  through  A \  are  present. 


fn/.i  —  j  'lu 


mnn  vars. 


outlet 
flow  angle 
conditions 
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equations 
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A  /  block 


Figure  D.6:  Structure  of  exit  blocks. 


The  upper  sub-blocks  implement  the  linearized  exit  slope  boundary  conditions.  Density 
variables  are  not  defined  at  i  =  I  so  the  lower  sub-blocks  contain  dummy  equations. 
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D.1.3  Matrix  Blocks  with  Boundary  Layer 


For  viscous  flows  each  surface  boundary  layer  is  represented  using  integral  boundary 
layer  equations  which  are  solved  simultaneously  with  the  inviscid  outer  flow.  The  details 
of  the  boundary  layer  discretization  are  best  covered  in  Reference  [9j.  Modifications  to 
the  ISES  methodology  for  quasi-3D  flow  are  discussed  in  Appendix  B. 

The  boundary  layer  adds  three  equations  and  three  unknowns  to  the  Newton  system 
for  every  node  on  each  surface  streamline.  The  new  unknowns  correspond  to  6C!r1' ,  69 
and  6(5')  for  the  wall  streamline.  The  boundary  layer  solution  is  driven  by  the  edge 
velocity  and  density,  defined  in  terms  of  the  unknowns  for  the  streamtube  at  the  wall. 
The  domain  of  dependence  ‘molecule’  for  the  viscous  variables  is  shown  in  Figure  D.7. 


Figure  D.7:  Domain  of  dependence  of  boundary  layer  unknowns,  from  Drela. 
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An  A,  block  with  boundary  layer  coupling  conditions  takes  the  form 
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On  the  periodic  inlet  streamline,  where  no  boundary  layer  is  present,  the  boundary  layer 
entries  are  cleared  and  the  diagonal  entries  set  to  one. 

Wake  boundary  layers  are  handled  somewhat  differently  -  the  upper  and  lower  sur¬ 
face  boundary  layer  unknowns  are  merged  at  the  trailing  edge  and  the  wake  is  treated 
as  a  single  boundary  layer.  Dummy  equations  are  used  as  placeholders  in  the  unused 
BL  rows.  The  pressure  and  geometric  continuity  conditions,  from  D.3,  are  imposed  with 
the  geometric  continuity  condition  modified  to  include  the  wake  thickness  between  the 
periodic  streamlines. 
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D.2  Solution  of  Linear  System 


The  linear  system  D.l  is  solved  by  treating  it  as  two  lineaT  problems,  a  banded  block 
system  and  a  smaller,  full  linear  system  for  global  unknowns.  The  original  system  can 
be  rewritten  as 
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where  D.2  can  be  efficiently  solved  using  a  specialized  block  Gauss  elimination  algo¬ 
rithm,  as  described  in  Giles  (15].  The  banded  system  is  solved  for  both  the  right-hand 
side  residual  vector  R  and  the  additional  right-hand  columns  of  global  influences  on 
flow  equations,  H.  The  resulting  partially-determined  solution  is  given  by  the  vectors 
6  and  h  such  that  the  complete  solution  will  be 

6,  =  6,  -  h ,6g  (D.4) 

when  Sg  is  determined.  The  small  linear  system  is  obtained  by  substituting  D.4  into 
the  global  equation  D.3.  This  gives  a  linear  system  in  the  global  variables,  typically  of 
order  2  or  3  but  up  to  ~  12  for  modal  design  applications. 

[Gg  -  £Gth,l*g  =  -Rg  -  £<M,  (D.5) 

1  I 

Substitution  of  Sg  into  D.4  gives  the  fully-determined  solution. 
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D.2.1  Solution  Method  for  Choked  Flow 


The  banded  matrix  system  in  D.2  becomes  increasingly  ill-conditioned  as  the  flow  in 
the  passage  becomes  choked,  requiring  modifications  to  the  boundary  conditions  and 
solution  procedure.  As  shown  by  Giles  [15],  the  well-posed  boundary  condition  for 
choked  flow  specifies  the  gradient  of  inlet  stagnation  density  (which  is  equivalent  to 
specifying  the  inlet  vorticity)  and  the  exit  pressure.  One  of  the  remaining  inlet  boundary 
conditions  must  be  dropped  —  either  the  mass  flow  or  the  inlet  stagnation  density.  Since 
mass  flow  has  been  implicitly  written  into  the  momentum  equations  the  inlet  boundary 
condition  on  stagnation  density  is  dropped. 

The  discrete  form  for  the  new  inlet  condition  specifies  the  gradient  of  inlet  stagna¬ 
tion  density,  usually  set  to  zero.  This  is  implemented  by  simply  subtracting  the  inlet 
stagnation  density  conditions,  equation  3.4,  of  neighboring  streamtubes,  i.e. 

/>",.>  -  >*,.  =  0  (D-6) 

This  produces  7-1  inlet  density  equations,  the  missing  equation  consists  of  a  mass- 
averaged  exit  pressure  condition 


;=i 


To  maintain  the  banded  structure  of  the  matrix  this  equation  must  be  moved  to  the  7-1 
block.  A  space  is  cleared  in  the  7-1  block  by  shifting  one  of  the  S-momentum  equations, 
at  j  =  jnhi/ti  UP  into  the  preceding  block  row,  see  Figure  D.8.  The  corresponding 
equation  in  the  7-2  row  is  shifted  up  into  the  7-3  row,  etc.  until  the  slot  in  the 
7=1  block  is  filled.  Observing  the  A  and  C  blocks  in  Figure  D.3,  note  that  only 
three  new  entries  in  the  shifted  row  fall  outside  the  existing  diagonal  block  structure  in 
V  — ►  C.  These  entries  are  treated  by  the  solver  essentially  as  three  extra  right-hand 
sides,  incurring  a  10%  penalty  in  solution  cost. 

Unfortunately  the  choked  solution  determines  the  inlet  stagnation  density  for  a  given 
mass  flow  and  exit  pressure,  rather  than  the  mass  flow  for  a  given  inlet  stagnation  density 
and  exit  pressure.  To  re-normalize  the  inlet  stagnation  density  to  the  reference  value 
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D.3  Multiple  Blades 


The  application  to  multiple  blades  adds  additional  complexity,  but  no  profound  changes 
to  the  basic  system.  The  effect  of  adding  one  or  more  additional  blades  is  handled  by 
treating  the  grid  as  a  single  large  grid  domain  with  embedded  stagnation  streamlines, 
which  in  effect  separate  it  into  N  passages,  where  N  is  the  total  number  of  blades.  Each 
stagnation  streamline  corresponds  to  two  streamlines  that  define  a  null  streamtube  and 
whose  points  are  co-located  up  to  the  leading  edge  of  the  blade,  see  Figure  D.9. 

Str*amtubes 


Streamlines 

Figure  D.9:  Layout  of  streamtubes  and  streamlines  in  multiple  blade  grid 


The  basic  arrangement  of  nodes,  streamlines,  equations  and  variables  is  unchanged, 
except  that  dummy  equations  replace  the  flow  equations  defining  the  densities  in  each 
null  streamtube.  Additional  periodicity  conditions  must  be  enforced  for  pressure  and 
geometric  continuity  across  the  collapsed  streamtube.  The  domain  of  dependence 
‘molecule’  for  the  equations  and  variables  across  the  null  streamtube  is  shown  in  Fig¬ 
ure  D.10.  The  S-momentum  equation  operates  on  a  single  streamtube  and  is  unchanged. 
Note  that  the  N-mom.  equation  spans  three  streamtubes  in  j. 
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Figure  D.10:  Domain  of  dependence  for  reduced  N-momentum  equation  across  stagna¬ 
tion  streamline. 

D.3.1  Block  Structure  for  Multiple  blades 


The  basic  structure  for  the  Jacobian  matrix  blocks  is  modified  to  establish  periodic  con¬ 
ditions  and  dummy  equations  for  the  interior  stagnation  streamlines.  The  new  structure 
for  the  A,  matrix  block  is 


XX  XX 

XXX 

XXX 

X  X  X  X 

XXX 

XXX 

I  -  \ 


i  —  y-i 

•  • 
•  • 

•  • 

•  • 

•  • 

•  • 


*—  side  1,  periodic  pressure 

•—  side  3,  periodic  geometry 
«—  side  4,  periodic  pressure 

<—  side  2,  periodic  geometry 

<—  null  streamtube 

|  BL  rows 


This  block  structure  corresponds  to  the  grid  in  Figure  D.9  where  streamtube  4  is  the 

s 

null  streamtube,  nestled  around  the  intermediate  blade.  Note  that  this  is  a  A,  block 
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for  a  station  without  boundary  layers.  Matrix  blocks  V*  through  C,  are  defined  anal¬ 
ogously,  without  the  geometric  continuity  conditions.  The  periodic  pressure  condition 
for  the  intermediate  blade  is  effectively  the  same  as  that  used  for  the  j  =  1  and  j  =  J 
streamlines,  but  it  applies  across  the  null  streamtube. 


A  case  with  N  blades  uses  6N  boundary  layer  equations,  increasing  size  of  the  block 
matrix.  This  increases  the  solution  cost  -  a  two  blade  solution  is  approximately  40% 
more  work  than  a  single  blade  case  (for  J  —  20  streamlines). 
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D.3.2  Multiple  Blade  Choked  Flow  Solution 


The  choked  flow  solution  procedure  is  modified  for  multiple  blades  by  independently 
treating  each  passage.  A  shifted  stream  tube  and  pexu  condition  D.7  are  used  for  each 
passage.  A  mass  flow  or  differential  mass  flow  DOF  is  added  for  each  passage  and 
constrained  by  inlet  stagnation  density  conditions  D.8.  The  Kutta  condition  for  the 
intermediate  blade,  usually  constraining  the  differential  mass  flow,  must  be  dropped  as 
an  explicit  constraint  and  applied  implicitly  in  the  banded  matrix  by  equating  the  upper 
and  lower  II  pressures  at  the  trailing  edge. 
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D.4  Offset-Periodic  Matrix 


A  special  grid  topology  using  an  offset-periodic  inlet  grid  was  introduced  in  Chapter  3 
to  reduce  grid  skew  problems  for  supersonic  Mach  numbers.  The  inlet  portion  of  the 
grid  was  constructed  with  quasi-normal  lines  roughly  normal  to  the  inflow,  resulting  in 
a  ‘stairstep’  inlet  boundary,  see  Figure  D.ll.  Periodic  conditions  must  now  be  enforced 
on  matching  points  on  the  j  =  1  and  j  =  J  streamlines  which  are  offset  in  streamwise 
index  by  a  factor,  i„jj  (note  that  i0jj  =  6  in  the  figure).  The  application  of  periodic 
pressure  and  geometry  conditions  leads  to  three  pressure  terms  and  one  geometry  term 
moving  to  off-diagonal  blocks  offset  by  i0jj. 


Figure  D.ll:  Offset-periodic  inlet  geometry  and  indexing. 
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The  block  structure  of  the  matrix  for  an  offset-periodic  grid  has  the  form 
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The  offset  for  the  matrix  above  corresponds  to  the  inlet  grid  in  Figure  D.ll.  The  global 
equation  rows  and  global  column  vectors  on  the  right-hand  side  are  not  shown  but  have 
the  same  form  as  D.2.  The  shifting  affects  only  the  first  part  of  the  matrix,  up  to  the 
leading  edge  of  the  main  blade,  where  the  offset  stops  (at  t  =  7  in  the  figure).  The 
remainder  of  the  grid,  including  the  exit,  uses  regular  unshifted  periodic  conditions. 

The  offset-periodic  pressure  and  geometry  terms  define  two  new  matrix  blocks,  both 
extremely  sparse.  The  three  pressure  terms  have  been  shifted  left  by  i0jj  into  the  lower 

*  /  3  / 

triangle  and  form  the  V,  through  C,  blocks,  each  corresponding  to  terms  shifted  out  of 
one  of  the  principal  blocks  near  the  diagonal.  The  geometric  periodicity  term  has  been 
shifted  right  by  i0/j  into  the  upper  triangle.  This  shifting  configuration  entails  the  least 
‘fill-in’  during  the  elimination  process. 
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D.4.1  Offset  Block  Structure 


The  blocks  V,  through  C,-  in  the  matrix  lower  triangle  contain  the  shifted  pressure 
periodicity  terms.  These  are  shown  on  the  left  side  of  Figure  D.12.  The  blocks  V, 
through  C,  near  the  diagonal  are  modified  as  well  -  all  the  terms  corresponding  to 
those  in  the  offset  blocks  are  cleared.  The  C,  block  differs  slightly  because  it  contains 
no  6p  terms.  The  block  U,,  consisting  of  a  single  non-zero  term  for  the  geometry 


l  —  J 


bpi-m,  1  —  7-1 


x  x 


+  +  I —•  J-\ 


*—  offset  pressure 


*—  offset 


geometry 


V',,  Z',,  B',,  A',,  C',  blocks 


U',  block 


Figure  D.12:  Structure  of  lower  and  upper  triangle  blocks  for  offset-periodic  matrix 


matching,  is  shifted  into  the  upper  triangle  of  the  matrix,  This  block  is  shown  in  the 
right  side  of  Figure  D.12.  The  corresponding  term  from  A,  is  cleared. 


D.4.2  Offset  Matrix  Solution 

The  offset-periodic  matrix  is  solved  using  a  modification  of  the  block  Gauss-elimination 
solver  from  ISES  [15].  The  elimination  process  itself  is  straight-forward,  but  becomes 
quite  complex,  particularly  in  the  offset  region,  due  to  the  need  for  data  structures  to 
minimize  storage  of  the  upper  triangular  terms.  The  lower  triangle  blocks  V  -*  and 

=  3  -  //  =  // 

V  -*  B  are  eliminated  using  the  factored  upper  triangle  blocks  C  and  U  .  The  factored 
entries  for  U  are  stored  in  an  array  of  column  vectors  whose  columns  fill  during  the 

a  /  a 

elimination,  from  the  original  block  U,  down  to  the  C,+In//  block.  The  general  form  of 
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the  factored  matrix  with  upper  triangle  offset  columns  is 
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The  extra  cost  associated  with  the  offset  solver  depends  on  the  total  number  of  offset 
blocks,  each  of  which  adds  a  portion  of  the  computational  cost  of  an  extra  right-hand 
side  (~3%  per  RHS).  The  fraction  of  the  extra  cost  compared  to  a  RHS  is  dependent 
on  the  relative  size  of  the  offset  region.  In  practice,  an  offset  solution  is  2-3  times  slower 
than  a  comparably  sized  standard  solution. 
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